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ABSTRACT 


The  development  of  a  unified  picture  of  short-pulse  high-intensity  multiphoton  processes, 
embracing  atoms,  molecules,  and  solids,  appears  possible  through  the  study  of  clusters.  Of 
particular  significance  are  possible  intra-cluster  processes  that  can  influence  the  mechanism  of 
ionization  and  lead  to  the  production  of  inner-shell  vacancies.  Inner-shell  excitation  leading  to 
prompt  X-ray  emission  is  specifically  considered  and  the  treatment  leads  to  the  definition  of  a 
critical  cluster  size  n^,  representing  the  achievement  of  maximal  X-ray  emission  from  the 
ensemble.  These  results  suggest  the  possibility  of  designing  a  new  class  of  molecular  materials 
optimized  for  the  efficient  production  and  amplification  of  X-rays.  Recent  experiments, 
specifically  designed  to  evaluate  certain  essential  predictions  of  the  cluster  picture,  have 
furnished  considerable  supporting  evidence  for  the  main  features  of  the  proposed  model.  These 
studies,  which  involved  Kr  and  Xe,  (1)  established  the  crucial  role  of  cluster  formation  in  the 
generation  of  soft  X-ray  (~  80  -  150  eV)  emission,  (2)  verified  the  scaling  of  this  new 
phenomenon  in  to  the  kilovolt  (~  800  -  2100  eV)  region,  (3)  provided  important  insights 
concerning  the  dynamical  character  of  the  strong-field  coupling  to  the  clusters,  (4)  demonstrated 
the  production  of  "hollow"  atoms,  and  (5)  combined  the  cluster  excitation  with  a  recently 
discovered  mode  of  channeled  propagation.  These  results  strongly  indicate  that  the  predicted 
amplification  can  be  achieved. 
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I.  INTRODUCTION 

'  Our  recent  work  has  been  aimed  at  the  study  of  the  properties  of  X-ray  emission  from 
multiphoton  excited  molecules  and  clusters.  As  discussed  in  this  report,  the  basic  scaling 
relations  stemming  from  a  general  analysis  [1]  of  prompt  X-ray  emission  from  ensembles  of 
closely  grouped  atoms  indicate  that  clusters  of  heavy  atoms  and/or  suitably  designed  molecules 
containing  heavy  atoms  have  unusually  favorable  properties  for  the  production  and  amplification 
of  X-rays  in  the  kilovolt  range.  The  scaling  laws  found  indicate  a  very  strong  dependence  on 
certain  atomic  and  molecular  properties  which  favors  heavy  materials  in  a  potent  and  obvious 
way.  In  this  analysis,  explicit  contact  is  made  with  six  independent  experimental  observations 
involving  three  materials  (Ar,  Kr,  and  Xe),  all  either  unexplained  or  anomalous  findings 
concerning  the  generation  of  short  wavelength  (X  <  200  A)  radiation.  Consistent  agreement  of 
these  data  with  the  analysis  is  found  in  every  case.  It  is  considered  significant  that  this  synthesis 
can  be  achieved  with  a  treatment  that  only  involves  known  atomic  and  molecular  properties 
without  the  incorporation  of  additional  free  parameters. 

Four  new  experimental  results  [2,3]  concerning  multiphoton  induced  emissions  (1)  from 
Xe  in  the  145  -  152  A  range,  (2)  from  Kr  in  the  72  -  140  A  region,  (3)  from  Xe  in  the  11  -  17 
A  range,  and  (4)  from  Kr  in  the  7  A  region,  arising  from  experiments  explicitly  designed  to  test 
the  cluster  model,  provide  strong  confirming  evidence  for  the  physical  picture  presented.  These 
results  support  the  concept  of  a  molecular  design  of  specific  systems  for  the  production  and 
zunplification  of  X-rays. 

The  prospects  of  X-ray  amplification  are  enormously  enhanced  by  knowledge  stemming 
from  previous  studies  [4,5]  of  confined  propagation  in  which  channel  formation  has  been 
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observed  [4],  This  physical  mechanism  enables  the  spatial  distribution  of  the  excited  material 
to  be  efficiently  produced  in  a  geometry  suitable  for  amplification  [6-8]. 

The  overall  goal  of  our  program  is  the  union  of  these  two  physical  processes,  namely, 
(1)  the  cluster  excitation  and  (2)  the  channeled  propagation,  in  order  to  produce  X-ray  sources 
of  exceptional  brightness  for  3-dimensional  microimaging. 

n.  DISCUSSION  or  research 

A.  Multiphoton  X-Ray  Generation  from  Molecules  and  Clusters 

■#  ■ 

A  unified  picture  of  short-pulse  high-intensity  mi^^^oton  processes,  embracing  atoms, 
molecu^^^d  solids  has  been  proposed  in  relatioft^to  llw  's^y  of  the  behavior  of  ekistd^  [1]. 
This  model  is  described  in  Appendix  A.  Specific  experimental  teits  of  the  cluster  concept  have 
also  been  performed,  the  most  informative  of  which  were  conductetfin  the  kilovolt  region  and 
involved  the  observaticm  of  Kr  L-shell  and  Xe  M-shell  emis^on^gliesejindings  are  presented 
in  Appendices  B  and  C.  The  production  of  "hollow"  atorps  was  also  defflonstrated  with  the 
observation  of  Xe(L)  emission  in  the  4  -  5  keV  range,  a  result  presented  in  Appendix  D. 

Important  scaling  relationships  have  been  established  by  these  studies,  particularly 
concerning  the  energy  transfer  rate.  These  results  are  given  in  Appendix  E.  The  main 
conclusion  is  that  energy  deposition  rates  exceeding  ~  1  W/atom  are  feasible  in  appropriately 
designed  clusters  or  molecules  incorporating  heavy  atoms. 

B.  Channeled  Propagation 

The  spatial  control  of  the  energy  deposited  for  excitation  plays  a  crucial  role  in 
establishing  the  conditions  for  X-ray  amplification.  The  existence  of  a  confined  mode  of 
propagation,  which  occurs  under  conditions  identical  with  those  suitable  for  excitation  of  the 
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clusters,  is  now  supported  by  both  theoretical  [6,9]  analysis  and  experimental  [10]  evidence,  as 
described  in  Appendices  F  and  G,  respectively.  This  phenomenon  has  an  enormous  influence 
on  the  scaling  relationships  [8]  governing  amplification,  the  analysis  of  which  is  presented  in 
Appendix  H.  The  main  conclusion  is  that  -  1  J  of  excitation  energy  is  sufficient  to  produce 
strong  amplification  at  quantum  energies  well  into  the  kilovolt  range, 
m.  RESEARCH  FINDINGS 

The  goal  of  our  work  is  the  union  of  the  material  excitation,  discussed  in  Section  II. A, 
with  the  confined  propagation,  described  in  Section  H.B.  This  step  would  bring  into  conjunction 
the  two  basic  requirements  for  efficient  X-ray  amplification,  namely,  (1)  a  highly  excited 
nonequilibrium  medium  and  (2)  a  spatial  distribution  configured  for  directional  amplification. 

Preliminary  experimental  findings  obtained  by  imaging  the  X-ray  emitting  zone  with  a 
pin-hole  camera  indicate  that  this  conjunction  is  possible.  The  results  involve  Xe  M-shell 
emission  (-  1  keV)  produced  under  the  conditions  described  in  Appendix  C.  The  pin-hole 
camera  was  positioned  so  that  a  transverse  view  of  the  excited  region  was  observed.  The 
camera  was  equipped  with  either  a  100  or  5  /im  pin-hole  and  a  filter  that  eliminated  radiation 
with  a  quantum  energy  less  than  ~  500  eV.  The  signal  was  recorded  on  a  two-dimension  array 
and  stored  in  a  microprocessor.  The  exposures  were  sufficiently  strong  that  single-pulse  images 
were  recorded. 

The  image  of  the  Xe  M-shell  kilovolt  radiation  recorded  with  the  100  fim  pin-hole  is 
illustrated  in  Fig.  (1).  The  z-coordinate  represents  the  direction  of  propagation  (right  to  left) 
and  the  transverse  direction  is  designated  by  the  r-coordinate.  The  spatial  resolution  of  the 
image  is  ~  100  fivn,  the  dimension  of  the  pin-hole.  A  line-out  along  the  axis  [(r,z)  =  (0,z)] 
is  shown  in  Fig.  (2).  A  clear  asymmetry  about  the  focal  position  of  the  248  nm  radiation  is 
present.  Significandy,  the  brighter  zone  is  on  the  M  side  of  the  focus  and  extends  for  a  distance 
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Direction  of  248  nm  Propagation 


Fig.  (1).  Pin-hole  X-ray  camera  image  in  transverse  direction  indicating  the  spatial 
distribution  of  Xe  M-shell  radiation  (~  1  keV).  The  pin-hole  had  a 
diameter  of  -  100  /xm.  The  stagnation  pressure  of  the  pulsed-valve  was 
30  psi  and  the  temperature  was  293  K. 
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Intensity  (orb.  units) 


Rayleigh  Range 


^ _  Direction  of 

248  nm  Propagation 


Fig.  (2):  Line-out  [(r,z)  =  (0,z)]  of  the  data  shown  in  Fig.  (1).  The  position  of  the 

248  nm  focal  point  (z  s  625  nm)  is  indicated  along  with  a  marker 
indicating  the  Rayleigh  range  (~  30  /xm)  corresponding  to  the  focusing 
system  used  for  the  248  nm  radiation.  The  length  5  represents  estimated 
length  of  the  region  in  which  material  from  the  pulse-valve  is  present. 
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present.  Significantly,  the  brighter  zone  is  on  the  side  of  the  focus  and  extends  for  a  distance 
(200  or  300  /im),  a  value  many  times  the  Rayleigh  range  (~  30  /xm)  of  the  optical  system  used 
for  focusing.  Figure  (3)  illustrates  the  image  obtained  with  a  5  /tm  pin-hole  while  Fig.  (4) 
shows  the  corresponding  axitd  line-out.  The  asymmetry  is  clearly  present  in  the  latter  as  in  Fig. 
(2).  Careful  inspection  of  the  data  from  other  exposures  reveals  the  existence  of  a  narrow 
filament  with  a  length  of  —  200  ^m  and  a  radial  extent  limited  by  the  resolution  of  the  imaging 
system  (-20  ;tm).  This  radial  profile  is  shown  in  Fig.  (5),  which  illustrates  a  transverse  line- 
out  of  the  filamentary  region. 

The  images  show  that  the  power  density  of  the  incident  248  nm  radiation  is  confined  in 
the  region  bevond  the  focus.  This  type  of  behavior  would  be  expected  if  the  channeling 
mechanism  discussed  in  Appendices  F  and  G  were  influencing  the  propagation.  All  the 
experimental  data  indicate  that  this  is  the  case. 

IV.  CONCLUSIONS 

The  union  of  a  (1)  multiphoton  cluster/molecule  mechanism  for  the  production  of  X-rays 
and  (2)  a  relativistic/charge-displacement  process  for  the  establishment  of  confined  propagation 
appears  possible.  This  combination  intrinsically  has  exceptionally  favorable  characteristics  for 
the  efficient  amplification  of  X-rays  in  the  kilovolt  spectral  region. 
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Direction  of  248  nm  Radiation 


Fig.  (3):  Pin-hole  X-ray  camera  image  in  transverse  direction  indicating  the  spatial 

distribution  of  Xe  M-shell  radiation  (—  1  keV).  The  pin-hole  had  a 
diameter  of  5  ^m.  The  stagnation  pressure  of  the  pulsed-valve  was  30  psi 
and  the  temperature  was  273  K. 
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Rayleigh  Range 


Line-out  [{r,z)  —  (0,z)]  of  the  data  shown  in  Fig.  (3).  The  approximate 
position  of  the  248  nm  focal  point  (z  =  625  /xm)  is  indicated  along  with 
a  marker  indicating  the  Rayleigh  range  (~  30  /xm)  corresponding  to  the 
focusing  system  used  for  the  248  nm  radiation. 


Intensity  (orb.  units) 
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Abstract.  The  development  of  a  unified  picture  of  short- 
pulse  high-intensity  multiphoion  processes,  embracing 
atoms,  molecules,  and  solids,  appears  possible  through  the 
study  of  clusters.  Of  panicular  significance  are  possible 
intra-cluster  processes  that  can  influence  the  mechanism  of 
ionization  and  lead  to  the  production  of  inner-shell  vacan¬ 
cies.  Inner-shell  excitation  leading  to  prompt  X-ray  emission 
is  specifically  considered  and  the  treatment  leads  to  the  def¬ 
inition  of  a  critical  cluster  size  representing  the  achieve¬ 
ment  of  maximal  X-ray  emission  from  the  ensemble.  These 
results  suggest  the  possibility  of  designing  a  new  class  of 
molecular  materials  optimized  for  the  efficient  production 
and  amplification  of  X-rays. 

PACS:  97.85.  33.00,  36.00,  42.55.Vc 


This  anicle  introduces  and  explores  the  possibility  of  de¬ 
veloping  a  unified  picture  of  short-pulse  (<100fs)  high- 
imensity  (>  10'^  W/cm*)  multiphoton  processes  embracing 
atoms,  molecules,  and  solids,  through  the  study  of  clusters 
[1].  Rare-gas  clusters  {RG)^,  on  account  of  their  simplic¬ 
ity  and  close  connection  with  known  atomic  and  molecu¬ 
lar  properties,  appear  panicularly  suitable  for  this  purpose. 
Of  panicular  significance  are  possible  intra-cluster  processes 
that  can  influence  the  mechanism  of  ionization.  This  can  ( 1 ) 
lead  to  the  production  of  inner-shell  vacancies  and  (2)  the 
generation  of  prompt  X-rays  with  a  potential  for  strong  am¬ 
plification.  The  analysis  developed  below  gives  simple  rela¬ 
tions  concerning  inner-shell  excitation  and  X-ray  emission 
applicable  over  a  laj^e  range  of  cluster  sizes  (1  <  n  <  10^). 
TTie  treatment  also  leads  to  the  definition  of  a  critical  clus¬ 
ter  size  rig  representing  the  achievement  of  maximal  X-ray 
emission  from  the  ensemble,  a  parameter  that  is  found  to  be 
directly  analogous  to  the  concept  of  nuclear  critical  mass. 

Explicit  contact  is  made  with  six  indep)endent  experimen¬ 
tal  observations  [2-4]  involving  three  materials  (Ar,  Kr,  Xe). 
all  either  unexplained  or  anomalous  findings  concerning  the 
generation  of  shon  wavelength  (A  <  200  A)  radiation.  Con¬ 


sistent  agreement  with  the  analysis  is  found  in  each  case. 
Significantly,  this  synthesis  is  achieved  with  the  exclusive 
use  of  previously  established  atomic  and  molecular  prop¬ 
erties.  A  connectipn  with  studies  of  X-ray  generation  from 
solid  surfaces  (BaF2)  is  also  made  in  conjunction  with  this 
analysis. 


1  Coupling  to  Clusters 
1. 1  Simple  Initial  Formulation 

The  rare-gas  cluster  is  considered  as  a  spherical  group  of 
(n)  atoms  having  mutual  interatomic  spacings  and  an 
overall  radius  R  for  the  ensemble  given  by  ^ 

(n  >  3),  and  R  =  rQ(n/2),  (n  =  I,  2),  in  analogy  with 
simple  descriptions  [5,6]  of  nuclei.  Consider  irradiation  of 
the  cluster  at  wavelength  A  with  the  stipulation  (a)  that 
the  characteristic  dimension  /Z  of  a  fully  ionized  cluster 
be  less  than  the  skin  depth  6  and  (b)  the  assumption  that 
the  atoms,  which  are  weakly  bound  by  van  der  Waals 
forces,  respond  to  the  external  optical  field  as  free  atoms 
whose  ionization  can  be  estimated  with  a  tunneling  picture. 
Condition  (a)  enables  all  atoms  in  the  cluster  simultaneously 
to  experience  the  same  field  strength,  since  R  <  8  ^ 
Consistent  with  these  assumptions,  we  will  set  the  maximum 
dimension  R^  20-30  A,  or  equivalently,  a  maximum 
cluster  size  of  10^.  Assumption  (b)  enables  a  simple 

experimentally  evaluated  procedure  to  be  used  to  estimate 
the  level  of  atomic  ionization  produced.  This  gives  the 
threshold  intensity  for  ionization  [7, 8]  as 

4i  =c£7/1287re*Z^  (1) 

where  is  the  ionization  potential,  Z  the  resulting  ionic 
charge,  c  the  speed  of  light,  and  e  the  electronic  charge. 

1 2  Intra-Cluster  Processes 

The  model  of  a  rare-gas  cluster  as  a  sum  of  independent 
atoms,  held  weakly  together  by  van  der  Waals  forces  and 
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responding  as  individual  systems  to  the  external  optical 
field,  is  expected  to  become  invalid  as  the  cluster  becomes 
sufficiently  large  and/or  the  external  electromagnetic  field 
becomes  sufficiently  strong.  This  limit  could  arise  through 
the  action  of  inelastic  scattering  processes  produced  by 
ionized  electrons  originating  from  atoms  in  the  cluster  that 
are  accelerated  by  the  external  optical  field  and  collide  with 
other  atoms  or  ions  in  the  system.  Roughly,  such  intra-cluster 
processes  can  be  considered  as 

(2) 

which  produces  excitation,  or  the  reaction 

^  e-  +  e-  +  ,  (3) 

which  causes  additional  ionization. 

An  approximate  estimate  of  the  number  of  such 
excitations  or  ionizations  can  be  written  as  the  product  of  the 
number  of  electrons  produced  in  the  cluster  (nZ\  the  atomic 
density  in  the  cluster  the  cross  section  for  the  inelastic 
channel  under  consideration  (e.g.  and  the  characteristic 
scale  length  of  the  cluster  (/Z).  Therefore,  for  reaction  (3), 

^  (nZ)«Te,rJ‘^ron'''^  =  ^  ,  (n  >  3) ,  (4) 

with  Z  determined  through  (1)  for  the  corresponding  inten¬ 
sity.  Such  a  mechanism  would  cause  the  pattern  of  ionization 
to  depend  upon  the  cluster  size  (n),  thereby  yielding  an  ob¬ 
served  spectrum  of  ionic  charge  states  deviating  from  that 
characteristic  of  the  free  atom  (n  —  1),  an  easily  measured 
feature  in  direct  conflict  with  the  picture  of  atomic  indepen¬ 
dence  assumed  above  in  the  simplest  initial  formulation. 

An  imp>onani  immediate  consequence  of  (4)  is  the  exis¬ 
tence  of  a  critical  size 


for  inner-shell  excitation 


STran^/^  7^7*5  ’ 


(n  >  3) , 


(6) 


with  a  =  e^fhc.  Since  this  simple  picture  assumes  that  the 
acceleration  occurs  in  a  time  comparable  to  or  less  than  a 
period  of  the  wave,  validity  of  this  approximation  demands 
that  the  cluster  neither  be  too  large  nor  the  external  field  too 
small:  namely,  it  requires  that  the  amplitude  of  the  driven 
excursion  of  the  electron  caused  by  the  field  be  comparable 
to  or  larger  than  the  cluster  dimension  iZ.  If  we  set  the 
limiting  situations  as  =  iZ,  assume  that  the  electron  is 
free,  and  put  x^  =  eElmuP-,  with  m  denoting  the  electron 
mass  and  uj  the  angular  frequency  of  the  external  optical 
field,  we  can  define  another  limiting  intensity  as 


W  = 


27r^ 

a 


c 


(n>3), 


(7) 


with  A  denoting  the  wavelength  of  the  optical  field  and 
being  the  electron  Compton  wavelength.  In  this  picture,  an 
inner-shell  electron  with  binding  energy  can  be  ionized  if 
the  intensity  exceeds  a  lower  bound  given  by  /q  or 
whichever  is  greatest. 


2  Prompt  X-Ray  Emission 

Excitation  of  the  type  discussed  above  can  lead  to  the 
emission  of  prompt  inner  0  -1)  shell  radiation  in  the  X- 
ray  range.  Generally,  in  addition  to  the  production  of  a  hole 
in  the  0  —  U  shell,  the  possibility  of  prompt  j  j  —  \ 
transitions  requires  retention  of  at  least  one  electron  in  the 
j -shell  during  the  course  of  iiradiauon.  This  determines  an 
upper  bound  on  the  intensity  of  irradiation  through 

(1)  which  is,  in  the  simple  model  considered,  independent 
of  the  cluster  size. 


defined  by  =  n.  Therefore,  if  corresponds  to  an  inner- 
shell  ionization,  with  n  >  an  inner-shell  vacancy  would 
be  produced  in  every  atom.  Furthermore,  since  the  critical 
size  of  the  cluster  depends  inversely  on  the  square  of  the 
atomic  density  inversely  on  the  cube  of  the 

cross  section  for  excitation  (cr~^),  is  a  direct  analogue  of 
the  nuclear  concept  of  a  critical  mass,  with  replacing  the 
neutron  fission  cross  section  cr^  and  Z  playing  the  role  of 
the  average  neutron  yield  u. 

Certain  conditions  must  be  met  in  order  for  intra-cluster 
processes  of  the  type  expressed  by  reactions  (2)  and  (3) 
to  occur  for  an  inner-shell.  Specifically,  such  an  excitation 
can  only  take  place  if  sufficient  kinetic  energy  is  possessed 
by  the  colliding  electron.  To  estimate  this  energy,  and  the 
conditions  necessary  for  its  transfer,  we  will  assume  that  the 
electrons  experience  an  acceleration  over  the  characteristic 
dimension  of  (iZ)  of  the  cluster  by  the  external  optical  field 
(E)  in  a  time  shon  compared  to  the  period  of  the  oscillating 
wave.  Funhermore.  once  the  electrons  leave  the  vicinity  of 
the  cluster  we  will  assume  that  they  represent  ionization 
and  are  lost  from  the  system.  This  simple  procedure  yields  a 
characteristic  energy  given  by  e,  S  eER,  a  result  that  can 
be  expressed  as  a  characteristic  minimum  intensity  needed 


2.J  Krypton 


Considering  Kr  A/-shell  emission  as  an  example,  from 
the  threshold  ionization  curve  [7-9]  derived  from  (1)  and 
exhibited  in  Fig.  1.  we  find  that  I^{N)  S  2  x  10'*  W/cm^. 
By  combining  the  conditions  stemming  from  (6)  and  (7) 
^8-  2  illustrates  the  conditions  on  intensity 
and  cluster  size  for  which  the  model  allows  prompt  M-shell 
emissions.  The  region  corresponding  to  prompt  L-shell  [10] 
emissions  =  1.67  keV)  is  also  represented  in  Fig.  2. 

Several  imponant  results  are  manifest  in  the  form  of 
Fig.  2.  First,  the  crossing  points  np(P)  and  found 

by  equating  (1)  with  (6)  and  (7),  respectively,  yield  the 
condition  for  the  existence  of  an  allowed  zone  for  prompt 
emission  as 


(ris)'®’© 


>  1 


(8) 


an  expression  which  we  observe  is  independent  of  Tq. 
Furthermore,  (8)  strongly  favors  the  siniation  in  which  the 
magnitudes  of  the  characteristic  energies,  and  hence,  the 
corresponding  levels  of  ionization  Z.  are  high.  This  is 
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Classical  MPI  of  KR 


Charge  Stote 


Fig.  1.  Threshold  intensity  for  ionization  shown  as  a  function  of 
resulting  charge  state  Z  calculated  for  Kr  with  the  classical  method 
described  in  [7,8].  The  respective  atomic  shells  (/C,  L,  A/,  .V)  are 
designated.  The  maximum  intensities  and  associated 

with  leaving  some  A/-sheII  and  iV-shcll  populations,  respectively,  are 
indicated.  Also  shown  is  the  intensity  (--5  x  10'^  W/cm“)  at  which  Kr^'*' 
transitions  at  -^-ICWA  were  observed  in  fluorescence  in  experiments 
reponed  in  [4,  11] 


Cluster  Size  (n) 


O  Kr  M-ahell  (-100  •V)  ©UVM.  Exp«Tm«n 
H  Kr  L-aheil  K*V)  Q  Forttddan  Zone 

Fig.  2.  Representation  of  limiting  conditions  applicable  for  production 
of  prompt  Kr  L-shell  and  A/-shcll  emission  from  Kr  clusters  as  a 
function  of  intensity  of  irradiation  and  cluster  size.  The  expcnmental 
datum  (LAND  corresponding  to  Lee  et  al.  [2]  and  the  crossing  points 
P  and  Q  are  shown.  The  additional  condition  7^(248  nm),  stemming 
from  ( 7).  is  also  represented.  See  text  for  other  details 

readily  seen  by  letting  be  approximated  by  hydrogenic 
binding  (E^:xZ^),  keeping  '-^constant,  and  noting 

that  follows.  As  a  natural  consequence  of  the 

sysiematics  of  atomic  binding  energies,  we  observe  that  large 


Wavelength  (A) 


Fig.  3a.  b.  Spectra  of  Kr  in  the  region  near  100  A.  a  Multiphoton 
produced  spectrum  of  Kr  produced  in  a  gas  jet.  as  described  in 
[11],  in  the  ^85-100  A  region  with  a  peak  intensity  of  irradiation 
of  ^5  X  I0**W/cm^.  The  4p  — ►  3d  features  from  Kr’"*'  appear 
prominently,  b  Spectrum  of  Kr  observed  in  a  theta  pinch  plasma  by 
Jones  and  Kailne  [12].  The  Kr^*^  4p  —  3d  transitions  are  indicated.  The 
correspondence  between  the  muitiphoton  and  the  theta  pinch  produced 
spectra  involving  Kr’'*'  in  the  100  A  region  is  marked 


allowed  zones  are  associated  with  level  structures  in  the 
kilovolt  range,  a  feature  apparent  from  comparison  of  the 
regions  of  M-shell  and  L-shell  emission  shown  in  Fig.  2. 

The  information  in  Fig.  2  can  be  directly  related  to  rwo  in¬ 
dependent  experimental  studies,  one  [2]  involving  the  search 
for  prompt  multiphoton-induced  X-ray  emission  from  free 
Kx  atoms  (n  =  1)  and  another  concerning  the  measurement 
of  fluorescence  produced  in  a  pulsed  gas  target  undergoing 
multiphoton  excitation  [4, 11].  The  datum  corresponding  to 
the  former  [2]  (/  S  3  x  10*^  W/cm*,  n  =  1),  designated 
as  LANL,  appears  on  Fig.  2.  Although  the  experiment  ac¬ 
tually  represents  all  intensities  I  <  10*^  W/enr  due  to  the 
distribution  of  intensity  occurring  in  the  focal  region,  we  see 
from  Fig.  2  that  there  is  no  overlap  of  the  experimental  con¬ 
ditions  with  either  of  the  zones  for  prompt  L-  or  M-shell 
emissions.  The  experimental  region  lies  well  below  the  L- 
shell  region  and  is  excluded  from  the  M-shell  zone  by  the 
crossing  P.  Hence,  prompt  X-ray  emission  is  forbidden,  a 
prediction  fully  in  agreement  with  the  experimental  obser¬ 
vation  of  none  [2]. 

The  second  experiment  involves  the  anomalous  detec¬ 
tion  [4. 11]  of  intense  Kr  M-shell  fluorescence  from  Kr^'*' 
4p  —  3d  transitions  at  ^  100  A,  as  illustrated  in  Fig.  3a, 
arising  from  material  produced  in  a  pulsed  gas  jet  and  irra- 


15 


340 


A.  McPhenon  et  ai. 


•diated  at  a  peak  intensity  of  -^5  x  10*^  W/cm"  with  subpi¬ 
cosecond  248  nm  radiation.  The  corresponding  spectral  re¬ 
gion  observed  from  a  theta  pinch  by  Jones  and  Kallne  [12] 
is  shown  in  Fig.  3b. 

The  observation  of  any  radiation  by  the  multiphoton 
excitation,  either  prompt  or  delayed,  immediately  presents 
a  paradox.  As  shown  in  Fig.  I,  with  a  peak  intensity  of 
~5  X  10^^  W/cm-,  no  production  of  is  possible,  a 

result  that  is  experimentally  based  on  earlier  measurements 
of  ion  production  [8]  involving  free  krypton  atoms  under 
identical  conditions  of  irradiation.  Yet  the  Kr^”^  array  at 
103  A  was  among  the  brightest  features  observed  throughout 
the  entire  region  snidied  [4, 11].  However,  Fig.  2  shows  that 
the  formation  of  small  clusters  (n  >  2-3)  would  permit  both 
anomalous  ionization  (Kr^*^  production)  and  the  emission 
of  Kr^-^  4  —  3  radiation  of  the  kind  illustrated  in  Fig.  3, 
thereby  reconciling  the  paradox.  Funhermore,  a  review  of 
the  hydrodynamic  propenies  and  operating  conditions  [4, 1 1] 
of  the  pulsed  valve,  particularly  in  comparison  to  other 
work  [13-15]  in  which  Kr2  and  Kr3  species  were  formed 
and  directly  measured  [14],  showed  that  the  hydrodynamic 
conditions  were  essentially  the  same  in  both  snidies.  Indeed, 
the  observed  pressure  dependence  of  the  Kr^"'“  signal  was 
close  to  that  seen  in  the  earlier  studies  [14]  involving  Kr-, 
and  Kr3  formation.  Both  of  these  considerations  lead  firmly 
to  the  conclusion  that  significant  production  of  small  Kr^ 
clusters  was  inevitable  in  the  multiphoton  experiments  in 
which  the  fCr^"^  transitions  were  observed. 


At 


Fig.  4.  Representation  of  limiting  conditions  for  excitation  of  prompt 
Ar  L-shell  radiation  from  Ar  clusters  as  a  function  of  intensity  of 
irradiation  and  cluster  size.  The  .U-shell  limit  is  taken  from  Rg.  7. 
The  allowed  zone  is  shaded.  Points  P  and  Q  designate  the  crossings 
described  in  the  texL  The  conditions  /n(248  nm)  and  7^(765  nm), 
stemming  from  (7),  arc  also  indicated.  The  experimental  datum  (LANL) 
of  Lee  et  ai.  [2]  is  shown  along  with  the  maximum  intensity  {^5  x 
10*^  W/cm“)  used  in  fluorescence  snidies  of  Ar  discussed  further  in  [3] 


22  Argon  and  Xenon 

t  Additional  observations  [2^]  with  Ar  and  Xe  can  be  ex¬ 
plained  mutatis  mutandis  and  the  corresponding  diagrams 
appear  for  these  materials  in  Figs.  4  and  5,  respectively.  The 
results  [2]  of  Lee  et  al.  (LANL  data),  which  show  no  prompt 
emission,  are  consistent  in  all  cases  illustrated  for  both  ma¬ 
terials.  Note  the  exclusion  of  the  atomic  region  (n  =  1)  for 
Ar  L-shell  and  Xe  N -shell  emissions  by  the  crossings. 

Again,  however,  the  exp)eriments  examining  the  multi¬ 
photon-induced  fluorescence  with  the  pulsed  jet  indicate  very 
different  behavior.  Fluorescence  [3]  from  Ar,  observed  with 
a  maximum  intensity  of  irradiation  of  ^>^5  x  10*^  W/cm-, 
detected  the  clear  presence  of  radiation  at  ^  165  A  on  the 
well-known  doublet  2s2p^  —  2r2p^  transition  [16]  of  Ar^“^. 
This  finding,  as  shown  in  Rg.  6,  is  another  paradoxical  re¬ 
sult,  since  Fig.  7  indicates  that  the  threshold  intensity  of 
Ar^*^  production  is  10**  W/cm-,  approximately  20-foId 
higher  than  that  used  in  the  experiment.  Furthermore,  spec¬ 
troscopic  analysis  of  the  plasma  [3]  led  to  the  conclusion 
that  electron  collisions  could  not  have  produced  the  Ar^'^ 
species.  In  this  connection,  we  note  that  the  cross  section 
for  electron  collisional  ionization  of  Ar****  is  quite  small 
[17].  However,  as  shown  in  Rg.  4,  small  cluster  forma¬ 
tion  (n  =  2. 3)  would  allow  the  excitation  of  an  L-sheil 
electron.  Since  the  maximum  exp)erimenial  intensity  is  also 
above  the  ;V/-shell  limit,  this  interaction  could  readily  gen¬ 
erate  the  Ar^"**  species  seen.  Therefore,  small  cluster  forma¬ 
tion  removes  the  anomaly  of  the  Ar^***  production  and  the 
observed  radiation  on  the  165  A  doublet  could  then  be  gen¬ 
erated  in  the  manner  previously  described  [3]. 


Cluster  Size  (n) 

Qxa  (S-6  k«V) 

S  X«  M  (CXa  -  IX)  k«V)  LAM.  Datun 

■  x«  K60  -  00  mV)  [jFffttdd^  Zonn 

^x*  aolld  stats  studlsal 

Rg.  5.  Representation  of  limiting  conditions  applicable  for  excitation  of 
prompt  Xe  AT-shcll.  Af -shell  and  I-shcll  radiation  from  Xe  clusters  as 
a  function  of  intensity  of  irradiation  and  cluster  size.  The  experimental 
datum  (LANL)  corresponding  to  Lee  et  al.  [2],  which  falls  in  the 
forbidden  zone,  is  shown.  The  boundaries  are  established  as  previously 
described  for  Kr  and  Ar  in  Rgs.  2  and  4,  respectively.  The  points  A^, 
By,  and  correspond  approximately  to  the  conditions  pcnaining  to 
the  data  in  Rgs.  1 1, 12,  and  13,  respectively.  Sec  text  for  details 
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Fig.  6.  Transitions  observed  in  Ar  with  a  puised-gas  jet  at  an  intensity 
of  --5  X  W/cm-  at  248  nm  as  reponed  in  [3].  The  Zs^2p^- 
2s^2p6  doublet  is  clearly  visible  along  with  two  unidentified  (UI)  lines. 
The  latter  may  be  associated  with  satellite  transitions  in  Ar®"*" 


Xe 

BOS  X©^  4d*5s5p-*4d^5s 
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Fig-  8.  Emission  spectrum  of  Xe  in  the  144-174  A  region  observed 
in  a  puised-gas  jet  at  a  peak  intensity  of  -5  x  10‘*  W/cm*  with 
subpicosecond  248  nm  in^adon.  Features  associated  with  Xe^"^. 
Xe*"^,  and  Xe^'*"  are  identified 


Xe*^  4d*5p-**4d^ 
i/z/A  Xe*^  4d*5p^4d*  3 


Classical  MPl  of  AR 


Fig.  7.  Threshold  intensity  for  ionization  shown  as  a  function  of 
resulting  charge  state  Z  calculated  for  Ar  with  the  classical  method  de- 
senbed  in  [7].  The  respective  atomic  shells  (L  and  \f)  are  designated. 
The  maximum  intensity  /^(Af)  associated  with  leaving  some  A/-shelI 
population  is  indicated.  Also  shown  is  the  threshold  for  Ar^^  produc¬ 
tion  10**  W/cm-)  and  the  maximum  intensity  (~5  x  10**  W/cm-) 
utilized  in  fluorescence  studies  [3]  of  Ar  using  a  puised-gas  jet  in  which 
the  165  A  2s2p*  —  2s'^2p^  transition  in  Ar^*^  was  observed 


Classical  MP!  of  X£ 


Rg-  9.  Threshold  intensity  shown  as  a  function  of  resulting  charge 
state  Z  calculated  for  Xe  with  the  classical  method  described  in  [7]. 
The  respective  atomic  shells  (/C,  L,  Af,  iV.  O)  are  designated.  The 
maximum  intensities  /^(Af),  1^{N)  and  I^{0)  associated  with 
leaving  some  .\f-sheH.  .V-shell  and  O-shcll  populations,  respectively, 
are  indicated.  The  experimental  intensity  (/^  St  5  x  10**  W/cm*) 
used  in  studies  of  fluorescence  from  Xe  using  a  pulscd-gas  jet  is 
indicated.  This  intensity  is  just  sufficient  to  lead  to  observable  Xe^'*’  in 
ion  measurements 


A  similar  analysis  can  be  applied  to  Xe  I-shell.  .Vf-shell. 
and  .V-shell  emissions  which  yields  conclusions  analogous 
to  those  found  for  Ar  and  Kr.  The  corresponding  allowed 
zones  for  prompt  emission  from  those  shells  are  indicated 
in  Fig.  5. 

The  multiphoton-induced  fluorescence  observed  [11. 18] 
from  Xe,  panicularly  from  levels  involving  4d  excitation. 


suggests  additional  important  insights.  Transitions  [19-21] 
originating  from  states  involving  a  Ad  hole  were  observed 
in  three  species,  Xe^+,  Xe*-*",  and  Xe’+.  The  identified 
configurations  of  the  lines  shown  in  Fig.  8  were  ^5s5v 
4d‘05s(Xe^*),  4d’5p  ^  4d“>(Xe®+).  and  4d«5p  - 

4d’(Xe’'^).  The  maximum  intensity  used  at  248  nm  was 
~5  X  10'*  W/cm-  which  is  just  sufficient,  as  seen  from  Fig.  9, 
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to  lead  to  observable  Xe^‘'“  production  in  ion  charge  state 
measurements  [18]. 

Experimentally,  these  transitions  presented  two  obviously 
anomious  features,  (1)  their  exceptional  strength  and  (2)  the 
nature  of  the  radiating  configurations  delected.  The  Xe^*^ 
line  at  165.3  A  was  by  far  the  most  intense  feature  seen  in 
any  of  our  studies  of  fluorescence,  indicating  that  copious 
4d  excitation  was  certainly  being  produced.  Funhermore, 
consider  the  configuration  [Ad^5p]  of  the  excited  state  in 
Xe^"^  producing  the  array  in  the  ">^150A  region.  This  level 
has  two  4d  holes,  but  retains  a  relatively  weakly  bound  5p 
electron.  Production  of  this  level  by  recombination  from 
Xe^°'^(4d®)  is  not  possible,  since  the  intensity  was  too  low 
for  its  production  [18].  Collisional  excitation  by  sufficiently 
energetic  electrons  in  the  plasma  formed  is  possible,  but 
calls  for  an  unusually  high  efficiency  given  the  very  high 
observed  intensity  of  the  xenon  transitions.  Furthermore,  it 
was  not  possible  to  reconcile  in  any  reasonable  way  the 
observed  density  dependence  of  the  emission  with  collisions 
arising  from  plasma  electrons.  In  fact,  as  in  the  case  of  Kr^"** 
radiation  discussed  above,  the  observed  pressure  dependence 
was  in  far  closer  accord  with  that  associated  with  the 
formation  of  small  rare  gas  clusters  [14,22]  arising  from 
the  complex  hydrodynamics  associated  with  the  flow. 

The  apparently  anomalous  4d  excitation  implied  by  the 
Xe  fluorescence  would  be  expected,  however,  if  the  cluster 
process  discussed  above  in  relation  to  Ar  and  Kr  is  operative. 
In  Fig.  5  the  triangular  zone  representing  the  allowed  region 
for  prompt  .V-shell  emission,  as  shown  enlarged  in  Fig.  10, 
indicates  that  dimer  (Xe^)  formation  would  be  sufficient 
to  permit  prompt  emission  to  occur  in  the  regions  of  the 
focus  for  which  the  intensity  corresponds  to  the  correct  range 
[/  <  Eniission  from  Xe^*^  and  Xe®'‘“  species  having 

a  4d  excitation  could  be  expected  to  originate  from  the 


Fig.  10.  Expansion  of  Fig.  5  showing  the  parameiers  concerning  4<f- 
subshell  excitation  in  Xc.  The  expenmentai  intensity  (/„p  5  x 

lO**W/cm“)  used  in  the  studies  of  fluorescence  is  indicated.  The 
limiting  curve  /,j(248nm)  is  shown.  The  O-sheil  limit  is  taken  from 
Fig.  9  and  the  allowed  zone  for  prompt  .V-shcll  emission  is  shown. 
The  critical  cluster  size  (n^  S  5)  corresponding  to  4d-shcU  excitation 
for  the  prevailing  experimental  conditions  is  indicated 


small  clusters  (2  <  n  <  20)  corresponding  to  the  indicated 
allowed  zone.  Anomalous  ionization  would  also  be  expected, 
which  could  possibly  produce  Xe^^(4<i*),  and  consequently, 
through  recombination  in  the  plasma,  the  Xe^  4d*5p  — ►  ^ 
lines  at  ~  150  A.  However,  the  information  shown  in  Fig.  10 
suggests  the  possibility  of  another  mechanism  of  production 
for  the  Xe^'‘’(4d*5p)  state. 

For  electrons  with  a  kinetic  energy  of  ~100eV,  the 
cross  section  for  4d  ionization  in  Xe  is  unusually  high,  a 
feature  arising  from  the  shape  resonance  involving  the  5p, 
5s,  and  4d  shells  well  known  from  studies  of  pholoionization 
and  electron  collisional  ionization  in  that  atom  [23—30]. 
The  measurements  establish  a  maximum  cross  section  of 

=  1.6  X  10“*^  cm*  for  4d  electrons  [24].  A  comparably 
high  value  [23,31]  also  holds  for  ions  such  as  Xe**^.  Since 
the  value  Z  =  7  corresponds  to  the  production  of  Xe®'*',  we 
can  now  evaluate  the  critical  cluster  size  through  (2).  With 
[32]  Tq  =  4.4  A,  we  obtain  =  5,  a  remaricably  small 
value  and  one  that  matches  the  allowed  region  shown  in 
Fig.  10.  Therefore,  since  the  critical  cluster  size  is  so  small 
in  this  case,  the  production  of  single  4d  vacancies  can  be 
nearly  saturated  even  in  relatively  small  clusters,  and  the 
production  of  double  vacancies  could  become  significant. 
Hence,  direct  double  vacancy  production  induced  by  the 
intra-cluster  mechanism  becomes  a  possibility  for  the  origin 
of  the  Xe^‘’"(4d*5p)  species. 


23  Solid-State  Studies 

The  results  for  Xe  presented  in  Fig.  5,  with  (a)  the  ne¬ 
glect  of  the  small  difference  in  atomic  number  separating 
Xe(Z  =  54)  and  Ba(Z  =  56),  (b)  the  further  neglect  of  the 
relatively  small  difference  in  atomic  spacings,  and  (c)  the  as¬ 
sumption  that  71  10^  corresponds  to  a  cluster  having  prop¬ 

erties  approximately  those  of  bulk  solid  matter,  enable  us 
to  make  contact  with  data  recently  obtained  in  experiments 
studying  X-ray  generation  from  BaFj  surfaces.  The  spectra 
of  interest  are  shown  in  Figs.  11, 12  and  13.  The  conditions 
penaining  to  the  data  shown  in  Fig.  1 1  are  represented  ap¬ 
proximately  by  the  point  in  Fig.  5.  Since  this  point  lies 
at  the  border  of  the  M-shell  allowed  zone,  it  carries  the 
strong  implication  that  anomalous  features  stemming  from 
the  type  of  direct  coupling  discussed  above  may  actually  be 
observable  in  the  spectrum  obtained  with  solid  BaF2. 

Transitions  involving  the  Ba  M-shell  arising  from  this 
form  of  interaction  would  be  expected  (1)  to  appear  as 
relatively  narrow  features,  (2)  exhibit  observable  satellite 
structure  due  to  the  distribution  of  electrons  presumably 
produced  in  the  iV-shell  of  the  system  during  the  interaction, 
and  (3)  display  a  considerably  difrerent  dependence  on 
the  conditions  of  irradiation  than  the  broad  quasiconiinuum 
discussed  in  [33]  and  [34].  Indeed,  an  anomalous  feature 
(A^  =  13.40 ±0.05  A)  exhibiting  all  three  of  these  expected 
characteristics  is  seen  in  Fig.  11.  The  transition  appears  to 
be  associated  with  weak  satellite  structure  and  it  was  found 
experimentally  to  be  quite  sensitive  to  the  conditions  of 
irradiation,  a  feature  illustrated  by  the  comparison  of  Fig.  1 1 
with  Figs.  12  and  13.  Higher  intensity,  corresponding  to 
point  Bj.  in  Fig.  5,  gives  relatively  less  radiation  at  A^.,  and 
irradiation  with  a  long  (~2  ns)  low  intensity  (<  10*^  W/cm*) 
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Fig.  11.  Spectrum  of  BaF^  irradiated  with  subpicosecond  f'>-270fs) 
ultraviolet  (248  nm)  radiation  at  a  peak  intensity  of  -^8  x  10'*  W/cm- 
observed  with  a  von  Hamos  curved  crystal  spectrometer  on  a  single 
exposure.  Occumng  along  with  the  Ba  .Vf -shell  continuum,  the 
Lq(14.98A)  line  and  an  unidentified  feature  at  13.40  ±  0.05  A  are 
seen.  These  data  were  recorded  in  Chicago  by  A.  McPherson  and 
A.P.  Sheverko.  We  acknowledge  the  assistance  of  P.G.  Burkhalter 
and  D.A.  Newman  of  the  Naval  Research  Laboratory  in  obtaining  the 
densitometer  trace  of  these  data.  Assuming  n  S'  10^,  the  conditions  of 
these  data  are  represented  on  Fig.  5  approximately  by  the  point 


11.00  Wavelength  [A]  16.00 


Fig-  12.  Spectrum  of  BaF-.  irradiated  as  in  Rg.  1 1  except  with  an 
incident  intensity  of  10’^  W/cm*.  The  F*"^  transition  appears 

stronger,  but  the  unknown  feature  has  become  relatively  weakened. 
These  data  correspond  approximately  to  point  in  Rg.  5 


0  to  11  12  IJ  U  IS 
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Fig-  13.  Spectrum  of  BaF^  irradiated  by  1.06  pm  radiation  with 
a  2  ns  pulse  length  at  an  intensity  <  10'^  W/cm*.  The  .V/-shell 
quasicontinuum  is  present,  but  the  unknown  A^  feature  is  absent. 
These  data,  which  were  recorded  at  the  Lebedev  Physical  Institute, 
Moscow,  by  A.  Shevel’ko  and  A.  MePhenon  in  April  1993,  correspond 
approximately  to  the  region  below  point  in  Rg.  5 
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pulse,  corresponding  to  a  position  far  removed  from  the 
Af -shell  allowed  zone  in  Fig.  5  (region  below  point  C^), 
leads  to  a  vanishing  of  the  line.  Finally,  theoretical 
calculations  performed  by  the  Moscow  group  (A.  Shevel’ko 
and  co workers)  indicate  that  A^  corresponds  best  to  the 
3d*®45“  —  3(^4$Hf  transition,  a  transition  whose  nature 
could  be  expected  with  the  proposed  method  of  coupling 
by  3d  ^  4/  excitation  through  reaction  (2). 

Two  additional  comments  relevant  to  the  data  shown 
in  Fig.  11  can  be  made.  First,  the  unknown  feature  A^  is 
definitely  not  weak,  since  the  9-13  A  continuum  [33,34] 
against  which  it  must  compete,  is  very  strong.  Second,  if  the 
transition  is  associated  with  a  direct  process  of  excitation 
of  the  nature  discussed  above,  then  it  can  only  occur  in 
the  early  stages  of  the  development  of  the  hot  plasma, 
since  the  high  temperature  plasma  would  quickly  lead  to 
full  ionization  of  the  iV-shell.  On  the  basis  of  LASNEX 
simulation  [34]  of  the  experiments  [33]  on  Ba  in  which  the 
electron  temperature  is  calculated,  the  relevant  time  appears 
to  be  approximately  the  initial  ~50-l(X)fs  of  irradiation.  It 
is  further  implied  that  the  excitation  would  also  be  confined 
to  a  region  quite  close  to  the  surface,  perhaps  in  the  first 
skin  depth.  Hence,  such  transitions  might  be  able  to  serve 
as  highly  spatially  localized  probes  of  the  interaction  for 
times  sufficiently  short  that  equilibrium  conditions  cannot 
be  established.  Of  course,  these  observations  only  give  a 
hint  at  the  possibilities,  and  much  further  study  is  needed  to 
verify  or  deny  these  hypotheses  and  possibilities. 

3  Implications  for  X-Ray  Amplification 

Overall,  six  previously  unexplained  or  anomalous  experi¬ 
mental  observations  involving  three  materials  (Ar,  Kr,  and 
Xe)  can  be  accounted  for  by  the  formation  of  small  clus¬ 
ters.  In  addition,  a  possible  connection  with  X-ray  generation 
from  solids  has  been  made.  Furthermore,  it  can  be  similarly 
shown  [35]  that  both  the  observed  [36]  production  of  2a ^ 
excitation  in  Nj,  involving  subsequent  emission  at  55.8  nm, 
and  the  anomalous  production  [37]  of  from  UF^  are 
consistent  with  the  picture  developed. 

In  aggregate,  these  results  have  two  imponant  implica¬ 
tions.  One  involves  the  possibility  of  developing  a  unified 
picture  of  high-intensity  multiphoion  processes  applicable 
for  sufficiently  shon  pulses  to  atoms,  molecules,  and  solids. 
The  second  concerns  the  concept  of  a  molecular  design  for 
systems  optimized  for  efficient  X-ray  production  and  am¬ 
plification.  The  latter  results  from  the  potential  confluence 
of  several  favorable  factors.  They  are  (I)  the  inner-shell 
specificity  of  the  excitation  expressed  in  Figs.  2, 4, 5,  and 
10,  (2)  the  high  efficiency  of  coupling  potentially  indicated 
(n  ^  n^),  and  (3)  the  ability  to  combine  simultaneously  the 
propitious  qualities  of  both  high  and  low  density  in  the  dy¬ 
namic  of  the  medium.  The  selectivity  (1)  and  the  efficiency 

(2)  are  consequences  of  the  fact  that  the  electronic  motions 
involved  in  the  coupling  occupy  only  a  relatively  small  frac¬ 
tion  of  the  available  phase  space.  With  respect  to  the  laner 

(3) ,  the  use  of  clusters  achieves  this  paradoxical  combi¬ 
nation  of  opposites  automatically.  Namely,  clusters  enable 
atomic  excitation  rates  to  achieve  the  extraordinarily  high 
values  characteristic  of  condensed  maner  while  simultane- 
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^usly  eiiminaiing,  through  local  rapid  expansion  to  low  par¬ 
ticle  density,  undesirable  collisional  interactions  that  would 
establish  equilibrium  conditions  in  the  excited  material. 

The  interaction  also  leads  naturally  to  the  production  of 
inverted  population  densities.  Consider  the  case  of  Ar.  as 
illustratedjn  Fig.  7,  irradiated  at  an  intensity  that  is  just 
below  In  this  circumstance,  seven  electrons  are 

removed  from  the  M -shell,  producing  with  the  config¬ 
uration  ls“2s-2p^3s,  but  no  Ar^"^  having  the  configuration 
l5-2s-2p^  can  be  formed.  However,  as  shown  in  Fig.  4,  in 
(Ar)^  clusters  with  n  >  3,  an  I-shell  electron  can  be  addi¬ 
tionally  liberated  at  the  same  intensity  with  the  outcome 

that  Ar^"^  can  be  produced,  but  only  with  the  excited  con¬ 
figuration  ls“2s-2p^3s.  Funhermore,  if  we  have  n  > 
every  atom  in  the  cluster  is  produced  in  the  excited  (Ar^*^)* 
state  and  no  ground  configuration  Ar^"*"  species  is  formed. 
Basically,  the  generation  of  the  ground  Ai^"^  state  is  totally 
blocked,  but  the  channel  leading  to  the  excited  level  is  fully 
open.  Therefore,  in  addition  to  the  selective  and  efficient 
production  of  radiation  on  the  3s  2p  transition  of  Ar^*^, 
these  dynamics  also  automatically  generate  inverted  popula¬ 
tion  distributions  with  the  capacity  for  amplification. 

Intra-cluster  collisional  ionization  of  the  M-shell  elec¬ 
trons  by  reaction  (3)  can  potentially  destroy  the  inversion 
by  enabling  the  generation  of  Ar^"^  ls-2s-2p^  ground  state 
species  in  the  example  considered  above.  However,  the 
use  of  heavy  materials,  which  exhibit  shape  resonances  for 
inner-shell  excitation  (e.g.  Xe  4{f)  and  which  have  large  cor¬ 
responding  cross  sections  ~  10“*^  cm*),  appears  to 
sharply  reduce  the  significance  of  such  losses.  Estimates 
of  the  cross  section  for  the  loss  channel  noted  above  (A/- 
shell  collisional  ionization),  computed  for  relevant  condi¬ 
tions  using  the  Bethe  approximation  [38],  give  an  estimate 
of  -  10"’**  cm*,  so  the  shape  resonances  correspond¬ 

ing  to  inner-shell  excitation  can  still  dominate.  Basically,  this 
question  reduces  to  the  issue  of  whether  ionization  cross 
sections,  as  a  function  of  energy,  can  exhibit  significant 
edges  with  respect  to  inner-shell  excitation.  As  shown  in 
Figs.  14  and  15,  the  experimental  evidence  demonstrates  that 
when  the  inner-shell  interaction  involves  a  shape  resonance 
[3 1, 39],  such  edges  can  be  prominent.  However,  the  dynam- 


Fig.  14.  Single  ionization  cross  sections  of  Xc-^  as  a  function  of 
electron  energy  {solid  circles).  Also  shown  are  expen  mental  results 
of  Griffin  et  ai.  (X)  and  Achenbach  ct  al.  (-f-).  The  curve  is  the 
semiempirical  calculation  using  the  Lotz  formula.  Rgure  and  caption 
taken  from  (31)  which  should  be  consulted  for  funher  details.  Figure 
used  with  permission 
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Fig.  15.  Absolute  cross  sections  for  single  ionization  of  Ba-^  ions. 
Open  squares  represent  measurements  employing  a  “cold”  Penning  ion 
source  (low  content  of  metastable  ions  in  the  parent  beam);  the  open 
circles  are  measurements  employing  a  “hot”  ECR  ion  source  (producing 
a  substantial  amount  of  metasiable  ions  in  the  parent  beam).  The 
statistical  uncenainties  are  smaller  than  the  size  of  the  symbols.  The 
solid  line  is  a  distorted-wave  calculation  for  direct  ionization  of  the  5p 
and  5s  subshells  of  ions  in  the  5s-5p^  ground-state  configurauon.  The 
chain  curve  is  a  distoned-wave  calculation  for  direct  ionization  of  the 
5p  and  5s  subshells  of  ions  in  the  5s-5p^5d  metasuble  configuration. 
Figure  and  caption  taken  from  [39]  which  should  be  consulted  for 
funher  details.  Figure  used  with  permission 


ics  of  these  collisional  amplitudes,  which  will  (3ccur  when 
the  strong  radiative  field  is  present,  can  only  be  reliably 
tested  by  experiment. 

Auger  processes  can  also  generate  losses.  In  this  case,  the 
Auger  rate  depends  quite  sensitively  on  the  configuration 
[40],  a  propeny  that  can  be  used  to  suppress  strongly 
the  Auger  decay.  For  the  short  pulse  excitation  ('^-lOOfs) 
being  considered,  Ar  ions  with  the  excited  state  configura¬ 
tions  (l5*2s*2p^3s-),  Ar^-" 

(ls“25*2p^3s3p),  and  Ar**^  (ls*2s‘2p^3s)  appear  to  have 
sufficiently  reduced  Auger  rates.  Overall,  the  figures  on 
both  collisional  and  Auger  losses  suggest  that,  in  properly 
selected  systems  undergoing  irradiation  at  an  appropriate 
intensity,  the  loss  channels  may  be  a  small  fraction  (<  10%) 
of  the  total  excitation. 

This  situation  could  lead  to  very  significant  levels  of 
amplification.  If  we  assume  that  the  amplifiying  X-ray 
transition  is  Doppler  broadened,  on  account  of  the  ionic 
velocity  arising  from  the  coulomb  explosion,  ignore  the 
collisional  and  Auger  loss  channels,  and  take  the  cluster  size 
n  =  n^.,  it  can  be  shown  [35]  that  the  gain  constant  per 
unit  electron  density  at  wavelength  is  given  by 


n,  SnZ\/2  c  r5  V  /  Vooy 


(9) 


for  homonuclear  clusters  composed  of  atoms  with  mass  M 
having  an  X-ray  spontaneous  emission  rate  .4;^.  In  (9)  we 
have  assumed,  in  order  to  determine  the  velocity  associated 
with  the  Doppler  width,  that  the  ion  kinetic  energy  is  given 
by  the  expression 


17.  =  £^£32  =  „2/322  £  £o 
^0  Ofl  '■q  ’ 


(10) 


a  value  that  can  be  shown  to  be  the  upper  bound  of  possibil¬ 
ities  for  a  uniformly  ionized  spherical  cluster  [35],  Table  1 
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Table  1.  Estimated  quantities  associated  with  the  dynamics  of  multiphoton  excited  clusters  and  soft  X-ray  amplification.  The  dominance  of  Doppler 
broadening  is  assumed  and  n  =  n^..  Sec  text  for  definitions  and  details 
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presents  the  estimated  quantities  for  a  selected  group  of  ma¬ 
terials.  It  is  immediately  seen  that  very  favorable  conditions 
for  amplification  {g^  ~  10^  cm'*)  may  be  achievable  for 
electron  densities  <  10*®  cm'^  at  wavelengths  as  short 
as  =  100  A  with  the  use  of  subpicosecond  excitation. 


4  Concept  for  .Molecular  Design 

A  constraint  governing  the  molecular  design  of  a  system 
optimized  for  X-ray  production  can  be  derived  from  the 
conditions  expressed  in  (I),  (5).  {6)»  and  (7).  In  general,  the 
allowed  zone  for  prompt  X-ray  emission  from  the  j-shell 
of  atoms  in  a  homonuclear  cluster  appears  as  a  triangular 
region  of  the  kind  illustrated  in  Fig.  16.  If  we  identify  the 


log  n 

Fig.  16.  General  representation  of  the  allowed  zone  for  prompt  j -shell 
emission  from  a  homonuclear  cluster  as  a  function  of  intensity  (f) 
and  cluster  size  n.  The  three  constraints  ansing  from  (1).  (6).  and  (7) 
are  shown  along  with  the  venex  3  having  coordinates  (n^,  Z^).  The 
diagram  is  similar  to  that  shown  for  the  Xc  .V-shell  in  Fig.  10 


optimization  of  the, system  with  (1)  the  utilization  of  the 
lowest  possible  intensity  and  (2)  excitation  of  the  cluster 
under  conditions  of  maximal  amplification  per  unit  density, 
then  the  venex  13  defines  that  optimum  subject  to  the 
subsidiary  condition  tiq  =  An  explicit  statement  for  this 
situation  can  be  easily  found.  From  equality  of  (6)  and  (7), 
we  find 


(11) 


an  expression,  which,  when  combined  with  the  condition 
rig  =  with  the  use  of  (5),  gives  the  result 


mc^ 


m) 


=  1. 


(12) 


Finally,  if  the  factor  Z  is  represented  by  setting 
from  (1),  then  the  overall  expression 


1287ra^  \mc^  J  x\hlg  mcP-  V2;rr^y 


=  1 


(13) 


is  found,  a  statement  which  relates  (1)  the  atomic  properties 
(^c  (2)  the  molecular  structure  constant  (Tq).  and  (3) 

the  parameters  describing  the  irradiation  {Iq,  A).  Although 
not  fully  establishing  any  panicular  molecular  configuration. 
(13)  explicitly  illustrates  the  relative  dependencies  of  these 
imponant  variables  with  respect  to  the  optimum  condition. 


5  Conclusions 

A  unified  picture  of  high-intensity  multiphoton  processes, 
involving  atoms,  molecules,  and  solids,  may  be  developed 
through  the  study  of  clusters.  Of  panicular  significance  are 
possible  inua-cluster  processes  that  can  influence  the  level 
of  ionization  produced  and  lead  to  the  excitation  of  inner- 
shell  states  with  the  subsequent  generation  of  prompt  X-rays. 
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The  behavior  of  shape  resonances  [41]  characteristic  of  in¬ 
elastic  electron  scattering  processes,  panicuiarly  for  ionized 
systems  [42]  for  which  rather  little  is  known,  feature  promi¬ 
nently  in  the  dynamics.  Overall,  at  least  six  previously  un¬ 
explained  or  anomalous  experimental  observations  involving 
three  materials  tind  interpretations  consistent  with  the  phys¬ 
ical  picture  developed.  It  is  considered  significant  that  this 
synthesis  can  be  achieved  with  the  exclusive  use  of  known 
atomic  and  molecular  properties. 

The  definition  of  a  critical  cluster  size  =  (r^/Zcr^,)^ 
for  inner-shell  excitation  highlights  the  complexity  of  the 
relationships  governing  the  coupling  by  combining  (1)  a 
static  structural  factor  (fq),  (2)  a  dynamical  consideration 
involving  inelastic  electron  collisional  processes  (0*^^),  and 
(3)  the  nonlinearity  of  the  interaction  with  the  strength  of 
the  radiation  field,  which  is  implicitly  represented  by  the 
factor  Z. 

The  general  result  is  that  molecular  systems  and  clusters 
intrinsically  have  properties  making  them  considerably  more 
suitable  than  atoms  for  the  produc'  on  of  X-ray  emission 
through  the  mechanism  of  multiphoion  coupling  at  high  field 
strengths.  This  outcome  creates  the  possibility  of  designing 
molecules  [43]  optimized  for  efficient  X-ray  production  at 
selected  wavelengths  and  provides  a  general  method  for  the 
generation  of  nonequilibrium  distributions  of  excited  states 
suitable  for  X-ray  amplification.  If  this  can  be  accomplished, 
a  class  of  presently  unknown  molecular  materials  is  likely 
to  emerge.  Presently,  since  the  available  clues  indicate 
that  heavy  atoms  are  strongly  favored,  such  materials  may 
possibly  be  related  to  metallic  clusters  [1],  giant  molecules 
[44]  involving  atoms  such  ,  Au.  Pb,  Th,  and  U,  and  systems 
incorporating  rare  ean^  ...oms  [45,46]  or  transition  metals 
[47]. 
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Muitiphoton  Induced  X-Ray  Emission  from  Kr  Ousters  on  3/-SheII  (  — 100  A) 

and  ^-Shell  ( — 6  A)  Transitions 
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University  of  Illinois  at  Chicago.  845  W.  Taylor.  Room  2136.  Chicago.  Illinois  60607-7059 
(Received  2  September  1993;  revised  manuscript  received  28  December  1993) 

Experiments  demonstrating  the  role  of  cluster  formation  on  multiphoton*induced  x-ray  emission  and 
the  scaling  of  this  phenomenon  into  the  kilovolt  range  have  been  performed  on  Kr.  For  the  Kr  M  shell, 
augmentation  of  Kr,  formation  leads  to  a  large  increase  in  Kr*"^  (4p—  Id)  emission  (— 100  A)  and  the 
appearance  of  a  strong  band  at  —90  A.  The  observation  of  L-shell  transitions  {— 5-7.5  A)  demon¬ 
strates  the  scaling  of  this  phenomenon  into  the  kilovolt  region  and  leads  to  the  conclusion  that  the  in¬ 
teraction  produces  direct  inner-shell  excitation  with  the  emission  of  prompt  x  rays. 

PACS  numbers;  32.80.Wr.  32.30.Rj 


It  has  been  postulated  [|]  that  a  unified  picture  of 
high-intensity  short-pulse  multiphoton  processes,  embrac¬ 
ing  atoms,  molecules,  and  solids,  appears  possible  through 
the  study  of  atomic  and  molecular  clusters.  Of  particular 
significance  in  this  analysis  are  intracluster  processes  in¬ 
volving  inelastic  electron  collisions  of  the  form 

e~+X^e-+e-  +  {X-^)*  (1) 

that  can  lead  to  the  production  of  inner-shell  excited 
species  capable  of  prompt  x-ray  emission. 

Essential  predictions  of  the  proposed  model  have  been 
explicitly  tested  in  new  experimental  studies  involving  five 
;aomic  shells  [Kr(A/),  Kr(I).  Xe(A'),  XefA/).- and 
Xe(Z.)].  This  Letter  establishes  the  crucial  role  of  cluster 
formation  in  the  generation  of  the  anomalous  emissions 
previously  observed  (2-4)  from  the  Kr  M  shell  (  —  110 
eV)  and  demonstrates  the  scaling  of  this  new  phenome¬ 
non  into  the  kilovolt  range  (-800-2100  eV)  through 
measurements  of  Kr  I-shell  radiation. 

The  detection  of  intense  Kr’"*"  (4p— 3d)  emission  at 
—  100  A  in  earlier  spectroscopic  studies  [2.3]  using 
pulsed-gas  targets  presented  an  obvious  paradox.  Mea¬ 
surements  of  ion  production  [5]  using  tenuous  gas  targets, 
conducted  under  identical  conditions  of  irradiation,  had 
demonstrated  that  no  production  of  Kr’"*"  was  possible 
from  free  Kr  atoms.  Hence,  the  identification  [6,7]  of  the 
Kr*'*’  emissions,  produced  by  either  a  prompt  or  a  de¬ 
layed  mechanism,  was  fully  anomalous.  However,  if 
sufficient  cluster  formation  were  produced  in  the  hydro- 
dynamic  flow  of  the  jet  furnishing  the  target  [2-4].  the 
cluster  model  showed  that  this  paradox  could  be  entirely 
resolved;  the  anomalous  emissions  found  a  simple  ex¬ 
planation  which  applied  generally  to  all  cases  considered 
[1]. 

The  multiphoton-induced  production  of  prompt  short 
wavelength  radiation  from  clusters  requires  that  certain 
conditions  be  met  [|].  These  requirements  establish  rela¬ 
tionships  between  the  intensity  of  irradiation  (/)  and  the 
cluster  size  (n)  which  define  allowed  zones  for  the  emis¬ 


sion  to  occur.  Such  emission  can  only  take  place  if 
sufficient  kinetic  energy  is  possessed  by  the  colliding  elec¬ 
tron  in  reaction  (1)  to  excite  an  inner  electron.  This  con¬ 
dition,  based  on  a  simple  formulation  of  the  work  done  by 
the  external  field  on  a  free  electron  in  the  cluster,  leads  to 
the  specification  of  two  limiting  intensities  given  by 


/o 


itcan^^hri 


/,(A.)  — 


2 

mc^ 

c 

n^^ 

a 

X 

X^ 

(2) 


in  which  ro,  a,  X,  m,  X,,  and  c  denote  the  inner-shell 
binding  energy  [8],  the  interatomic  spacing,  the  fine 
structure  constant,  the  wavelength  of  irradiation,  the 
electron  mass,  the  electron  Compton  wavelength,  and  the 
speed  of  light,  respectively.  Hence,  an  inner-shell  j  —  1 
electron  with  binding  energy  f,  can  be  ionized  if  the  in¬ 
tensity  exceeds  a  loy/er  bound  given  by  either  /©  or  /,(X), 
whichever  is  greatest.  In  addition  to  the  production  of  a 
vacancy  in  the  j  —  1  shell,  prompt  y— j  —  1  emission  re¬ 
quires  the  retention  of  at  least  one  electron  in  the  outer  j 
shell  during  the  course  of  irradiation.  This  determines  an 
upper  bound  on  the  intensity  which,  with  the  use  of  the 
tunneling  ionization  model  [9],  leads  to  a  limiting  intensi¬ 
ty  tmMxij)  given  by 


/m.*(y)  - cE/fy )/l  line ‘Z ^  , 


(3) 


where  Epij)  is  the  ionization  potential  of  the  most  tightly 
bound  y-shell  electron,  Z  the  resulting  ionic  charge,  and  e 
the  electronic  charge.  Hence,  the  allowed  zones  for 
prompt  x-ray  emission  from  atoms  of  a  homonuclear  clus¬ 
ter  generally  appear  in  the  form  shown  for  the  Kr  M  and 
L  shells  in  Hg.  I. 

This  picture  leads  to  the  prediction  that  an  increase  in 
the  cluster  density  and  average  cluster  size  should  lead 
to  both  a  corresponding  intensification  and  a  spectral 
modification  of  the  detected  emission.  Since  Kr,  cluster 
formation  can  be  strongly  augmented  by  cooling  the  flow 
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in  the  jet  [10, 1 1],  in  light  of  the  fonn  of  the  allowed  zone 
for  the  M  shell  in  Fig.  1,  a  comparison  of  spectra  ob¬ 
served  at  comparable  pressures  (densities)  and  substan¬ 
tially  different  temperatures  should  indicate  the  impor¬ 
tance  of  the  role  of  cluster  species  in  the  production  of  the 
radiating  ions. 

The  experimental  apparatus  used  has  been  described 
elsewhere  [31.  The  laser  [I2j  employed  in  the  present 
study,  which  had  a  wavelength  of  248  nm,  a  pulse  width 
of  —300  fs,  and  a  power  of  —0.7  TW,  was  focused  with 
an  //lO  CaF:  lens  into  the  gas  target  and  provided  a 
maximum  focal  intensity  of  (0.5-1.0)  x  |0'’  W/cm^  con¬ 
ditions  essentially  the  same  as  in  the  previous  work  [2-4j. 
The  gas  target  was  produced  by  a  high-pressure  pulsed 
valve  fitted  with  a  circular  sonic  nozzle  having  a  diameter 
of  0.5  mm.  The  gas  flow  was  cooled  with  free-flowing  dry 
nitrogen  which  had  passed  through  a  liquid  nitrogen  bath. 

A  thermocouple  attached  to  the  valve  body  measured  the 
temperature  of  the  gas  and  a  grazing  incidence  spectrom¬ 
eter  was  used  to  record  the  spectra  [2-4l. 

The  relevant  spectral  comparison  for  Kr  in  the  72-1 12 
A  region  is  presented  in  Fig.  2,  The  Kr  spectrum  illus¬ 
trated  in  Fig.  2(a)  was  obtained  with  a  stagnation  pres¬ 
sure  of  1 15  psia  at  293  K,  conditions  very  similar  to  those 
of  the  earlier  work  [2,3j.  Figure  2(b)  presents  the  Kr 
spectrum  corresponding  to  a  significantly  reduced  nozzle 
temperature  (238  K)  and  a  stagnation  pressure  (130 
psia).  nearly  the  same  as  that  pertaining  to  the  data  in 
Fig.  2(a). 

Three  principal  differences  distinguish  the  spectra 
shown  in  Figs.  2(a)  and  2(b).  They  are  (I)  a  large  in¬ 
crease  in  the  observed  Kr’"^  (4p— 3<f)  signal  upon 
reduction  of  the  temperature  (e.g.,  from  —30  to  —450  at 
-~I03  A),  (2)  the  appearance  of  a  strong  highly  struc- 
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le  tured  band  at  90  A  at  the  lower  temperature,  and  (3)  the 
)-  large  relative  reduction  in  Fig.  2(b)  of  the  narrow 
I-  features  appearing  in  the  72-83  A  region  of  Fig.  2(a). 

Spectral  data  [6,7]  indicate  that  the  features  observed  in 
e  Fig.  2(b)  in  the  86-97  A  region  arise  from  Kr'®"*" 
^<^)  transitions  in  the  -86-91  A  range  and  Kr’"*" 
I  (4j4/7— 4s3</)  lines  in  the  -90-97  A  interval.  A  con¬ 
tribution  from  3/7  excitation  in  both  Kr*'*’  and  Kr'®'*'  may 
occur  in  this  same  region.  The  spectrum  also  gives  evi¬ 
dence  for  the  presence  of  a  broad  unresolved  band  in  the 
-76-85  A  region  composed  of  dense  arrays  (6J  originat¬ 
ing  from  Kr'®"*"  (4j4p— 4j3</)  and  Kr"'*'  (4p— 3d) 
transitions. 

These  two  signatures,  a  sharply  increased  signal 
strength  and  the  appearance  of  new  categories  of  transi¬ 
tions  and  higher  charge  states,  are  precisely  the  outcomes 
anticipated  from  enhanced  cluster  formation  [ij.  One 
expects  that  a  higher  cluster  (Kr„)  density  for  fixed  n 
would  cause  a  corresponding  increase  in  the  emission  and 
that  the  formation  of  larger  clusters  would  lead  to  a  fur¬ 
ther  intensification  of'the  emission  and  to  new  stages  of 
ionization  and  classes  of  transitions.  Since  the  90  A  band 
only  appeared  in  the  cooled  spectrum,  we  believe  that  it  is 
mainly  associated  with  the  formation  of  the  larger  clus¬ 
ters  produced  at  the  lower  temperature  (238  K). 

Figure  2(a)  also  exhibits  the  clear  visibility  of  relative¬ 
ly  weak  narrow  transitions  in  the  -72-83  A  region 
which  can  be  associated  with  4/  —  3d  transitions  [7]  in 
Kr  .  In  distinct  contrast  to  the  lines  in  the  85-104  A 
range,  these  features  do  not  grow  dramatically  as  the 
temperature  is  decreased.  They  remain  weak  and  become 
relatively  insignificant  as  the  temperature  is  reduced. 
This  sharply  contrasting  behavior  indicates  the  action  of 
a  different  mechanism  for  their  excitation.  Since  the  two 
spectra  shown  in  Fig.  2  were  generated  under  essentially 
the  same  average  plasma  conditions,  conventional  col- 
lisional  mechanisms,  such  as  recombination,  should  occur 
in  largely  the  same  fashion  in  both  cases.  This  expecta¬ 
tion  matches  the  observed  behavior  of  the  transitions  as¬ 
sociated  with  4/  orbital  excitation  and  we  attribute  those 
weak  features  to  recombination. 

The  cluster  model  predicts  the  existence  of  large  al¬ 
lowed  regions  associated  with  level  structures  in  the  kilo¬ 
volt  range,  a  property  manifest  in  Fig.  I  through  compar¬ 
ison  of  the  zones  associated  with  the  M  and  L  shells 
Specifically,  the  model  [ll  predicts  that  the  irradiation  of 
Kr  clusters  with  intensities  above  —10'*  W/cm*.  for 
nearly  an  arbitrary  cluster  distribution,  should  result  in 
emission  of  prompt  x  rays  from  L-shell  transitions  [7 
13.141  in  Kr  ions  (—7  A). 

Experiments  testing  this  scaling  behavior  have  been 
conducted  with  a  peak  intensity  of  — 8x  lO'*  W/cm^.  In 
this  case,  the  focusing  was  provided  by  an  //3  off-axis 
parabolic  mirror  and  a  crystal  spectrometer,  using  Kodak 
Industrex  film  and  PET  for  the  4.5-7.5  range,  was  used. 

The  L-shell  emission  spectrum  of  Kr  shown  in  Fig  3 
which  involves  transitions  [7.13.14J  from  Kr’"*"  (^-24- 
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FIG.  2,  Comparison  of  Kr  spectra  in  the  72-112  A  region 
produced  at  different  temperatures  by  multiphoion  excitation 
with  subpicosecond  irradiation  at  248  nm  with  a  maximum  in¬ 
tensity  of  (0.5-!.0)x|0'’  W/cm^.  The  vertical  scales  in  both 
spectra  arc  absolute  values  giving  a  valid  measure  of  the  com¬ 
parative  signal  strengths.  The  spectral  resolution  is  0.9  ±0.1  A. 
(a)  Nozzle  temperature  293  K  and  stagnation  pressure  115 
r  psia.  UI  designates  unidentified  transitions.  The  base  line  cor¬ 
responds  10—17  counts/shot,  (b)  Nozzle  temperature  238  K 
and  stagnation  pressure  130  psia.  Unresolved  arrays  from 
Kr  and  may  contribute  to  the  signal  in  the  —76-85  A 
region. 


27)  in  the  4.5-7.5  A  region,  was  recorded  with  720  laser 
pulses.  This  spectrum  exhibits  several  salient  properties: 
(Da  group  of  strong  3/—*  21*  lines  in  the  6.3-7.7  A  re¬ 
gion,  (2)  a  small  set  of  far  weaker  4/, 5/—  2p  transitions 
in  the  4.8-5. 6  A  region,  (3)  lines  with  significant  strength 
having  2s  excitation  (e.g.,  F),  and  (4)  a  substantial  num¬ 
ber  of  transitions  involving  satellite  lines  (l4l  [e.g., 
Rk^^^([//2p^3/3/'—  2p^3d3l3r)  and  Kr^^'^(U/2p^3I 
^  2p^3d3l)  shown  in  the  inset]. 

Several  considerations  bear  directly  on  the  interpreta¬ 
tion  of  these  data.  First,  the  radiating  charge  states  are 
anomalous.  In  analogy  with  the  Kr’"**  transitions  dis¬ 
cussed  above,  the  presence  of  lines  from 
“24-27)  is  equally  anomalous,  since  the  tunneling  ion- 
•  izaiion  picture  [5,9]  predicts  no  charge  stales  higher  than 
Kr“  .  Second,  the  radiative  transitions  are  fast.  The 
lifetimes  [15]  of  the  3j— 2p,  3d^2p.  and  3p— 2j 
transitions  occurring  in  the  -6-7  A  region  fail  largely  in 
the  10-1000  fs  range  and  the  2p^  2j(  — 200  eV)  transi¬ 
tion  [7]  has  a  lifetime  [16,17]  less  than  10  ps.  Third, 
recombination  and  plasma  collisional  processes  are  rela¬ 
tively  slow.  Considering  the  properties  of  the  flow,  the 
relative  positions  of  the  focal  zone  and  the  nozzle,  and  the 
level  of  ionization  (Zs25),  the  maximum  average  elec¬ 
tron  density  [4]  is  —  5 x  I0^®  cm  This  gives  a  rate  less 

1812 


Wavdength  (A) 

FIG,  3.  Kr  Z.-shclI  emission  spectrum  produced  by  muliipho- 
ton  excitation  with  248  nm  radiation  at  a  peak  intensity  of 
— 8x|0'*  W/cm^.  The  stagnation  pressure  of  the  pulsed  jet 
was  115  psia  and  the  gas  temperature  was  243  K.  Transitions 
from  Kr^**"  (^—24-27)  are  identified.  The  inset,  presenting 
data  from  another  exposure,  details  the  6,5-7.0  A  region  and 
shows  the  prominent  Kr^^’*'(l^)  and  Kr^‘*'(i/)  satellite  lines  ob¬ 
served  in  these  experiments, 

than  10*®  s"**  for  both  recombination  [18]  and  2s -^2p 
collisional  excitation  [16], 

In  light  of  these  considerations,  the  experimental  obser¬ 
vations  lead  to  the  following  interpretations.  The  five 
weak  transitions  iA^E)  associated  with  Kr^*'*’  ail  involve 
transitions  of  the  form  nl  —  2p  with  n  >  4.  Since  the  .V- 
shell  limit  [/m*x(A^)]  falls  entirely  below  the  allowed  L- 
shell  zone  in  Fig.  I,  the  postulated  cluster  mechanism 
cannot  account  for  the  presence  of  these  lines.  In  analogy 
with  the  weak  4/—  3d  features  occurring  in  Fig.  2(a), 
these  lines  are  attributed  to  recombination  from  a 

species  identified  by  transitions  {HJJ.N.O), 

The  strong  lines  observed  in  the  6.3-7.7  A  region  uni¬ 
formly  involve  transitions  with  3s,  3p,  or  3d  electrons. 
Using  the  most  weakly  bound  of  these  (3d),  and  combin¬ 
ing  the  information  in  Fig.  I  with  the  conditions  for  tun¬ 
neling  ionization  [1,9],  we  can  estimate  the  minimum 
cluster  size  necessary  to  observe  L-shcll  ionization  and 
still  retain  a  3d  electron.  Since  the  maximum  intensity 
consistent  with  the  retention  of  a  3d  electron  is  —  2x  I0** 
W/cm^  the  L-shell  zone  indicates  that  clusters  with 
n>3  would  suffice.  Evidence  [10,11]  on  Kr,,  formation 
indicates  that  an  appreciable  density  with  sizes  at  least  as 
large  as  /r— 10  is  formed  under  our  experimental  condi¬ 
tions. 

Lines  F  and  G  involve  2s  excitation  in  Kr^"*".  Since 
2p^2s  radiative  decay  is  fast  and  electron  collisional 
communication  between  the  2s  and  2p  states  is  slow,  the 
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formation  of  2slp*‘Zp  levels  by  recombination 

from  excited  2s2p^  Kr^^'*'  ions  would  have  to  be  rapid. 
However,  since  the  plasma  conditions  cannot  support  a 
recombination  rate  comparable  to  the  radiative  decay,  it 
follows  that  the  Kr**'*’  excited  states  are  directly  pro¬ 
duced  in  the  excitation  of  the  cluster.  The  same  reason¬ 
ing  and  conclusion  hold  in  relation  to  the  2d  —  Ip  satel¬ 
lite  lines  fe.g.,  Kr^^^lK)  and  Kr^'^{£/)I. 

The  comparison  of  the  spectrum  illustrated  in  Fig.  3 
with  corresponding  data  obtained  with  a  Z  pinch  [13]  re¬ 
veals  a  relative  abundance  of  satellite  lines  in  the  multi¬ 
photon  spectrum  and  significant  reversals  in  relative  line 
intensities.  Both  observations  are  in  natural  accord  with 
the  cluster  model.  The  strong  appearance  of  satellite 
transitions  in  the  multiphoton  spectrum  [e.g.,  Kr^^'*’(K)], 
at  the  relatively  low  average  electron  density  characteris¬ 
tic  of  this  experiment,  is  again  indicative  of  a  direct 
mechanism  of  production  [l]  in  the  cluster  which  can  be 
associated  with  the  existence  of  the  allowed  L-sheil  re¬ 
gion.  In  the  6.3-7.7  A  region,  the  multiphoton  (my) 
3j  —  2p  transitions  consistently  exhibit  the  anomalous 
property  of  having  a  greater  intensity  than  the  corre¬ 
sponding  2d—2p  lines,  while  the  reverse  holds  for 
the  Z-pinch  (Z)  spectrum  [I3J.  Specifically,  for  Kr^*"^ 
R>L  (my).  R<i  (z)  and  for  Kr^^'*'  N>J  (my). 

^  ^Z).  Recombining  collisionally  dominated  optical¬ 

ly  thin  plasmas  characteristically  produce  emission  with 
nd—2p  transitions  dominating  [I9]  the  corresponding 
ns  —  2p  lines. 

The  anomalous  strength  of  the  3j  transitions  congenial¬ 
ly  matches  the  expectations  of  the  cluster  picture.  Elec¬ 
tron  loss  can  occur  by  tunneling  ionization  and  intraclus¬ 
ter  inelastic  electron  collisions.  With  a  rnaximum  intensi¬ 
ty  of  8x  10'*  W/cm’,  3j  electrons  cannot  be  removed  by 
the  tunneling  process  [9],  although  the  full  3</ shell  is  ion¬ 
ized  for  intensities  >  2x  io'»  W/cml  Collisionally.  since 
the  2d  ionization  cross  section  [20.21]  is  approximately 
tenfold  the  value  for  the  3j  level,  retention  of  the  2s  elec¬ 
tron  is  again  favored.  It  follows  that  in  the  region  corre¬ 
sponding  to  L-shell  excitation,  particularly  for  intensities 
above  '-10'*  W/cm^  a  higher  probability  exists  for  re¬ 
tention  of  3j  electrons  than  for  2d.  This  tendency  should 
also  be  enhanced  for  the  higher  charge  states  seen,  as¬ 
suming  a  greater  intensity  for  their  production.  Namely, 
should  be  greater  for  Kr*^'*’  than  for 
Kr*  ,  a  feature  that  is  reflected  in  the  muitiphoton  data 
(.N/J>R/L). 

In  summary,  experimental  studies  of  multiphoton- 
induced  emissions  in  the  5-7.5  A  region  from  Kr  have  es¬ 
tablished  that  the  cluster  mechanism  demonstrated  in  the 
-80- 1 10  eV  range  scales  into  the  kilovolt  zone.  The  ex¬ 
perimental  evidence  supports  the  conclusion  that  muiti¬ 
photon  coupling  to  clusters  can  directly  generate  inner- 
shell  vacancies  and  the  production  of  prompt  x  rays.  Fi¬ 
nally,  recent  measurements  [22]  demonstrating  the  multi¬ 
photon  production  of  Xe(L)  emission  further  corroborate 
this  conclusion. 
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Appendix  C;  "Multiphoton  X-Ray  Generation  from  Clusters:  Analysis  of  Energy  Transfer 
Scaling" 
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ABSTRACT 

A  preliminary  analysis  is  given  concerning  the  scaling  relations  for  multiphoton 
production  of  inner-shell  vacancies  and  prompt  X-ray  emission  from  molecules  and  clusters. 
The  discussion,  which  is  based  on  recent  studies  of  X-ray  emission  (5-16  A)  observed  from  Kr 
and  Xe  clusters,  indicates  that  energy  deposition  rates  exceeding  -  1  W/atom  are  feasible  in 
appropriately  designed  molecules  incorporating  heavy  atoms. 

1.  Introduction 

A  cluster/molecule  mechanism^  has  been  postulated  for  the  multiphoton  production  of 
X-rays.  Of  particular  significance  in  this  analysis  are  intra-cluster  processes  of  the  form 

e-  *  X  —  e-  *  e-  ^  (XT  0) 

which  can  lead  to  the  excitation  of  inner-shell  electrons  and  prompt  X-ray  emission.  This 
concept  represents  the  cluster/ molecule  analog  of  an  interaction  previously  considered  for 
multiphoton  excitation  of  free  atoms”. 

Recent  experiments,^’  speciiically  designed  to  evaluate  certain  essential  predictions  of 
the  cluster  picture,  have  furnished  considerable  supporting  evidence  for  the  main  features  of 
the  proposed  model.  These  studies,  which  involved  Krand  Xe,  (1)  estabUshed  the  crucial  role 
of  cluster  formation^  in  the  generation  of  soft  X-ray  (-  80-150  eV)  emission,  (2)  verified  the 
scaling  of  this  new  phenomenon  into  the  kilovolt  (-  800-2100  eV)  region,**  and  (3)  provided 
important  insights  concerning  the  dvnamical  character  of  the  strong-field  coupling  to  the 
clusters.  Although  these  findings  are  only  the  initial  studies  of  these  interactions  and  much 
further  examination  is  necessary ,  there  currently  exists  a  sufficient  body  of  information  on  which 
to  base  an  analysis  of  the  consequences  of  these  results.  Therefore,  a  preHminary  assessment 
of  certain  key  scaling  relations  is  the  goal  of  the  discussion  presented  below.  Particular 
emphasis  is  given  to  an  exploration  of  the  upper  limits  of  power  and  energy  transfer,  since  these 
limits  will  govern  the  efficacy  of  multiphoton  processes  in  the  generation  of  bright  X-ray 
sources. 
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2.  Discussion 


2.1.  Backgwund 

ennnUnf^n^/^^  ^5  conccnung  limiting  cross  sections  for  multiphoton 

f  intensity  hmit  (>  lO^^  W/cm^)  for  heavy  materials.  TTie  analysis 
incor^rated  assumptions  concerSng  the  electronic  motions 
nmcnm?L^?h  conjastent  with  establishing  an  upper  bound  on  the  cross  section. 

TTte  outcome  was  the  jnediraon  that  an  energy  tracer  rate  in  the  range  of  O.I-LO  W/atom  was 

feasiole  with  an  intensity  of  inadiation  of  —  10^9  mZO  rtr/  -»  •  ,  , 

rsTA  «/  lO  -10‘“  w/cm-.  Tills  enormous  value  for  the 

™  of  tnt.ro  mnsfe  IS  mu^y  comptnblt  a  to  devtioptd  during  tht  o.hrit - - 

inadianon  of  solid  surfeces  at  high  intensity.^  ^  ^uupitasccona 

^  necessary  for  such  a  high  rate  of 

^  ^ms/  the  results  of  cuirmt  e:cpeiimental  studies  indcate  that 

SS^men1??n^-  t  "=  f  mes  of  that  magnitude.  SpeciTicany.  the 

Se  ^  fi^on  that  ^les  them  to 

average  rfMiti'rv  ^  nia^.  Hence,  a  gas  of  these  molecules  responds  as  a  solid  whose 

average  density  can  be  mdependentiy  controfled.  In  such  a  medium,  since  the  highly  ionized 
molecules  subsequently  undergo  a  rapid  coulomb  explosion,  the  system  quickly  convLs  it^f 
me°rg^  combining  a  relatively  low  average  panicle  density  within  unLuaUy  elevated 

2.2.  Conditions  for  duster  Excitation  and  X-Ray  Emission 

Consider  a  spherical  clu^  of  n  identical  atoms.  The  muldphoton  induced  production 
of  EmmEi  shon  waveien^  radiation  from'dusters  requires  that  certain  conditions  be  met.« 
pese  requmsmrau  est^Ush  relationships  between  the  intensity  of  inadiation^  and  the  cluster 

to  occur.  Certainly,  such  emission  can 
only  take  place  if  sufTiaent  lonenc  energy  is  possessed  by  the  coUiding  electron  in  Eq  fl)  to 

«c^e  an  inner  elecutm.  This  condition  leads  to  the  specification  of  two  limiting  intensities 


Sit  an 


(2) 


and 


/,(« ■ 

a  ' 


fr. 


I 


(3) 


m  which  <e.  To.  a,  X,  m,  and  c  denote  the  inner-shell  binding  energy,*  the  interatomic 
spacing  of  the  clus^,  the  fine  structure  constant,  the  wavelength  of  iiSiation,  the  elecmon 
the  el^n  Compton  waveiemgth,  and  the  speed  of  tight,  respectively.  In  this  picture, 
m  mn^-sheil  0-1)  electron  with  binding  energy  can  be  ionized  if  the  intensity  ^ 

‘’y  ^  whichever  is  greatest.  In  addition  to  the  p^uction  of 
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a  vacanq/  in  the  (j-I)-sfaeil,  groniEI  j  *  j-1  emission  requires  the  retention  of  at  least  one 
decOTn  m  the  outer  j-sheH  during  the  course  of  inadiaiion.  Tliis  determines  an  utmer  hnnnri 
on  the  intensity  which,  wi±  the  use  of  the  tunneling  ionization  model,^  to  a  limiting 
mtensiry  I^C)  given  by 


-  cE^oym 

where  is  the  io^don  potendai  '0  of  the  most  tightly  bound  j-shell  electron.  Z  the  resulting 

frnm  arnm?*n^  Charge.  Hcnce,  the  allowed  zone  for  prompt  X-ray  emission 

femaoms  of  a  homonuclear  cluster  appears  geneiaUy  as  a  triangular  region  of  the  form  shown 


(eg  n 


Fig.  l;  GeaenJ  repr««iuuon  of  the  iJIowed  zone  for  prompt  j-shcU  amaion  Horn  m  homonuclear  cluster  as  . 
fMcaoa  of  mtensity  (T)  and  cluster  sue  a.  Tie  three  coasnatais  arisiog  ftom  Eqi.  (2).  (3).  and  (4)  axe 
shown  liong  wiUi  the  corresponding  mtersecuons  P.  Q,  and  J, 

It  can  be  readily  shown  ^  that  an  approximate  estimate  of  the  number  N_  of  excitations 
occurring  through  reacuon  (1)  in  the  cluster  can  be  written  as  the  product  of  the  number  of 
decu-ons  produced  in  the  cluster  (nZ).  the  atomic  density  in  the  duster  the  cross  section 

for  the  inelasnc  channd  under  considemuon  and  the  characteristic  scale  length  R  of  the 
duster  given  approximatdv  by  a  ^  ui  uic 


(5) 


R  - 

Therefore,  for  reaction  (1), 

m  {nDa =  n*^Z^,  (6) 

wi±  z  determined  through  Eq.  (4)  for  the  corresponding  intensity.  Equation  (6)  is  a 
conservative  estimate,  since  it  assumes  that  the  electrons  driven  by  the  Kctemal  field  experience 
only  a  single  transit  through  the  duster  and  are  subsequently  lost  to  ionization. 

An  important  immediate  consequence  of  Eq.  is  the  existence  of  a  critical  size 


defined  by  the  condition  *  n.  Therefore,  if  corresponds  to  an  inner-shell  ionization,  with 
n  S  ng  an  inner-shell  vacancy  would  be  produced  in  every  atom.  Furthermore,  since  the  cridcal 
^e  of  the  duster  depends  inversdy  on  the  square  of  the  atomic  density  ■  r|)  and 
inversdy  on  the  cube  of  the  cross  section  for  exdtation  n^  is  a  direct  analogue  of  the 

nuclear  concept  of  a  cridcal  mass,  with  r^ladng  the  neutron  fission  cross  secdon  and 
Z  playing  the  role  of  the  average  neutron  yield  v.  We  note  that  the  deflnidon  of  a  cridcal 
duster  size  n^  for  inner-shell  exdtadon  highlights  the  complexity  of  the  rdadonships  governing 
the  coupling  by  combining  (1)  a  siadc  structural  factor  (r^),  (2)  a  dynamical  consideradon 
involving  inelasdc  dectron  coilisional  processes  (ffgj),  and  (3)  the  nonlinearity  of  the  interacdon 
with  the  strength  of  the  radiadon  fidd,  which  is  impiicidy  represented  by  the  factor  Z. 

2.3.  Synopsis  o/Experimensal  Findings 

The  zeneral  picture  presented  above  has  been  spedficailv  tested  in  a  recent  series  of 
experiments.^’'*  Specifically,  these  studies  have  examined  the  Xe(N)  and  Kr(M)  shells^  in  the 
80- 150  eV  region  and  the  Xe(M)  and  Kr(L)  shells^  in  the  SOO-2100  eV  range.  The  grecu-oscop- 
ic  results  obtained  in  all  four  cases  were  found  to  be  generally  in  agreement^’**  with  the 
expectadons  of  the  duster  picture.  *  The  studies  of  the  Xe(N)  *and  KrfM)  shell  emissions^ 
established  the  crucial  role  of  duster  formadon  in  the  generadon  of  the  observed  specua.  The 
corresponding  examinadons  of  the  .Xe(M)  and  Kr(L)  shell  radiadon  established  important  aspect 
of  die  scaling  of  this  new  phenomenon  into  the  Jdlovolt  spectral  range. 

3.  PRELIMINARY  ANALYSIS  OF  SCALING  RELATIONSHIPS 

Important  dues  concerning  the  mechanism  of  interaction  and  the  scaling  of  the  energy 
transfer  rate  arose  from  the  analysis^  of  the  specaa  in  the  5-16  A  range  involving  the  Kr  L-shdl 
and  Xe  M-shdl  emissions.  Spedfically,  these  data  provided  detailed  informauon  concerning  the 
dynamics  of  the  energy  transfer  and  indicated  (I)  the  significance  of  muldple  transits  of  the 
driven  electrons  through  the  duster  and  (2)  the  possible  influence  of  coherent  dectron  modons 
on  the  magnitude  of  the  muldphoton  coupling. 
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3.1.  Esiimaxe  of  Multiple  Transia 


The  specirnscopic  evidence'*  indicates  that  fh^  ^  ^  .  u 

interval  which  corresponds  to  the  durarion  of  the  ni,  •  a  confined  to  a  shon  tune 

of  diis  time  from  the  data  gave  a  time  r  -  10  of  the  duster.  An  estimate 

range  anticipated  from  the  dynamics  of  a  couioiS  wiu^  fits  reasonahly  wdl  with  the 

pe.r:od  of  the  incident  radiation  is  siimificantlv  ^  ^ 

external  field  can  experience  than  this  time,  the  dectrons  driven  by  the 

incident  wave.  The  number  m  of  such  transi^"^  by 

2tc 

IT'  (8) 


m 


S.2.  Scales  o/B^ny  Transftr  Sase  ami  PaaUbU  SaU  ^Cohmam  Uowma 

It  was  possible  to  establish  a  lower  bound  on  the  cmerTTnmtnn  j 

rate  per  atom  P„  from  the  Kr  data.**  The 

Kceptionally  high  magnitude.  Moreover,  this  experimental  *“  240  mW/arora,  an 

hish^  than  the  iffiggr  bound  that  could  be  ^Snificantly 

based  on  the  inelastic  scatterine  orocess  «riv«,  ^  ?  the  basts  of  Bq.  (6),  a  statement  that  is 
indtoK  that  U.e  » 

correcuon  to  Eq.  (6)  i,  required.  ^  convqttmtly,  dial  an  iinponanr 

compandon“S^i,  caperi^al  figure  waa  (bund  In 

iraurfer  cro«  recUon 

ihe  maximum  crea  recdon  .r„  war  die  unLrapS  Imre.  Bna  bound  on 

®«  -  87c*f  -  3.6  X  10-“cin^.  (9) 

In  combination  with  the  intensity  used  (L  —  8  x  inI8  w//.«,2\  •  .a.  t  -t  i  a 

value  gave  a  corresponding  thSredcal  ^ure  P  .  °  W/c^)  m‘*«^°volt  studies,'*  this 

which  appropriately  exceeds  the  experimental  v^"e.  ^  mW/atom,  a  magnitude 

mechanism  of  Sianced  cLplL^wSch  L^!rivS*anS^S decSo ^  ^ 

external  field.  Without  this  feature  in^r^^Z  ^i^  eie^mc  moaon^^  mduced  by  the 

in  clear  disagreement  with  the  experimentally  indicat^*4'S%h^h‘*“^'^ 
stems  physically  from  constructive  interfeLce  TL 
scarte.-mg  process  leading  to  atomic  excitation.  Rouehlv 

into  the  analysis  of  the  dusters  leads  to  a  modifican^^i  consideranon 

(10) 

reumare^rtsTiimte  of°Ldti»r  V 

rmreuidd.  for  .V„ 

(8)  and  (10) .  -omoinco  tacior  r  «  2n:Z(X  from  the  product  of  Eqs. 
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(II) 


a  result  which  leads  to  a  correspondingly  changed  cridcai  duster  size 


« I — -7-1  -  r-v. 

2tcZ^o, 


(12) 


-  . .  ^  ^  ^  ^  “  generally  valid  for  the  physical  conditions  of  interest,  the  introduction 

of  this  modificanon  has  four  important  consequences.  They  are  (I)  the  probability  of  ezdtadon 
IS  sufastan^y  raised,  (2)  the  critical  duster  size  is  very  greailv  reduced,  (3)  a  significant 
^veiength  dependent  is  mtroduc^  into  both  of  these  quantities,  iid  (4)  dusters  composed  of 
heavy  atoms  having  high  levels  of  ionization  (Z)  are  potently  favored.  Furthermore,  since  the 
magnitude  of  T  can  be  quite  large,  very  substantial  modifications  of  the  strengths 

expenment^y  produced  can  be  expected.  For  example,  on  the  basis  of  experimental  data 
concerning  Kr  L-shell  exdtation,^  Z  a  25  and  m  a  20,  givi  “ 


case  of  the  Kr  L-shell,  the  ratio  of  Eq.  (6)  to  Eq.  (12)  is 


giving  r  a  500.  Therefore,  in  the 


•  125  X  10*, 


(13) 


an  enormous  reduction  in  the  cridcai  size  which  gives  n'  a  33,  a  value  in  a  range  that  can  be 
pracdcally  achieved. 

5.J.  Molecular  Design  Concept 


The  modificanons  to  the  scaling  reiatianships  discussed  above  alter  the  molecular  design 
condidons  for  X-ray  generation.  A  constraint  governing  the  mofectifar  desiVn  of  a  system 
ggnmjggd  for  x-ray  produedon  can  be  derived  from  the  condidons  expressed  in  Eqs.  (2),  (3) 
(4),  and  (12).  We  recall  that  ^  allowed  zone  for  prompt  X-ray  emission  appears  as  a  triangular 
region  of  the  kind  niusnated  in  Fig.  (1).  If  we  identify  the  opdmizadon  of  the  system  with  (I) 
the  udlizadon  of  the  lowest  possible  intensity  and  (2)  exdtadon  of  the  cluster  under  condidons 
of  mammal  amplificadon  per  unit  density,  ^  then  the  venex  B  defines  that  opdmum  subjea  to  the 
subsidiary  randidon  n^  »  n'g.  An  explicit  statement  for  this  situadon  can  be  easily  found. 
From  equality  of  Eq.  (2)  and  Eq.  (3),  we  find 


an  expression,  which,  when  combined  with  the  condidon  n^  »  n'^  with  the  use  of  Eo  (12) 
gives  the  result  ^ 


f  Hniily,  if  a.e  teor  is  rspitsemsd 

by  semng  *  l^Q)  from  Eq.  (4),  we  obtain 


,  128Ra3 „^2  J  ^  j  |  ^^2  I  j 


1. 


(16) 


an  expression  which  relates  (1)  the  atomic  pnjDcnies  Te  E_m  *  7  m  th^  . 

^r^-s  (r  r)  snd  (3)  the  inffitsity  of 

parncuiar  molecular  configuianon,  Ea.  fl61  ffmltrrtiv  .T-.  ^ _ Aisnmg  any 

these  important  variables  with  respect  of  the  TrlanYg.jcpgndgnciCS  of 

J.^e  ICr  L-Shell  System 

to  the  ““  infontutive  end  led 

to  the  need  for  the  modificauons  expressed  by  the  fector 

IxeZ 
X  ' 


r  ■ 


(17) 


^jmsMces  to  may  allow  ampJificanon.  ^  We  now  explore  the  consequences  of  this 
S!!SI^th“m^n  A  on  3  —  2  tiansiiions  k  Kr  ions  in  ^ 

•Ti  I  assembles  the  parameters  relevant  to  die  exdtatton  of  the  Kr  L-shell  system.  We 

wi^  s^ificaily  consider  the  Kr-6-K  tnnsinon^^’^  ^  lp2iihp^2d7Tt70X 

whicnhasaiad:^ve  lifetimel6of-  10  fs.  An  estimatt  of  tte  Doppler  bxoad^g  to  wo^^ 
occur  as  a  rt^lt  of  the  coulomb  explosion  of  the  cluster  indicates  to  Doppler^broadening 

^  use  a  stimulated  emission  cross  sec^  reflecting  thJ 
^er  wxd±.  Using  the  results  developed  elsewhere,!  and  incorporating  the  enhancement  given 
by  Eq.  (17).  we  can  wnte  the  X-ray  gam  constant  g^  per  unit  dectron  density  n,  as  ^ 


gi 


”8jc~(  e  i  .2  I  i  J“3~]  ^1 

\  •  J  /  \  •/ 


(18) 


fw  homogeneous  clusters  composed  of  atoms  having  mass  M  and  an  X-ray  transition  with  a 

w  indeoendenr  of  the  level  of  atomic 
S^e  '  wavelength  of  exdmion  X.  If  we  assume  further 

that  ±e  greatest  penn^le  dectron  density  of  the  medium  is  the  critical  dectron  densiw 
cor^ponam^  to  wavdength  X.  dien  the  maximum  value  of  the  gain  constant  g,  adi^k 
scales  as  X  ,  a  dependence  heavay  favoring  shorter  wavelengths.  ^ 


Table  I: 


PMmetM  to  the  Kr  L-«heU  Kt2fi+  {lr2j?2p^ - l*=2p22p^3d)  tnnsitioa  u  6.7 

with  excmnaii  u  248  ma.  Tie  coinnm  (fc/tO  aasasae*  thax  the  fall  ionxc  populeoan  tesides 
the  ^<*»^*<* 


Mucriai  j  r,(A)  Tnnitlsa 

1  UM 

«e 

(tV> 

I*„rW/eBrj 

1  rt&) 

"*  1  'S 

1  <» 

f4/n,(enr) 

Kr  1  j 

-  IfiTO 

26  1 

-  •  X  10‘* 

1  4  J  X  10*“  1  -  10 

1  Op)  1 

-  33  1  -  10‘< 

2.6  X  10*“ 

4.  Conclusions 

Recent  experimental  findings  have  indicatrd  that  Kr  and  Xe  dusters  underEoine 
mulUpnoton  examon  in  strong  fields  produce  intense  X-xay  emission  from  directly  erntfd  ionic 
mtes.  An  analysis  of  the  sca^g  rclarionsfaip  involved  in  the  muUiphoton  coupling,  in  light  of 
these  new  expenmmtd  indicate  that  energy  deposition  rates  exceeding  -  1  W/atom  are 
feasible  in  appropriately  designed  molecules  incorporating  heavy  atoms. 
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Multiphoton-induced  X-ray 
emission  at  4-5  keV  from 
Xe  atoms  with  muitipie 
core  vacancies 

A.  McPherson,  B.  D.  Thompson,  A.  B.  Borisov, 

K.  Boyer  ft  C.  K.  Rhodes 

Department  of  Physics  (M/C  273),  University  of  Illinois  at  Chicago. 

845  West  Taylor  Street  Chicago,  Illinois  60607-7059.  USA 

Several  recent  experimentaJ  findings'^  have  pointed  to  a  possible 
route  for  making  an  X-ray  laser,  which  could  in  principle  provide 
an  imaging  system  capable  of  molecular  resolution^.  The  method 
involves  the  multiphoton  excitation  of  atoms  in  van  der  Waals 
clusters  or  in  molecules  to  yield  ions  with  core-electron 
vacancies'*^,  which  can  then  decay  by  emission  of  X-rays,  in  con¬ 
junction  with  a  self-chanelling  propagation  mode  of  electromag¬ 
netic  radiation^  The  multiphoton  excitation  may  be  stimulated  by 
ultrahigh-brightness,  subpicosecond  pulses  of  laser  light^.  We  have 
previously  observed^  emission  of  X-rays  from  L-shell  transitions 
in  core-excited  krypton  atoms  using  this  approach.  Here  we  report 
the  multiphoton  production  of  X-rays  of  wavelength  2-3  A  from 
highly  ionized  xenon  atoms  which  possess  a  large  number  of  inner- 
shell  vacancies  whUe  retaining  several  electrons  in  relatively 
weakly  bound  outer  orbitals.  Atoms  with  this  inverted'  electronic 
configuration  are  designated  ‘hollow  atoms’^’.  We  find  that  gener¬ 
ation  of  hollow  atoms  can  become  the  dominant  excitation  mo^ 


for  such  systems,  making  their  exploitation  in  an  X-ray  laser  a 
real  possibility. 

In  the  cluster/molecule  approach  proposed  previously  for 
X-ray  amplification,  a  molecule  or  a  loosely  bound  cluster  of 
atoms  is  irradiated  by  an  intense  (>10‘^  W  cm“‘)  source  with  a 
wavelength  matched  to  the  excitation  conditions  of  the  atomic 
constituents.  An  initially  ionized  free  electron  in  the  cluster  or 
molecule  may  then  lead  to  further  ionization  via  intra-cluster 
inelastic  electron-atom  collisions: 

By  adjusting  the  excitation  intensity  depending  on  the  inter¬ 
atomic  spacings  and  inner-shell  binding  energies,  processes  such 
as  these  can  lead  to  the  selective  generation  of  inner-shell  excited 
species.  We  have  tested  some  of  the  predictions  of  this  model 
previously  for  the  production  of  Kr(L)  {5-7-A)  emission^  A 
significant  feature  of  the  Kr(L)  spectrum  was  the  anomalous 
strength  of  the  satellite  lines,  a  property  that  indicates  the  occur¬ 
rence  of  inner-shell  (rt  =  2)  excitation  while  outer-shell  (n  =  3) 
electrons  remain  bound  to  the  atom. 

We  have  now  obtained  exp)erimental  evidence  for  the  direct 
multiphoton  production  in  gas-phase  xenon  clusters  of  hollow 
atoms  having  multiple  L-shell  vacancies.  Specifically,  we  have 
observed  Xe(L)  emission  in  the  2-3  A  range  produced  by  excita¬ 
tion  of  the  clusters  with  ultra-intense  subpicosecond  ultraviolet 
radiation.  One  measurement  involves  the  recording  of  X-ray 
‘pinhole’  CCD  (charge-coupled  device)  camera  images  of  the 
radiating  zone  with  a  spatial  resolution  of  -^100  pm.  Another 
entails  spectral  measurements  obtained  with  a  von  Hamos 
curved  crystal  spectrometer®.  Both  studies  were  performed  with 
the  apparatus  and  methods  previously  described^.  The  clusters 

Laser 


Propagation  (248  nm !  270  fe) 


Rayleigh  range  (28J  pm) 

hJ _ I _ I _ I _ I — I — i — > — I _ »— J _ « _ I _ I _ L 

13  14  13  12  11  10  9  8  7  6  5  4  3  2  1 

Z  (100  pm) 


L 


Approximate  position 
of  focal  zone 


_L 

0 


FIG.  1  X-ray  Images  produced  by  a  single  exposure  with  the  double 
pinhole  camera.  The  upper  and  lower  pinholes  correspond  to  the  Xe(L) 
and  Xe(M)  emissions,  respectively.  The  pinholes  had  a  diameter 
^100  pm.  limiting  the  spatial  resolution  to  that  value.  The  stagnation 
pressure.of  the  pulsed  jet  was  75  pounds  per  square  inch  above  atmos¬ 
pheric  pressure  and  the  temperature  of  the  nozzle  was  233  K.  The 
approximate  location  of  the  focal  region  is  indicated  along  with  the 
corresponding  Rayleigh  range  of  the  focusing  system  ( f/3).  The  more 
restricted  spatial  distribution  of  the  shorter  wavelength  emission 
(^fi>x>3.5  keV)  is  expected,  as  significantly  higher  intensities  are 
needed  for  the  L-shell  excitation  relative  to  those  required  for  the  corre¬ 
sponding  M-shell  excitation^ Intensity  grey  scale  is  presented  in 
pseudocolour. 
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FIG.  2  a.  The  spectrum  of  the  short*wavelength  component  [fKs)x> 
3.5  keV)  shown  In  Fig.  1,  recorded  with  720  laser  pulses  by  a  von 
Hamos  spectrometer  operating  in  third  order.  The  spectrum  shows  a 
single  broad  feature  having  a  width  of  ^700  eV  and  an  asymmetrical 
wtng  extending  to  the  blue  side.  Additional  features  are  also  indicated 
and  identified.  Although  occurring  generally  in  the  expected  spectral 
region  for  L-shell  transitions  in  Xe,  the  measured  spectrum  is  entirely 
different  from  that  normally  anticipated  (for  example,  Ne-like).  Specifi¬ 
cally,  the  customary  L-shell  spectrum  at  a  low  average  plasma  density 
(He  ^10^°  cm is  dominated  by  a  few  sharp  lines.  In  contrast.  It  is 
evident  that  the  multiphoton-produced  spectrum  illustrated  is  predomi¬ 
nantly  composed  of  a  very  broad  single  feature.  The  film  was  read  with 
a  film  scanning  system  (Materials  Data  Inc.,  Livermore.  California)  using 
an  Agfa  (Studio  Scan)  scanner  for  negatives.  D  denotes  a  localized 
defect  in  the  film.  The  estimated  accuracy  of  the  abscissa  is  ±50  eV. 
b,  Xe{L)  emission  recorded  in  the  4-6  keV  range  following  the  impact 
of  slow  Xe'*®"  ions  on  Cu  surfaces  as  described  In  ref.  7.  The  features  in 
the  4-5  and  5-6  keV  regions  arise  from  3d-^2p  and  4d-*2p  transitions, 
respectively,  in  an  atom  that  possesses  a  large  number  of  n^3  elec¬ 
trons.  On  account  of  the  interactions  associated  with  the  presence  of 
the  outer  (n^3)  electrons,  the  allowed  transitions  become  distributed 
over  a  large  spectral  region  and  a  broad  band  of  emission  develops. 
The  correspondence  of  the  3d-»2p  emission  in  the  4-5  keV  range  with 
the  multiphoton-produced  spectrum  shown  in  a  is  marked.  The  abscissa 
has  been  adjusted  to  closely  match  that  of  a.  The  estimated  accuracy 
of  the  abscissa  is  ±50  eV.  (Shown  with  permission’’.)  c.  Corresponding 
Xe(L)  emission  recorded  following  the  impact  of  slow  Xe"*®"  ions  on  a 
Cu  surface’  which  also  manifests  properties  similar  to  the  spectrum  in 
a.  The  2p^  and  2p^  (shaded)  zones  overlap  with  a  shoulder  at  ^-5  keV 


were  produced  in  a  cooled  pulsed-gas  jet.  We  use  a  laser’  with 
a  wavelength  of  248  nm,  a  pulse  length  of  ^300  femtoseconds 
( fs,  s)  and  a  peak  power  of  ^0.7  lerawatts  (TW,  10'“  W). 
The  focusing  system  produces  a  maximum  intensity  of  ^8 
exawattscm"’  {few,  10‘®  W). 

The  X-ray  ‘pinhole’  camera  pictures  are  single-shot  exposures 
obtained  by  viewing,  in  a  direction  transverse  to  the  propagation 
of  the  248-nm  beam,  the  plasma  excited  by  the  focused  subpico¬ 
second  248-nm  radiation  directly  below  the  pulsed  jet.  The 
camera  has  two  pinholes  that  can  record  images  simultaneously. 
One  pinhole  has  a  stainless-steel  filter  with  a  thickness  of 
12.7  pm,  making  it  entirely  ‘blind’  to  radiation  with  a  quantum 
energy  ^a)x<3.5keV.  The  other  pinhole  has  a  filter  with  the 
composition  (2,000  A  A1+  1  pm  polycarbonate  +  10  pm  Be)  so 
that  all  radiation  with  a  quantum  energy  fio)\>  700  eV  is  trans¬ 
mitted.  Furthermore,  as  the  Xe(M)  emission  at  keV  was 
found  experimentally  to  be  considerably  brighter  than  the  Xe(L) 
radiation  at  -4  keV,  with  the  use  of  these  filters,  one  image 
represents  the  spatial  distribution  of  the  Xe(M)  radiation  while 
the  other  presents  the  corresponding  distribution  of  the  Xe(L) 
emission.  The  experimentally  observed  spatial  distributions  are 
shown  in  Fig.  1. 

In  Fig,  2  we  show  the  spectra  in  the  4-5  keV  range  obtained 
by  our  method  ( Fig.  2a),  and  those  for  Xe(  L )  emission  following 
the  impact  of  Xe'*^^  (Fig.  2b)  and  Xe'**"  (Fig.  2c)  ions  on  a 
copper  surface’,  which  were  shown  to  come  from  hollow-atom 


E(eV) 


Energy  (eV) 

and  feature  S  in  a.  The  estimated  accuracy  at  the  abscissa  is  ±50  eV. 
(Shown  with  permission’). 


states.  The  spectra  share  several  characteristics  including  general 
spectral  position,  breadth  and  particular  features.  The  overall 
correspondence  leads  to  the  conclusion  that  the  2d-*2p  emis¬ 
sions  observed  in  both  cases  arise  from  the  same  type  of  ionic 
states,  namely  hollow-atom  configurations  in  which  inner-shell 
vacancies  are  associated  with  the  presence  of  a  large  number 
(-10)  of  outer-shell  (n^3)  electrons. 

The  multiphoton  spectrum  shown  in  Fig.  2a  has  a  clear  asym¬ 
metry,  weak  narrow  components  superimposed  on  the  dominant 
broad  feature,  and  visible  modulations  in  the  intensity  that  can 
be  correlated  with  spectral  characteristics  observed  in  studies 
of  ion-surface  emission®*’,  laser-imploded  microballoons’  and 
radiatively  driven  capsules'*^.  The  asymmetry  of  the  multiphoton 
spectrum  can  be  explained  by  the  presence  of  multiple  L-shell 
vacancies.  For  example,  the  broad  feature  S  in  Fig.  2^  at  ^5  keV 
lies  at  a  position  representing  a  null  in  Fig.  2b,  but  this  region 
corresponds  to  a  significant  shoulder  in  Fig.  2c  which  overlaps 
the  multiple-vacancy  2p‘  and  2p^  zones  (shaded).  Three  transi¬ 
tions  in  Xe"*^"  (2.72-2.75  A),  whose  known  wavelengths  are  indi¬ 
cated  to  coincide  with  the  feature  B,  correspond  to  the  strongest 
lines  observed  in  the  microballoon’  and  capsule studies.  Addi¬ 
tionally,  two  narrow  transitions  are  identified  in  the  multiphoton 
spectrum  at  2.94  and  2.86  A,  attributed  to  3s-*2p  transitions  in 
Xe**^"  and  Xe'**’^,  respectively.  The  clear  appearance  of  the  Ne- 
41  like  Xe**^"  (35 -►2/7)  line  mirrors  the  observation  of  the  corre- 
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Charge  state 

FIG.  3  Threshold  intensities  for  ionization  of  electrons  from  bound 
states  of  xenon  as  derived  from  ref.  12.  Specific  electron  shells  and 
levels  are  indicated.  With  a  maximum  experimental  intensity 
/e.p-S  ^  10^®  W  cm'^.  the  Coulomb-suppression  picture  leads  simulta¬ 
neously  to  the  disappearance  of  4cf  states  and  the  retention  of  the  3d 
electrons.  Hence,  the  absence  in  the  multiphoton  spectrum  illustrated 
in  Fig.  2a  of  the  4cf— 2p  band  shown  in  Fig,  2b  and  c  is  another  feature 
of  experimental  agreement  between  the  results  of  the  multiphoton  exci¬ 
tation.  the  findings  of  the  ion-surface  studies^  and  the  expectations  of 
the  cluster  approach^'^. 

spending  Ne-like  (35 -♦2/7)  transition  as  the  strongest 

feature  in  multiphoton  produced  Kr{L)  emission^ 

The  presence  of  the  Xe'^^'^  species  with  the  excited-state  con¬ 
figuration  (2/7'*35)  and  the  overlap  of  the  feature  S  with  the  spec¬ 
trum  shown  in  Fig.  2c  arising  from  the  Xe"^*^  ions  confirm  the 
direct  production  of  double  vacancies  by  the  multiphoton  cou¬ 
pling  and  suggest  that  higher  levels  of  2p  shell  ionization  are 
also  present.  The  generation  of  the  Xe'*^’"  {2p'*2s)  stale  from  the 
neutral  atom  requires  an  energy  transfer"  exceeding  70  keV. 

The  mulliphoton  spectrum  in  Fig.  2a  does,  however,  exhibit 
a  significant  difference  from  those  produced  by  ion-surface 
collisions';  the  feature  at  ^5.5  keV  in  Fig.  2b  and  c,  which  corre¬ 
sponds  to  the  4d-^2p  component,  is  absent  in  Fig.  2a.  This  fact 


FIG.  4  Comparison  of  the  morphologies  of  the  multiphoton-produced 
emission  for  Kr(L)  (ref.  2)  and  Xe(L).  The  two  spectra  are  registered  by 
the  corresponding  (2p*-3s)  Ne-like  transitions  for  Kr^®"  and  Xe^***.  The 
conversion  from  a  line-dominated  spectrum  for  Kr  to  a  continuous  dis- 
tnbution  reflecting  the  characteristics  of  the  hollow-atom  configuration^2 
produced  with  Xe  is  clear. 


is  easily  explained  by  the  Coulomb-suppression  model  of  multi¬ 
photon  ionization'^,  the  threshold  conditions  for  which  are  illus¬ 
trated  in  Fig.  3. 

The  comparison  of  the  morphologies  of  the  multiphoton- 
produced  Kr(L)  and  Xe(L)  spectra  shown  in  Fig.  4  is  infor¬ 
mative.  If  the  spectra  are  registered  to  the  corresponding  2p^3s 
Ne-like  Kr”®"^  and  Xe^*^"^  features,  the  two  spectra  are  essentially 
isomorphic  with  the  major  exception  that  the  gaps  between  the 
dominant  sharp  features  of  the  Kr(L)  distribution  are  filled  in 
by  the  high  density  of  satellite  lines  seen  in  the  Xe{L)  emission. 
Hence  the  line  spectrum  for  charge  state  Z=  36  is  replaced  at 
Z  =  54  by  satellite  transitions  characteristic  of  the  hollow-atom 
spectrum.  Extrapolation  of  this  trend  leads  to  the  conclusion 
that,  for  heavy  atoms  in  clusters  and  appropriately  designed 
molecules,  the  multiphoton  excitation  of  hollow-atom  states  will 
be  the  predominant  excitation. 

A  comparison  of  the  spectrum  shown  in  Fig.  2a  with  a  spec¬ 
trum  thermally  produced  with  an  electron  temperature 
Tc^lA  keV  at  a  comparable  electron  density  {n^  %  10^"^  cm~^),  as 
shown  in  Fig.  5,  reveals  a  significant  difference.  The  thermal 
spectrum  deviates  grossly  from  the  multiphoton-generated  distri¬ 
bution.  Specifically,  the  former  is  essentially  null  between  2.8 
and  3.0  A,  the  region  corresponding  to  the  strong  maximum  in 
the  latter.  Conversely,  the  high  degree  of  spectral  correspondence 
shown  between  Fig.  2a,  b  and  c,  as  the  effective  electron  tempera¬ 
ture  associated  with  the  ion-surface  interaction  is  very  low 
( 5  eV),  indicates  that  the  multiphoton-produced  emission 
involves  a  negligible  contribution  from  thermal  excitation. 
Because  the  multiphoton  spectrum  must  involve  electron  colli- 


FIG.  5  Comparison  of  the  multiphoton-generated  Xe<L)  spectrum  (a) 
with  a  thermally  produced  spectrum^®  (b)  obtained  at  a  comparable 
electron  density  (n,  a:  10^^  cm"®).  The  thermal  spectrum  was  obtained 
at  a  temperature  estimated^®  to  be  1.4  keV.  In  contrast  to  the  spec¬ 
tral  similarity  exhibited  in  Fig.  2,  these  two  spectra  correspond  poorly, 
particularly  in  the  important  region  where  the  maximum  of  the  multi¬ 
photon  spectrum  occurs  (^2.8--3.0  A),  (b  shown  with  permission^®.)  S 
is  a  significant  shoulder  feature.  D  is  a  localized  film  defect,  B  is  a 
strong  line  (see  text),  a  is  an  identified  strong  feature  of  the  thermal 
spectrum^®  and  ^  is  a  feature  of  the  thermal  spectrum  corresponding 
to  B. 
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sions  in  the  kilovolt  range',  this  signifies  that  electron  motions 
produced  by  the  external  field  occupy  only  a  very  small  fraction 
of  the  available  phase  space.  Therefore,  we  can  conclude  that  the 
electronic  motions  causing  the  observed  multiphoton  excitation 
possess  a  very  signficant  level  of  order  and  are  not  appreciably 
perturbed  by  a  thermal  component'^  □ 
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A  simple  derivation  in  the  Coulomb  gauge  of  the  nonlinear  Schrddinger  equation  describing  propaga¬ 
tion  of  powerful  ultrashort  circularly  polarized  laser  pulses  in  underdense  cold  inhomogeneous  plasmais 
is  presented.  Numerical  solutions  are  given  for  the  two-dimensional  axisymmetric  case  for  both  initially 
homogeneous  plasmas  and  static  preformed  plasma  columns.  These  solutions  account  for  (i)  diffraction, 

(ii)  refraction  arising  from  variations  in  the  refractive  index  due  to  the  spatial  profile  of  the  electron  den¬ 
sity  distribution,  (iii)  the  relativistic  electronic  mass  shift,  and  (iv)  the  charge  displacement  resulting  from 
the  transverse  ponderomotive  force.  The  most  important  spatial  modes  of  propagation  corresponding  to 
(1)  purely  relativistic  focusing  and  (2)  the  combined  action  of  both  the  relativistic  and  charge- 
displacement  mechanisms  are  described.  The  latter  leads  to  the  formation  of  stable  confined  modes  of 
propagation  having  paraxially  localized  regions  of  high  intensity  and  corresponding  paraxially  situated 
cavitaied  channels  in  the  electron  density.  It  is  further  demonstrated  that  the  dynamical  solutions  of  the 
propagation  tend  asymptotically  to  the  lowest  eigenmodes  of  the  governing  nonlinear  Schrddinger  equa¬ 
tion.  Finally,  the  calculations  illustrate  the  dynamics  of  the  propagation  and  show  that  the  relativistic 
mechanism  promotes  the  initial  concentration  of  the  radiative  energy  and  that  the  subsequent  charge 
displacement  stabilizes  this  confinement  and  produces  waveguidelike  channels. 

PACS  number(s):  52.40.Nk,  42.65.Jx,  52.35.Mw,  52.40.Db 


1.  INTRODUCTION 

The  interaction  of  relativistically  intense  subpi¬ 
cosecond  laser  pulses  with  gaseous  media  has  been  an 
area  of  vigorous  research  for  the  past  several  years.  For 
ultraviolet  wavelengths  on  the  order  of  200-300  nm,  the 
intensity  region  of  interest,  in  which  relativistic  effects  be¬ 
come  important,  lies  above  ~  10*®  W/cm^.  The  propaga¬ 
tion  of  radiation  in  such  media,  for  intensities  greater 
than  ^10*^  W/cm^,  naturally  causes  strong  nonlinear 
ionization  in  ail  matter.  Hence,  the  pulse  itself,  even  in 
regions  where  the  intensity  is  relatively  low  compared  to 
the  peak  value,  removes  many  electrons  [1,2]  from  the 
atomic  or  molecular  constituents  creating  a  plasma 
column  in  which  the  main  high-intensity  component  of 
the  pulse  propagates.  Therefore,  in  a  reasonable  first  ap¬ 
proximation,  the  investigation  of  the  resulting  propaga¬ 
tion  can  be  divided  into  two  separate  and  distinct  areas. 
They  are  (1)  the  atomic  and  plasma  physics  occurring  in 


the  field  of  the  intense  electromagnetic  wave  leading  to 
the  ionization  and  (2)  the  subsequent  nonlinear  propaga¬ 
tion  of  the  radiation  in  the  plasma  that  is  generated.  The 
work  described  below  concerns  the  latter  issue. 

To  our  knowledge,  Akhiezer  and  Polovin  published  the 
first  treatment  of  high-intensity  electromagnetic  waves  in 
a  cold  plasma  [3].  They  derived  the  equations  describing 
the  propagation  as  a  function  of  the  single  canonical  ar¬ 
gument  icot—kz)  appropriate  for  one-dimensional  wave 
motion,  reduced  the  problem  to  Lagrangian  form  with 
two  integrals  of  motion,  and  presented  either  exact  or  ap¬ 
proximate  solutions  corresponding  to  particular  sim¬ 
plified  special  cases. 

Several  subsequent  treatments  have  been  devoted  to 
the  acceleration  of  charged  particles  in  either  relativistic 
beat  waves  or  in  the  tail  of  a  single  relativistic  pulse. 
Specifically,  Noble  [4]  applied  the  equations  derived  by 
Akhiezer  and  Polovin  to  the  study  of  single-  and  double¬ 
wave  propagation.  In  other  work  [5]  beat  waves  in  hot 
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plasmas  are  described  kinetically.  The  excitation  of  par¬ 
ticles  in  cold  plasmas  occurring  in  the  tail  of  a  single  in¬ 
tense  laser  pulse  has  also  been  investigated  [6]. 

Importantly,  this  previous  work  [6],  which  considered 
the  case  of  linear  polarization,  established  that  harmonics 
of  the  fundamental  frequency  are  not  produced  at  a 
significant  power,  if  the  plasma  is  underdense  and  the 
group  velocity  of  light  is  close  to  the  speed  of  light. 

Utilizing  the  quasistatic  approximation  for  the  fluid 
equations  of  cold  underdense  plasmas,  Sprangle,  Esarey, 
and  Ting  [7]  derived  the  first  nonlinear,  fully  self- 
consistent  set  of  equations  describing  the  propagation  of 
relativistic  laser  pulses  into  plasmas  for  a  one¬ 
dimensional  geometry.  This  formulation  was  used  to  ob¬ 
tain  insight  into  relativistic  self-focusing,  including  the 
self-consistent  electron-density  profiles,  wake-field  gen¬ 
eration,  optical  guiding,  and  second-harmonic  genera¬ 
tion.  However,  since  this  treatment  was  one-dimensional, 
the  multidimensional  case  remained  unsolved. 

An  important  finding  of  Sprangle,  Esarey,  and  Ting  [7] 
concerned  the  diffractional  erosion  of  the  leading  edge  of 
a  pulse  propagating  in  the  plasma.  However,  in  the  work 
described  below,  we  concentrate  on  the  evolution  of  the 
central  portion  of  the  pulse  as  it  propagates  in  the  plasma 
and  neglect  the  associated  erosion  of  the  leading  edge. 

Extensive  literature  exists  on  the  motion  of  electrons  in 
radiation  fields  of  certain  explicitly  given  forms.  The 
solution  for  the  case  of  plane  monochromatic  waves  is 
known  [8],  and  this  result  has  been  considerably  extended 
in  subsequent  work  [9,10].  Of  course,  the  well- 
established  Volkov  solution  for  the  Dirac  equation  also 
exists  [11].  In  particular,  Bardsley,  Penetrante,  and  Mit- 
tleman  [12]  have  numerically  simulated  the  relativistic 
dynamics  of  electrons  in  a  one-electron  picture  that  in¬ 
cludes  the  effects  of  space  charge  and  the  spatial  distribu¬ 
tion  of  the  radiation  field. 

The  first  treatment  of  the  relativistic  self-focusing  in 
plasmas  was  developed  by  Max,  Arons,  and  Langdon 
[13].  In  addition,  the  general  character  of  the  elec¬ 
tromagnetic  propagation  in  plasmas  has  undergone  con¬ 
siderable  analysis.  Schmidt  and  Horton  [14],  Hora  [15], 
and  Sprangle  et  al  [16,17]  have  evaluated  the  thresholds 
for  relativistic  self-focusing  using  analytic  methods. 
Especially  germane  is  the  work  of  Sun  et  al.  [18],  who, 
for  initially  homogeneous  plasmas,  derived  the  two- 
dimensional  (r,z)  nonlinear  Schrodinger  equation  govern¬ 
ing  propagation,  including  consideration  of  the  combined 
effect  of  the  relativistic  nonlinearity  and  charge  displace¬ 
ment.  This  work  presented  the  lowest  eigenmode  of  the 
nonlinear  Schrodinger  equation  and  included  numerical 
evaluation  of  the  threshold  of  relativistic  self-focusing, 
the  value  of  which  has  been  approximately  estimated  in 
earlier  work  [14~17].  Finally,  this  analysis  [18]  presented 
the  (r,2)  dynamics  of  the  propagation  for  the  perturbed 
lowest  eigenmode  of  the  nonlinear  Schrodinger  equation 
for  cases  not  involving  spatial  cavitation  of  the  electron 
density.  We  note  that  Kurki-Suonio,  Morrison,  and  Taji- 
ma  [19]  have  also  developed  the  stationary  analytic  solu¬ 
tions  to  this  equation  for  a  one-dimensional  geometry. 

Additional  related  works  can  be  cited  in  this  context. 
Relativistic  self-focusing  and  beat-wave  phase-velocity 


control  in  plasma  accelerators  are  kindred  subjects  [20]. 
Computations  involving  particle  simulations  revealing 
the  initial  process  of  self-focusing  and  subsequent  pon- 
deromotively  driven  electron  motions  have  been  per¬ 
formed  [21].  Other  calculations  analyzing  the  plasma  dy¬ 
namics  and  self-focusing  in  heat-wave  accelerators  [22], 
as  well  as  the  consideration  of  the  nonlinear  focusing  of 
coupled  waves  [23],  have  also  appeared. 

In  prior  publications,  we  have  (1)  investigated  the  gen¬ 
eral  behavior  of  two-dimensional  (r,z)  axisymmetric  rela¬ 
tivistic  self-focusing,  (2)  presented  results  on  the  process 
of  stabilization  of  laser  pulses  in  plasma  columns  [24],  (3) 
described,  with  the  use  of  an  analytical  model  [25],  the 
steady-state  characteristics  in  cavitated  channels  having 
overdense  walls,  (4)  reported  preliminary  results  of  calcu¬ 
lations  that  evaluated  the  combined  action  of  the  relativ¬ 
istic  and  charge-displacement  mechanisms  and  indicated 
the  formation  of  stable  confined  cavitated  modes  of  prop¬ 
agation  [26],  and  (5)  presented  the  experimental  evidence 
of  relativistic  and  charge-displacement  self-channeling  of 
intense  subpicosecond  ultraviolet  radiation  in  plasmas 
[27],  including  specific  comparison  with  the  results  of 
this  computational  model  [27]. 

The  current  work  presents  a  full  description  of  the 
theory  of  nonlinear  propagation  of  intense  axisymmetric 
ultrashort  laser  pulses  in  cold  underdense  plasmas.  In 
this  analysis  we  use  the  term  ultrashort  to  indicate  that 
the  duration  of  the  pulses  r  satisfies  the  inequality 
r,  »r»r^,  with  r,  and  designating  the  response 
times  of  the  ions  and  electrons,  respectively.  We  consider 
both  homogeneous  plasmas  and  preformed  plasma 
columns.  In  this  study,  no  attempt  is  made  to  establish 
consistency  between  the  local  ionization  state  and  the  lo¬ 
cal  laser  intensity,  the  issue  outlined  in  Sec.  1.  For 
sufflciently-low-Z  materials,  such  as  hydrogen  (H2),  this 
calculation  would  be  unimportant,  since  full  (maximal) 
ionization  would  be  achieved  even  in  rather  low  intensity 
(-lO^^-lO*^  W/cm^)  regions  [28].  Since  the  main  re¬ 
gime  of  interest  for  this  work  involves  intensities  greater 
than  - 10^®  W/cm^  we  believe  that  the  plasma  condi¬ 
tions  we  have  chosen  for  analysis  are  sufficiently  close  to 
the  true  self-consistent  state  to  be  adequate  for  the  in¬ 
tended  scope  of  this  study. 

The  calculations  discussed  below  have  been  performed 
with  the  specific  goal  of  exploring  the  characteristic  dy¬ 
namics  and  stability  of  the  propagation,  including  partic¬ 
ularly,  the  interplay  of  the  relativistic  and  charge- 
displacement  processes.  Therefore,  in  order  to  illustrate 
this  behavior,  numerical  results  are  presented  that  por¬ 
tray  the  propagation  and  self-focusing  action  as  a  conse¬ 
quence  of  (1)  the  purely  relativistic  nonlinearity  and  (2) 
the  combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms.  Although  the  relativistic 
influence  cannot  be  truly  physically  separated  from  the 
charge-displacement  process,  it  has  been  examined  sepa¬ 
rately  because  this  comparison  provides  insight  on  the 
dynamics  of  the  focusing  action,  specifically  the  process 
by  which  the  charge-displacement  and  the  confined  prop¬ 
agation  in  the  electronically  cavitated  channel  develop. 
Since  both  the  relativistic  effect  and  the  charge  displace- 
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ment  tend  to  locally  produce  a  reduction  in  the  plasma 
frequency  {(Op)  that  is  more  significant  in  the  high- 
intensity  regions,  both  of  these  effects  perturb  the  wave 
fronts  in  a  manner  that  encourages  convergence  of  the 
wave  and  the  formation  of  localized  high-intensity  zones. 
Furthermore,  for  both  the  purely  relativistic  and  the 
combined  cases,  we  have  calculated  the  stationary  eigen- 
modes  of  the  nonlinear  Schrodinger  equation  and  show 
their  relationship  to  the  modes  of  propagation  dynami¬ 
cally  developed  in  the  plasmas.  In  particular,  it  is  shown 
that  the  charge  displacement,  especially  that  resulting  in 
dynamical  cavitation,  has  a  very  strong  effect  on  the  spa¬ 
tial  character  of  the  propagation  [26]. 

The  calculations  indicate  that  the  efficiency  of 
confinement  of  the  propagating  energy  is  potentially 
high,  namely,  that  a  large  fraction  of  the  incident  power 
can  be  trapped  in  a  channeled  mode  by  the  combined  ac¬ 
tion  of  the  relativistic  and  charge-displacement  processes. 
Furthermore,  in  the  asymptotic  regime,  it  is  found  that 
the  channel  is  characterized  spatially  by  an  intensity 
profile  and  electron-density  distribution  corresponding  to 
the  lowest  z-independent  eigenstate  of  the  nonlinear 
Schrodinger  equation.  Overall,  the  principal  conclusion 
of  this  work  is  that,  under  appropriate  conditions,  a  new 
dynamical  mode  of  stable  highly  confined  propagation 
naturally  evolves  for  the  propagation  of  sufficiently  short 
(r,  »r»rg )  pulses  of  coherent  radiation  in  plasmas.  By 
new  dynamical  modes,  we  mean  the  self-channeling  of  the 
radiation  through  the  formation  of  stabilized  electroni¬ 
cally  cavitated  paraxial  modes,  which  result  from  the 
combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms.  Interestingly,  for  ultraviolet 
wavelengths  in  the  (200~300)-nm  range,  the  power  densi¬ 
ties  naturally  associated  with  these  environments  can  ap¬ 
proach  thermonuclear  values. 

In  Sec.  II  the  underlying  physical  concepts  are  present¬ 
ed  and  the  governing  nonlinear  Schrodinger  equation  is 
derived.  Section  III  presents  the  analysis  of  the  eigen- 
modes  of  this  equation.  Representative  results  of  numeri¬ 
cal  simulations  of  the  two-dimensional  axisymmetric 
propagation  are  given  in  Sec.  IV.  Section  IV  also  con¬ 
tains  a  comparison  with  the  purely  relativistic  case.  Fi¬ 
nally,  the  conclusions  are  summarized  in  Sec.  V. 

11.  GENERAL  CONSIDERATIONS 
A.  Physical  model 

Several  physical  phenomena  [1-7,14-19,24-26,29] 
play  a  role  in  the  nonlinear  dynamics  gove**ning  the  prop¬ 
agation  of  intense  coherent  radiation  under  the  condi¬ 
tions  being  examined  in  this  study.  They  are  the  follow¬ 
ing. 

(a)  The  creation  of  a  plasma  column  by  ionization  in 
the  temporally  early  region  of  the  laser  pulse. 

(b)  The  influence  of  the  spatial  variation  of  the  refrac¬ 
tive  index  arising  from  the  nonlinear  response  of  the 
dielectric  properties  of  the  medium.  Two  mechanisms 
are  related  to  the  electronic  component;  specifically,  the 
relativistic  shift  in  electron  mass  and  the  ponderomotive- 
ly  driven  electron  motion  which  tends  to  displace  elec;:^ 


trons  from  the  high-intensity  zone.  For  sufficiently  short 
pulses,  only  the  electrons  are  expelled  from  the  laser 
beam  and  the  more  massive  ions,  due  to  their  substantial¬ 
ly  greater  inertia,  are  regarded  as  motionless  [17,18,25]. 
A  third  mechanism  is  the  nonlinear  response  arising  from 
the  induced  dipoles  of  the  ions,  but  this  is  generally  small 
and  negligible  [25]. 

(c)  Defocusing  mechanisms,  caused  by  diffraction  from 
the  finite  aperture  and  refraction  by  the  transverse  inho¬ 
mogeneities  in  the  electron  density. 

(d)  Dissipation  of  laser-beam  energy  by  (i)  motion  of 
the  electrons,  (ii)  ionization  of  the  gas  atoms,  (iii)  genera¬ 
tion  of  harmonic  radiation,  (iv)  production  of  inverse 
bremsstrahlung,  (v)  Compton  scattering,  and  (vi)  other 
amplitudes  of  nonlinear  scattering. 

The  present  work  incorporates  four  phenomena:  (1) 
the  nonlinear  response  of  the  refractive  index  of  the  plas¬ 
ma  due  to  the  relativistic  increase  in  the  mass  of  the  free 
electrons,  (2)  the  refractive  index  variation  due  to  the  per¬ 
turbation  of  the  electron  density  by  the  ponderomotive 
force,  (3)  the  diffraction  caused  by  the  finite  aperture  of 
the  propagating  energy,  and  (4)  the  refraction  generated 
by  the  transversely  inhomogeneous  plasma  density  associ¬ 
ated  with  the  formation  of  a  plasma  column.  In  these 
calculations,  preformed  static  plasma  columns  were  used 
in  order  to  approximate  the  radial  distribution  of  ioniza¬ 
tion  that  is  expected,  if  the  incident  laser  pulses  were  pro¬ 
ducing  the  ionization  on  their  temporally  leading  edge. 
Finally,  the  calculations  were  performed  for  a  length  of 
propagation  that  is  much  shorter  than  the  characteristic 
length  for  dissipation  of  the  energy  of  the  pulse  through 
ionization  or  other  modes  of  energy  loss. 

It  should  be  noted  that  the  relativistic  intensities 
characteristic  of  the  phenomena  examined  in  this  work 
can  be  currently  obtained  experimentally  [30-32].  In 
particular,  the  experimental  parameters  characteristic  of 
the  ranges  that  would  apply  to  the  study  of  these  phe¬ 
nomena  are  presented  in  Table  I.  With  the  wavelength 
and  range  of  electron  densities  shown  in  Table  I,  the  plas¬ 
ma  is  always  underdense,  namely,  icop/oj)^  «\,  where 
cx)  =  2Trc/X  is  the  angular  frequency  of  the  laser  radiation 
and  (i)p—{4Tre^N^/mg)^^^  is  the  customary  plasma  fre¬ 
quency. 

B.  The  propagation  equation 

Consider  the  propagation  of  an  intense  ultrashort  laser 
pulse  in  a  plasma  with  an  initially  radially  inhomogene¬ 
ous  electron  density,  described  by  the  function  fir),  so 


TABLE  I.  Experimental  parameters  characteristic  of  the 
ranges  that  would  apply  for  the  study  of  relativistic  and  charge- 
displacement  self-focusing  of  laser  pulses  in  plasmas. _ 


Peak  intensity 

W/cm^ 

Pulse  length 

T==  100-1000  fs 

Initial  focal  spot  radius 

ro=sl-3  ^m 

Wavelength 

X  2:0. 248  nm  (KrF*  laser  light) 

Target  gas  densities 

p=10“-10^®  cm"^ 

Initial  unperturbed 
electron  density 

AT,  ==10’’- 10^'  cm-^ 
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that  max/(r)=l.  We  denote  the  vector 

and  the  scalar  potentials  of  the  electromagnetic  fields  as 
A  and  respectively,  and  the  corresponding  electric 
field  as  E.  Let  the  momentum  of  the  electrons,  the 
current  density,  and  the  charge  density  be  denoted  as  P*, 
j,  and  p,  respectively.  We  assume,  for  the  short  timescale 
of  interest,  that  the  ions  are  inertially  frozen  in  space 
[17,18,24,25].  Then, 

□  A=c"'V-^-(47r/c)j  ,  (1) 

0/ 

V^^=— 47rp, 

(V- A)=0  ,  (3) 

+c"'[VeX(VX  A)]  ,  (4) 

i=-eN,y,,  p=e(N'°'-N,)  ,  (5) 

v,=p,/m,, 

r  =  tl  +  IPei'/(m,,oC)']'"'  . 

In  the  set  of  statements  above,  Eqs.  (1)  and  (2),  are  the 
Maxwell  equations,  Eq.  (3)  is  the  Coulomb  gauge  condi¬ 
tion,  Eq.  (4)  is  the  equation  of  motion  of  the  electrons, 
Eqs.  (5)  are  the  definition  of  the  current  density  and  the 
charge  density  is  the  initial  charge  density,  while  N, 
is  the  dynamical  charge  density),  and  Eqs.  (6) -represent 
the  relativistic  relation  between  the  velocity  v,  and  the 
momentum  of  the  electrons.  In  Eqs.  (6),  m,  o  is  the  rest 

electron  mass,  □  =  — c  “  ^3^  /3/ 

It  is  convenient  to  normalize  the  values  in  the  equa¬ 
tions  presented  above  as  follows: 

A  =  (e/mjoc’)  A,  $={e /m^QC^)<l>  , 

E=(e/m,.oC^)E  ,  =  .  (7) 

V,  =yjc,  N^=  N, /JV,_o  . 


1  3fl 

3a 

"  3r  ■ 

3z 

both  space  and  time  is  much  greater  than  both  the  plas¬ 
ma  (2irc/(0pQ)  and  electromagnetic  (A.)  wavelengths. 
This  leads  to  the  validity  of  the  inequalities 


«|fca|,|fcpa| 


where  k^=k\—kp,  kQ—co/c,  and  k^—a^Q/c.  We  use 
the  notation  a)^.o=4ire^7V,,o/'«*.o  specificaUy  to  denote 
the  unperturbed'  plasma  frequency.  With  the  assump¬ 
tions  and  approximations  stated  above,  Eqs.  (1)— (8)  be¬ 
come 

□  A  =  fc/A^er"' A  ,  (iO) 

VH  =  kl[N,-f(r)]  ,  (11) 

V((f>-y)=0  ,  (12) 

r=(i+lpj')‘"',  P*=A.  (13) 

The  term  c“‘3V<^/3r  is  omitted  in  Eq.  (10),  since  y  and 
<l>,  according  to  Eqs.  (9),  (12),  and  (13),  do  not  have  a 
high-frequency  dependence.  Equations  (11)  and  (12)  re¬ 
sult  in  the  expression  for  the  electron  density, 

Ar,=max[0,/(r)+V^VV]  .  (14) 

The  logical  function  max(0,  )  provides  for  the  physical¬ 
ly  obvious  and  necessary  condition  Nf>0.  The  analo¬ 
gous  expression  for  the  electron  density  has  been  previ¬ 
ously  derived  by  Sun  et  al.  [18]  for  the  case  f(r)=\.  It 
should  be  noted  that  Eq.  (12)  states  the  condition  for  the 
balance  of  the  ponderomotive  and  the  electrostatic  forces 
for  the  relativistic  case.  Through  combining  Eqs.  (9), 
(10),  and  (14),  we  establish  the  equation  for  the  slowly 

varying  complex  amplitude  of  the  vector  potential 

a(r,z,t)  as 

'  1  3^3 

Vg  dt  3z 

+  -j(Vja+k^{l-y~'maxlO,/(r}  +  k-^Vly]}a) 


with  the  understanding  that,  henceforth,  the  tilde  sign 
will  be  suppressed.  Using  the  relations  (p, ‘Vlp, 
=  VipjV2-p,X(VXp,)  and  Vy  =  V|pjV2y,  Eq.  (4) 
becomes 

|-(p,- A)-v,X[VX(p,- A)]  =  V(<^-r)  •  (8) 

In  the  limit  t»2v/cj^,  the  expression  p,  =  A  is  approxi¬ 
mately  valid.  As  we  shall  see  below,  this  condition  means 
that  the  electron  response  can  be  regarded  as  adiabatic. 
Furthermore,  assuming  the  vector  potential  to  be  circu- 
lariy  polarized,  we  write 

A(r,z,t)  =  ||(e,+iey)a(r,r,/)exp[/(6>t-A:z)]+c.c.  ]  . 

(9) 

Consistent  with  the  statements  made  above,  we  em¬ 
phasize  the  use  of  the  assumption  that  the  pulse  length  in 


In  accord  with  the  previous  assumptions  used  in  deriving 
Eq.  (15),  we  have  neglected  the  second  z  and  t  derivatives 
in  this  expression.  In  Eq.  (15),  Vg  =€6^^  is  the  group  ve¬ 
locity  of  light  in  unperturbed  plasma,  €q=1—Io)j,o/ci))^ 
is  its  corresponding  dielectric  constant,  and 
y =(  1-1-  |a|2)i/2.  It  should  also  be  noted  that  the  electric 
field  vector  E  and  the  vector  potential  A  are  related 
through  the  approximate  relation 


E(r,z,t)=:-/A:o  A(r,z,f)  . 

The  calculation  of  the  propagating  wave  form  for  the 
central  high-intensity  region  of  the  pulse  is  accomplished 
by  considering  the  solutions  of  Eqs.  (15)  along  its  charac¬ 
teristics.  Changing  the  variable  t  to  q  =  t  —z/Vg,  Eq.  (15) 
4  becomes 
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■^o  +  -^(Vja+fcp^{l-r  'max[0,/(r)+;:p 

=0  ,  (16) 

for  which  we  seek  solutions  of  the  form  q  =const.  We  re¬ 
call  that  certain  solutions  of  Eq,  (16)  for  the  special  case 
/(r)=  1  have  already  been  presented  [18,19,33]. 

The  results  of  the  computations  presented  in  this  work 
are  given  with  respect  to  the  variables  (q,r,z). 
Specifically,  the  data  presented  in  the  figures  below  illus¬ 
trate  the  propagation  of  the  radiation  along  the  z  axis  for 
const.  Equation  (16)  describes  the  two-dimensional 
(r,z)  dynamics  of  propagation  for  coherent  circularly  po¬ 
larized  radiation  in  plasmas  having  an  inhomogeneous 
transverse  (r-dependent)  distribution  of  electron  density. 
The  basic  physical  phenomena  explicitly  embodied  in 
these  equations  have  been  outlined  in  Sec.  II  A. 

III.  STATIONARY  SELF-LOCALIZED  MODES 
OF  PROPAGATION 

Equation  (16),  in  the  case  of  homogeneous  plasmas 
[/( r)  =  1  ],  has  axisymmetric  partial  solutions  of  the  form 

a{r,z)=U^^{kpr)t\i[)[i{kj/2k){s—\)z]  ,  (17) 

where  s  is  a  real-valued  dimensionless  parameter  and  the 
real-valued  function  ^  obeys  the  ordinary  differential 
equation  [18] 

(18) 

In  this  discussion  the  dimensionless  argument  p  =  kpr  is 
used  and  V|  in  Eq.  (18)  denotes  the  Laplacian  with  the 
derivatives  designated  with  respect  to  this  variable  (p). 
The  nonlinear  term  in  Eq.  (18)  is 

„  =niax(0, 1 +Vjyj  „  )  ,  (19) 

The  natural  boundary  conditions  for  Eq.  (18)  are 

dU,„{0) 

- - =0,  C/,^(oo)-0.  (20) 

dp 

The  first  condition  assures  axial  symmetry  of  the  solu¬ 
tions,  while  the  second  is  necessary  for  the  property  of 
finite  energy  for  the  solutions.  We  designate  these  solu¬ 
tions  as  “eigenmodes”  and  their  significance  is  explored 
below. 

Consider  initially,  however,  the  eigenmodes  corre¬ 
sponding  to  the  purely  relativistic  case,  functions  which 
are  obtained  from  Eqs.  (18)  and  (19)  by  neglecting  the 
term  Since,  we  are  explicitly  eliminating  the  pon- 

deromotive  potential  and  its  influence  on  the  motion  of 
the  electrons,  in  this  situation  the  electron  density  will 
not  be  self-consistent.  In  this  restricted  case,  Eq.  (18)  can 
be  rewritten  in  the  form 


2 

+  V{U,„,s) 

=-p-' 

dp 

2 

5,  n  ^ 

‘^P 

where  K(C/,,„,5)=(s/2)C/^„-(l  +  C/,^„)‘''^  +  l. 

The  second  boundary  condition,  stated  in  Eq.  (20),  en¬ 
sures  finiteness  of  energy  only  in  the  situation  where  the 
dynamical  system  defined  by  the  last  two  expressions  has 
three  rest  points,  one  of  them  being  zero,  that  is  for  the 
case  0<5  <1.  Exactly  as  in  the  classic  theory  of  cubic 
media  [34,35]  for  values  of  the  parameter  s  belonging  to 
this  interval,  the  localized  eigenmodes  constitute  a  count¬ 
able  set  that  is  ordered  by  the  number  of  zeros -n  in  each 
of  its  members  for  finite  values  of  p.  Figure  1  illustrates 
the  first  four  eigenmodes  for  s  =0. 95. 

Consider  now  the  eigenmode  corresponding  to  the 
unrestricted  situation  that  includes  both  the  relativistic 
and  charge-displacement  mechanisms,  a  case  which  is 
found  to  be  fully  analogous.  The  eigenmodes  of  the  equa¬ 
tion  describing  the  combined  relativistic  and  charge- 
displacement  processes  for  the  interval  0<5<1  also 
make  a  countable  set  that  can  be  ordered  in  the  same  way 
as  the  one  associated  with  the  purely  relativistic  case. 
The  zeroth  (lowest)  eigenmode  U^Q(p)  with  s=0.95, 
along  with  the  corresponding  electron  density  N^q(p), 
given  by  the  second  formula  in  Eq.  (19),  are  depicted  in 
Fig.  2.  These  two  functions  have  been  developed  in  ear¬ 
lier  work  [18].  The  first  and  second  relativistic  and 
charge-displacement  eigenmodes  for  the  value  of  the  pa¬ 
rameter  5=0.95  and  the  corresponding  electron-density 
eigenmodes  ^5,i(p)  and  N^J<p)  are  depicted  in  Figs.  3  and 
4,  respectively.  These  higher  eigenmodes  exhibit  the  im¬ 
portant  feature  that  cavitation  may  occur  in  the  electron¬ 
ic  component  of  plasma  [see  Fig.  (4)],  even  if  it  does  not 
occur  in  the  case  of  the  lower  eigenmodes  [see  Figs.  2  and 
3], 

A  prior  study  [18]  treated  the  (r,z)  dynamics  of  the 
evolution  of  the  perturbed  lowest  eigenmodes  for  the  case 
in  which  no  cavitation  of  the  electronic  component  of 
plasma  occurs.  The  lowest  eigenmodes  were  perturbed 
by  multiplying  them  by  constants  in  the  neighborhood  of 
unity.  In  this  case,  the  known  regime  involving  oscillato¬ 
ry  propagation  was  observed.  However,  the  instability  of 
the  numerical  method  applied  in  that  work  [18]  prevent¬ 
ed  the  performance  of  similar  computations  for  cases  in¬ 
volving  cavitation.  Our  calculations  show  that  the  analo¬ 
gously  perturbed  lowest  eigenmodes  V^^oip)  also  propa¬ 
gate  in  the  same  oscillatory  fashion  when  electronic  cavi¬ 
tation  occurs,  the  corresponding  density  of  which  is  given 
by  N^Qip),  It  is  also  found  that  the  same  behavior  devel- 


-  ;  0  J 


FIG.  1.  Stationary  axially  symmetric  eigenmodes  with 
5  =0.95  corresponding  to  purely  relativistic  self-focusing. 
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FIG.  2.  The  zeroth  (lowest)  axially  symmetric  stationary 
eigenmode  with  5=0.95  corresponding  to  relativistic  and 
charge-displacement  self-focusing:  the  normalized  field  ampli¬ 
tude  distribution  and  the  normalized  electron  density 

distribution  N^Qip). 


ops  for  the  perturbed  first  eigenmodes,  which  have  one  or 
two  cavitated  channels  in  the  electronic  density  described 
by  N^  i(p).  The  results  of  these  simulations  can  be  re¬ 
garded  as  evidence  of  the  (r,2)  stability  of  axially  sym¬ 
metric  lowest  and  first  eigenmodes  against  small  pertur¬ 
bations  in  the  amplitude.  It  should  be  noted,  however, 
that  these  higher  eigenmodes  are  presumably  unstable 
against  small  azimuthal  perturbations,  the  nature  of 
which  violates  the  assumed  axial  symmetry  of  the  distri¬ 
butions. 

The  important  consequence  of  the  character  of  the 
solutions  discussed  above  is  the  fact  that  the  power 

contained  in  the  intensity  distribution  corresponding  to 
the  lowest  eigenmode  unlike  the  case  of  a  cubic 

medium,  depends  on  the  parameter  s.  Namely,  this 
power  decreases  as  s  increases. 

The  infimum  of  by  s,  in  the  internal  0<5  <  1,  is 
called  the  critical  power  (Pcr^  relativistic  and 

charge-displacement  self-focusing.  This  power  equals  the 
critical  power  of  the  purely  relativistic  self-focusing 
[18,36]  and 

^cr=  inf  ,^i=  ii'n^'P,  • 

0<s<l  i  — I-O 


FIG,  3.  The  first  axially  symmetric  stationary  eigenmode 
with  5=0.95  corresponding  to  relativistic  and  charge- 
displacement  propagation:  the  normalized  amplitude  distribu¬ 
tion  and  the  normalized  electron-density  distribution 
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FIG.  4.  The  second  axially  symmetric  stationary  eigenmode 
with  5=0.95  corresponding  to  relativistic  and  charge- 
displacement  propagation:  the  normalized  amplitude  distribu¬ 
tion  U^iip)  and  the  normalized  electron-density  distribution 
^..2(P). 

The  explicit  evaluation  of  this  critical  value  is  presented 
below.  As  U^Q(p)—*^0,  N^Q(p)-^  \  for  s— ►!— 0;  we  have 
[18]  therefore, 

«  (i  +  ca/o/2)~'-i  *  -uUn. 

j^l-0  ’  5-.I-0 


Furthermore, 

J7,.o(p)  ~  t/o(p), 

— 0 

where  U^ip)  is  the  positive,  monotonically  decreasing 
(reaching  zero  for  no  finite  value  of  p)  solution  to  the 
boundary-value  problem  [18,36] 

ViC/o-et/o+jC^^=0, 

dUo 

-^{0)=0,  C/o(oo)=0, 

dp 

with  e=  1  ’^5.  A  change  of  the  variables,  which  is  a  stan¬ 
dard  procedure  for  treating  the  case  of  the  cubic  non¬ 
linearity  [37],  shows  that  [36] 

l7o(p)=(2e)'^Vf‘^V)  . 

where  go  the  customary  Townes  mode,  i.e.,  the  posi¬ 
tive,  monotonically  vanishing  solution  to  the  following 
boundary-value  problem: 

^(0)-0,  go(oo)  =  0. 

dp 

From  the  definition  of  P<,r  relation  between  U^ip) 

and  go(p)>  it  follows  that 

f’cr=  inf  lim 

0<5<l  S  — 1“0 

=2irJJul{p)pdp 

=4w^^gl(p)pdp  =  lP„,  , 

where  P^r.c  =2^7 ^ Qg\(p)pdp  is  the  critical  power  of  the 
Kerr  self-focusing  in  cubic  media  [34,37]. 

Using  the  resulting  relation  between  the  normalized 
values  of  the  critical  power  of  the  relativistic  and  charge- 
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displacement  self-focusing  and  the  critical  power  of 
the  Kerr  self-focusing  ha\c  the  elementary  result 

that 

P  =2P 

Finally,  by  calculating  g^ip)  numerically,  we  find  the  nor¬ 
malized  value  of  Per  high  precision  [36]  to  be 

/>„  =  2i>„,,=4,r /^"g^(p)pJp^23.4018  . 

Hence,  we  establish  the  expression  for  the  critical  power 
[36]  as 

gl(p)pdp{co/o}pQ? 

=rl.6198XlO‘°(a>/a)p,o)^  W  . 

The  constant  factor  involved  in  this  statement  improves 
on  that  given  previously  [18]. 

The  relation  between  U^ip)  and  g^ip)  enables  the  de¬ 
velopment  of  the  asymptotic  expressions  (in  the  case 
6—0)  for  both  the  peak  value  of  the  amplitude  of  the  sta¬ 
tionary  eigenmode  and  for  the  radius  of  the  eigenmode  by 
using  the  appropriate  characteristics  of  the  Townes 
mode.  The  results  are 

t/,,o(0)^(2e)'^2go<0) 

and 


for  6—0. 

Special  computations  were  performed  in  order  to 
determine  how  close  the  values  for  U^qIO)  and 
(/of0)  =  {26)*^^go^0^  ^re  in  the  case  where  6«1.  The 
calculations  show  that  |  is 

2.026X10"^  and  2.096X  10“^  for  6=  10"^  and  10"^  re¬ 
spectively. 

IV.  THE  TWO-DIMENSIONAL  CASE 

It  is  convenient  to  treat  Eq.  (16)  numerically  using  the 
normalized  coordinates 

r,  =r/ro  ,  2,  =z/{2krl)  , 

)  =  aQ^a{r,z)  , 

where  is  a  characteristic  radius  of  the  initial  intensity 
profile  and  for  selected  value  of  q. 

For  simplicity,  we  put  f \{ri  )=/(r)  and  omit  below  the 
subscript  1  of  Tj  and  2|  for  brevity.  Thus,  the  mathemat¬ 
ical  statements  can  be  expressed  in  the  following  set  of 
equations: 

|^^  +  /Vju+(F(/,(r),|uP)u  =0,  r>0  (21) 

02 

u(r,0)  =  Uo(r),  max|uol  =  l  ^22) 

r 

4^(0, z)  =  0,  «(oo,z)=0.  (23) 

or 

The  nonlinear  term  F  is  the  real-valued  operator  c 


F(/„|)=a,(l-(l+a2^)-'^2 

Xmax[0,/,(r)+af'V^(l+a2S)'^^]i  , 

(24) 

and  the  dimensionless  parameters  a  ^ ,  ^2  are  defined  as 

(25) 

where  Iq  is  the  peak  intensity  at  the  entrance  plane 
(2=0)  of  the  medium.  The  parameter 
/^^m^^QCO^c^/(47re^)  is  known  as  the  relativistic  intensity 
[15]. 

The  ratio  of  the  power  of  the  beam  Pq  and  the  critical 
power  of  the  relativistic  and  charge-displacement  self- 
focusing  Fcr»  defined  in  Sec.  Ill,  is  an  important  parame¬ 
ter  characterizing  Eqs.  (21)-(24).  In  the  notations  of  the 
present  section,  the  value  of  Po/P^r  can  be  expressed  in 
the  following  way  [36]: 

Po/P„  =  {aia2/B)  fJ‘\uo(r)\^rdr  , 
with  dimensionless  constant  B  given  by  [36] 

5=2 fj’gi(p)pdp^3.724  5l  . 

When  the  initial  transverse-intensity  distribution  is 
Gaussian,  namely,  luo(r)P==exp(  —  r^),  we  have 

Po/P,,=a,a2/(2B)  . 

We  note  that  in  several  other  studies,  the  critical  power 
of  the  relativistic  self-focusing  is  defined  alternatively,  ba¬ 
sically  with  B=4.0.  In  the  present  discussion,  the  ex¬ 
pression  for  Pq/P^j.  involves  a  value  of  B  <4.0,  so  that 
the  corresponding  ratio  Pq/P^t  is  raised.  Therefore,  a 
pulse  having  a  Gaussian  initial  transverse-intensity  distri¬ 
bution  and  a  flat  initial  wave  front  with  the  parameters 
fli  =248.6192  and  a2“0*031,  undergoes  self-focusing.  In 
this  case,  the  ratio  Pq/P^^,  as  defined  above,  evaluates  to 
1.0347,  namely,  ^0  >  Per-  However,  using  the  value 
jB=4.0  would  give  Pq<P^^,  a  statement  contradicting 
the  results  of  the  computations. 

A.  The  initial  conditions 

In  this  section,  we  examine  the  self-focusing  of 
coherent  radiation  for  pulses  having  Gaussian  or  hyper- 
Gaussian  transverse  and  longitudinal  intensity  distribu¬ 
tions  [24]  of  the  form 

/|^=o=/o(^^)=^mexp[-(^/r)^‘-~(r/ro)^^]  , 

N^>2,  N2>2  (26) 

with  r  and  t  being  dimensional.  We  assume  [24]  that 

^/^»/*  =  10^^  Wem”^  with  /*  designating  the  ap¬ 
proximate  value  of  the  threshold  for  rapid  nonlinear  ion¬ 
ization  [1,2].  The  spatial  amplitude  distribution  of  the  in¬ 
cident  radiation,  defined  by  Eq.  (26)  for  the  case  of  a  flat 
incident  phase  front,  is  of  the  form 
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Uo(r)=exp(  — r^^/2)  ,  N2-2  .  (27) 

The  dimensionless  parameter  Oj  in  Eq.  (25),  correspond¬ 
ing  to  the  incident  pulse  intensity  on  the  axis  (r=0)  of 
the  beam  /off),  is 

a2=/o(f)//,  =  (/„//r)exp[-(f/r)‘'^']  .  (28) 


The  transverse  profile  of  the  plasma  column,  created  by 
the  temporally  leading  edge  of  the  pulse,  is  simulated  by 
the  hyper-Gaussian  function 

f{r)  =  exp[-(r/r,f^],  N2>2  ,  (29) 

where  r  is  the  radial  coordinate.  The  aperture  of  the 
plasma  column  can  be  estimated  [24]  with  the  use  of 
the  relation 

/o(f..'o)-'^o('o)exp[-('-,/'-o)^^]=-^*  • 

For  example,  in  the  case  of  the  Gaussian  transverse  inten¬ 
sity  distribution  (.N2=2),  the  aperture  of  the  plasma 
column  for  /*  =  10**  Wcm  \  /,.=:0.45X  10^°  Wcm 
and  /o(fo)  =  (0.1)/„  gives  r,  =2.47ro.  In  thi^situation, 
the  homogeneous  plasma  approximation  /(r)=l  is  valid 
[24].  In  contrast,  it  follows  from  Eq.  (30)  for  plateaulike 
incident-transverse-intensity  distributions  that  the  aper¬ 
ture  of  the  simulated  plasma  column  tends  to  the 
value  of  the  beam  aperture  Tq.  Therefore,  in  the  example 
given  above  for  iV2  — 8,  we  have  r,  —  1.25ro.  Thus,  de- 
focusing  of  the  beam,  which  is  significant  because  of  the 
near  coincidence  of  the  apertures  of  the  laser  beam  and 
the  plasma  column,  must  be  taken  into  account  when  the 
evolution  of  a  beam  with  a  plateaulike  incident- 
transverse-intensity  distribution  is  studied  [24]. 


B.  Relativistic  self-focusing 


Consider  the  two-dimensional  (r,2)  solutions  of  the  sys¬ 
tem  of  equations  embodied  in  Eqs.  (21),  (22),  and  (23)  for 
the  purely  relativistic  nonlinear  term,  the  form  of  which 
can  be  obtained  from  Eq.  (24)  by  disregarding  the  term 
involving  Vj.  In  this  case,  we  find 

F{fu^)=a,[\-fy{r){l+a2i)-^'^]. 


The  relativistic  self-focusing  mechanism  prevails  over  the 
charge-displacement  mechanism  outside  of  the  focal  spot 
under  the  conditions  fl]  »1,  <12  — (Discussion  of  the 
conditions  for  the  prevailing  of  the  relativistic  self- 
focusing  can  also  be  found  in  Ref.  [33]).  The  situation 
represented  by  Eqs.  (21),  (22),  and  (23),  with  the  above 
nonlinear  term,  describes  self-focusing  with  a  nondissipa- 
tive  saturation  of  the  nonlinearity.  The  properties  of  the 
solutions  for  this  case  depend  essentially  on  the  values  of 
two  conserved  integrals  given  specifically  by 

F,  =  r  '|u(r)|Vdr  (31) 

•>0 


and 


in 

dr 


— <^(r,  |uP) 


r  dr  , 


(32) 


where 

<^(r,i)=  f^F{r,T))dv 

=fl,f|-(2/fl2)/i(r)[(l+a2|)'^^-l]l  •  (33) 

Figure  5  illustrates  the  nature  of  the  calculated  solu¬ 
tions  corresponding  to  this  purely  relativistic  case.  The 
parameters  of  the  incident  radiation  and  plasma,  for  the 
examples  presented  in  Fig.  5,  are  A. =0.248  /xm, 
/,=0. 45X10^°  Wcm'^  /o  =  j/r  =  3X  10*^  Wcm"% 
ro  =  3  /xm,  and  iV.^o^'^.SX  10^°  cm"^  The  correspond¬ 
ing  values  of  the  associated  dimensionless  parameters  are 

a ,  2. 486  192 X 10^  and  G2  =  f 

Figure  5(a)  presents  the  result  for  purely  relativistic 
propagation  in  a  homogeneous  plasma  along  the  z  axis 
for  a  pulse  having  an  incident  Gaussian  transverse  inten¬ 
sity  distribution  and  a  flat  wave  front  [^2  =  2  in  Eq.  (27)] 
for  the  value  of  q  defined  in  Sec.  II B,  corresponding  to 
/o  =  }/^.  Figure  5(b)  presents  the  analogous  graph  for  a 
plateau  like  incident  transverse  intensity  distribution 
[iV2==8  in  Eq.  (27)].  We  conclude  from  these  results  that 
the  solution  is  critically  dependent  upon  the  initial  condi¬ 
tion  represented  by  the  detailed  character  of  the  incident 
transverse  intensity  distribution. 

The  ratio  of  the  beam  power  (Pq)  to  the  critical  power 
(Pgf)  of  relativistic  self-focusing  for  the  given  values 
of  the  parameters  Gj  and  Gj  yields  — 22,252  in  the 

case  of  the  Gaussian  initial  transverse  intensity  distribu¬ 
tion  [iV2  =  2  in  Eq.  (27),  Fig.  5(a)]  and  Po/^cr  =  20. 168  in 
the  latter  example  of  the  plateaulike  [iV2  =  8  in  Eq.  (27), 
Fig.  5(b)]  initial  transverse  intensity  distribution. 

Note  that  for  the  values  of  the  parameters  Ui,a2 
incident  wave  forms  studied,  we  have  from  Eq.  (32) 
Pj  <0.  In  this  situation,  the  following  inequality  is  valid: 

maxlu  (r,2)|^>  (4/G1G2  )|F2l 


namely,  with  respect  to  the  radial  coordinate  r,  the  max¬ 
imum  beam  intensity  has  a  positive  lower  bound  indepen¬ 
dent  of  2.  This  conclusion  can  be  established  in  the  same 
way  that  Zakharov,  Sobolev,  and  Synakh  [37]  demon¬ 
strated  the  analogous  result  in  their  earlier  work  on  self- 
focusing.  Therefore,  a  powerful  relativistic  beam,  which 
self-focuses  in  a  homogeneous  nonabsorbing  plasma 
(/x  =  0.  Ref.  [24]),  results  in  a  field  distribution  represent¬ 
ing  a  pulsing  waveguide  when  Pi  <  0.  Figures  5(a)  and 
5(b)  explicitly  illustrate  the  formation  of  such  a  regime. 
These  pulsing  waveguides  consist  of  alternations  of  ring 
structures  and  focal  spots  on  the  axis  of  the  beam.  The 
power  confined  in  these  complex  modes  represents  ap¬ 
proximately  50%  and  90%  of  the  total  incident  power  for 
the  cases  depicted  in  Figs.  5(a)  and  5(b),  respectively.  It 
should  be  noted  that  oscillating  periodic  solutions  of  this 
type  have  also  been  presented  by  other  workers  [17]. 

In  the  case  Pi  >  0,  an  inequality  comparable  to  Eq.  (34) 
cannot  be  obtained.  Following  the  method  of  Zakharov, 
Sobolev,  and  Synakh  [37],  we  simulate  the  case  Pi>0hy 
considering  a  beam  initially  focused  or  defocused  by  a 
lens  at  the  entrance  plane  (2  =  0)  of  the  medium.  Let  the 
5^ocusing  (or  defocusing)  length  of  this  lens  be  R  —krlRf 
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FIG.  5.  Purely  relativistic  propagation  with  /q  =  3X10‘'’ 
W/cm^,  To  =  3  /im,  A.=0,248  fim,  and  N^q  —  I^SX  10^°  cni”\  (a) 
The  formation  of  the  pulsing  waveguide  regime  in  the  case  of 
the  relativistic  self-focusing  of  a  pulse  with  a  flat  incident  wave 
front;  Gaussian  initial  transverse-intensity  distribution  [A'’2  =  2 
in  Eq.  (27)],  homogeneous  plasma,  (b)  The  formation  of  the 
pulsing  waveguide  regime  in  the  case  of  the  relativistic  self- 
focusing  of  a  pulse  with  a  flat  incident  wave  front;  hyper- 
Gaussian  initial  transverse-intensity  distribution  [^2  =  8  in  Eq, 
(27)],  homogeneous  plasma,  (c)  The  single-focus  regime  in  the 
case  of  the  relativistic  self-focusing  of  an  initially  focused  pulse 
q/2  in  Eq.  (35)];  Gaussian  incident-iransverse-intensity 
distribution  [Ni^l  in  Eq.  (35)],  homogeneous  plasma,  (d)  The 
formation  of  the  quasistabilized  regime  in  the  case  of  the  rela¬ 
tivistic  self-focusing  of  a  pulse  with  a  flat  incident  wave  front  in 
a  plasma  column  [^3=8,  r*=ro  in  Eq.  (29)];  hyper-Gaussian 
initial  transverse-intense  distribution  [V2  =8  in  Eq.  (27)]. 


{Rf  dimensionless).  The  condition  R^>0  signifies  that 
the  pulse  is  initially  focused  and  Rj-<Q  indicates  that  it  is 
initially  defocused.  Then,  the  corresponding  initial  con¬ 
dition  can  be  written  as 

«o(r)=exp[-/V2  +  /rV(2i?^)],  N^>2  .  (35) 

Furthermore,  let  P2 \Rf\=R j-q.  We  note  that  the 
case  P2>0  corresponds  physically  to  a  high  degree  of  the 
initial  focusing  or  defocusing:  Rj^q>  \Rj^\>0,  Moreover, 
when  Rj-q<\Rj-\S  <xiy  P2<0. 

Figures  5(a)  and  5(c)  display  the  calculated  intensity 
distributions  for  the  case  of  relativistic  self-focusing  in 
homogeneous  plasma  of  beams  with  Gaussian  incident- 
transverse-intensity  distributions  and  values  of  Ry  given 
by  Ry=oo  (flat  wave  front)  and  Rf^^/ly  respectively.  A 
more  detailed  study  of  the  dynamics  of  the  transition 
from  the  pulsing  waveguide  regime  [see  Fig.  5(a)],  to  the 
single-focus  regime  [see  Fig.  5(c)  representing  the  propa¬ 
gation  of  an  initially  sufficiently  sharply  focused  beam)] 
can  be  found  in  Ref.  [38].  Importantly,  the  computations 
show  that  the  value  P 2=0  is  not  the  threshold  separating 
these  two  regimes  of  relativistic  self-focusing.  This  tran¬ 
sition  occurs  as  the  first  focus  gains  power  and  is  shifted 
closer  to  the  entrance  of  the  medium  (z=0),  while  the 
remaining  foci  are  shifted  in  the  opposite  direction  and, 
in  the  limit  of  large  displacement,  become  diffused.  Note 
that  single-focus  regimes  of  propagation  have  also  been 
observed  by  Sprangle,  Tang,  and  Esarey  [17]. 

Defocused  beams  can  evolve  in  a  different  fashion.  Ini¬ 
tially  sufficiently  sharply  defocused  pulses  (Ry<0, 
\Rf\^RfQ/2y  P2>0)  propagating  in  the  relativistic  re¬ 
gime  monotonically  diffuse  on  the  radial  periphery  and 
do  not  exhibit  the  phenomenon  of  self-focusing. 

The  detailed  spatial  character  of  the  plasma  column 
can  have  a  strong  influence  on  the  evolution  of  the  propa¬ 
gation.  Figure  5(d)  illustrates  the  relativistic  propagation 
in  a  plasma  column  along  the  z  axis  for  a  pulse  corre¬ 
sponding  to  transverse  profile  of 

the  plasma  column  is  given  by  fir)  as  defined  by  Eq.  (29) 
with  A^3  =  8  and  r,  =ro.  The  comparison  between  Figs. 
5(b)  and  5(d)  demonstrates  that  the  defocusing  of  a 
hyper-Gaussian  beam  in  a  plasma  column,  with  an  aper¬ 
ture  close  to  the  radius  of  the  beam,  fundamentally  alters 
the  spatial  dynamics  of  the  propagating  energy  [24],  The 
defocusing  causes  a  fraction  of  the  beam  to  spread  away 
from  the  column,  but  the  remaining  energy  of  the  beam 
resolves  into  a  state  that  balances  the  relativistic  self- 
focusing,  defocusing,  and  diffraction.  The  power  trapped 
in  this  quasistabilized  state  so  formed  is  approximately 
25%  of  the  incident  power. 

Previous  analytical  estimates  [38]  have  shown  that  the 
relativistic  self-focusing  length  is  minimal  when  the  in¬ 
cident  pulse  intensity  on  the  axis  of  the  beam  satisfies  the 
condition  lQ(t)  =  2f.  Moreover,  specific  computations 
also  show  that  this  inference  remains  valid  for  the  case  of 
the  relativistic  self-focusing  of  pulses  with  flat  initial 
phase  fronts  and  Gaussian  incident-transverse-intensity 
distributions  in  homogeneous  plasmas. 

The  locus  of  the  first  focus  for  the  case  of  the  relativis¬ 
tic  self-focusing  has  been  presented  previously  [24].  In 
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that  study,  the  beam  initially  had  both  a  flat  phase  front 
and  a  Gaussian  transverse-intensity  distribution.  The  pa¬ 
rameters  were  the  same  as  those  applying  to  Fig.  5(a),  and 
the  duration  of  the  pulse  corresponded  to  t=0.5X  10 
s.  The  minimal  z  of  the  locus  is  reached  when  lQ(t)  =  2Ir- 
If  the  maximum  intensity  on  the  axis  of  the  beam  is  such 
that  >  2I„  then  the  trajectory  of  the  first  focus  exe¬ 
cutes  a  path  that  has  three  reversal  points.  Two  of  them, 
corresponding  to  an  identical  value  of  z,  are  due  to  the 
condition  /o(f)  =  2/,.  The  central  one,  corresponding  to 
a  greater  value  of  z,  arises  when  1^(1)= I In  principle, 
this  behavior  makes  it  possible  to  distinguish  the  process 
of  relativistic  self-focusing  from  that  associated  with  the 
Kerr  nonlinearity,  since  the  latter  results  in  a  locus  of  the 
first  focus  having  only  a  single  reversal  point. 

C.  Relativistic  and  charge-displacement  self-channeling 

This  section  is  devoted  to  the  description  of  the  self¬ 
channeling  occurring  when  both  the  relativistic  and 
charge-displacement  mechanisms  are  included  in  the  in¬ 
teraction.  This  situation  is  described  by  Eqs.  (21)- (24) 
and  includes  the  important  nonlinear  term  involving 
appearing  in  Eq.  (24).  The  numerical  results  show  that 
sufficiently  intense  short  duration  (ri»r»r^)  axisym- 
metric  pulses  readily  undergo  self-channeling  in  plasmas 
over  a  rather  broad  range  of  conditions.  Moreover,  it  is 
found  that  a  large  fraction  of  the  total  incident  power  of 
the  beam  can  be  trapped  in  a  stabilized  mode  confined  to 
the  axis  of  propagation. 

A  specific  example  is  informative  in  representing  the 
general  behavior  exhibited  by  the  propagation  in  the  re¬ 
gime  in  which  the  influence  of  relativity  and  charge  dis¬ 
placement  are  both  significant.  These  results  are  exhibit¬ 
ed  in  Fig.  6.  In  the  case  of  propagation  of  an  initially 
Gaussian  transverse  wave  form  incident  in  a  homogene¬ 
ous  plasma,  with  the  parameters  >.  =  0.248  fxm, 
/o  =  4-/,  =  3XlO'^  Wcm’%  ro  =  3/xm,  and  N,q  =  1.5 
X10^°  cm“^  (a,  =248.6192,  a2  =  T»  /^o/^cr  =  22.252), 
the  numerical  computations  show  that,  as  soon  as  the 
first  focus  on  the  axis  of  propagation  is  formed,  electronic 
cavitation  occurs.  Specifically,  this  leads  to  complete  ex- 
pulsion  of  the  electronic  component  of  plasma  from  the 
paraxial  domain  [Fig.  6(b)].  This  process  results  in  a 
quasistabilized  cavitated  channel  in  the  electron  distribu¬ 
tion  which  extends  along  the  entire  axis  of  propagation 
past  the  location  of  the  initial  focus  [26].  We  note  that 
some  of  the  lowest  stationary  solutions  corresponding  to 
the  relativistic  and  charge-displacement  problem,  includ¬ 
ing  cavitation,  were  developed  by  Sun  et  al.  [18].  The 
analysis  shows  that  the  first  focus  involves  about  45%  of 
the  total  incident  power  of  the  propagating  energy.  A 
fraction  of  the  remaining  power  is  dissipated  through 
diffraction  on  the  periphery,  while  another  component  is 
temporarily  involved  in  the  formation  of  a  pulsing  ring- 
shaped  structure  [Fig.  6(a)].  Subsequent  energy  exchange 
between  the  ringlike  feature  and  the  paraxial  zone  is  ob¬ 
served,  and  as  a  consequence  of  this  interaction,  a  certain 
part  of  the  energy  of  the  pulsing  ring  is  diffracted  away 
while  the  remaining  power  joins  the  paraxial  domain. 

In  finer  detail,  the  following  aspects  of  the  evolution  of 


the  radiative  energy  are  also  revealed  by  the  calculations. 
As  described  above,  after  the  formation  of  the  first  focus, 
considerable  power  is  transferred  from  the  paraxial  focal 
zone  to  the  ring-shaped  feature,  which,  at  this  stage  in 
the  evolution  of  the  pulse,  contains  approximately  68% 
of  its  total  initial  power.  This  intense  ring,  which  spreads 
away  from  the  paraxial  domain,  produces  a  correspond¬ 
ing  ring-shaped  cavity  in  the  electronic  distribution  [see 
Fig.  6(b)].  The  refraction  resulting  from  this  strongly 
perturbed  electron-density  profile,  together  with  the  rela¬ 
tivistic  self-focusing  mechanism,  causes  the  wave  to  re¬ 
turn  energy  to  the  core  of  the  beam.  Thus,  the  charge 
displacement  produces  a  potent  additional  self-focusing 
action,  which  leads  to  the  formation  of  a  confined  paraxi¬ 
al  mode  of  high  intensity  stabilized  along  the  axis  of 
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FIG.  6.  The  self-channeling  of  a  pulse  with  a  Gaussian  initial 
transverse-intensity  distribution  [Ni  =  2  in  Eq.  (27)]  and  a  flat 
incident  wave  front  in  initially  homogeneous  plasma  in  the  case 
of  the  relativistic  and  charge-displacement  propagation  with 
/o  =  3X10'’  W/cm^  ro  =  3  /im,  A.=0.248  fim,  and 
(Vo=(V,o=7.5X  10^°  cm~\  (a)  The  distribution  of  the  normal¬ 
ized  intensity,  (b)  The  distribution  of  the  normalized  electron 
density,  (c)  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  [/,(r)//o]'^^  and  the  normalized  elec¬ 
tron  density  Nsir)/No  for  s  =0.554. 
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propagation.  This  phenomenon  of  confined  propagation 
is  designated  as  self-channeling.  It  should  be  noted  that 
the  calculations  of  the  electrostatic  energy  associated 
with  the  charge  displacement,  which  is  given  by  the  ex¬ 
pression  f  ^Elr  dr,  with  the  electrostatic  field 

Ec  defined  by  the  equation  VE^  =  —  47rp  and  designat¬ 
ing  the  radial  extent  of  the  channel,  show  that  this  energy 
can  be  relatively  small.  Specifically,  for  the  conditions 
represented  in  Fig.  6(b)  [z=95.4  ^m],  the  electrostatic 
energy  (JV)  accounts  for  only  0.18%  of  the  total  energy 
of  the  laser  radiation  per  unit  length. 

The  essential  finding  of  these  calculations  is  that  the 
combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms  produces  a  strong  tendency  for 
the  generation  of  spatially  highly  confined  modes  of  prop¬ 
agation  which  are  stabilized  along  the  axis  of  propaga¬ 
tion.  Furthermore,  the  study  of  a  range  of  other  cases  in¬ 
dicates  that  these  modes  are  exceptionally  stable  and  that 
a  considerable  fraction  of  the  incident  power  can  be 
confined  in  them.  The  result  is  the  controlled  generation 
of  a  very  high  peak  intensity  in  these  channeled  modes 
with  values  reaching  ^  10^^  W/cm^  for  the  range  of  con¬ 
ditions  studied. 

Important  characteristics  of  the  asymptotic  behavior 
of  these  confined  modes  have  also  been  established.  It 
has  been  shown  that  the  distribution  of  the  amplitude 
u(r,2),  for  large  values  of  z,  tends  asymptotically  to  the 
lowest  eigenmode  of  the  nonlinear  Schrodinger  equation 
(see  Sec.  III).  Specifically,  for  the  example  discussed 
above,  the  computations  have  demonstrated  that  the 
asymptotic  radial  amplitude  distribution  corresponds  to 
the  lowest  eigenmode  with  5^0.554.  In^this  case,  the 
asymptotic  intensity  distribution  /^(r)— contains 
46%  of  the  total  incident  power. 

The  normalized  asymptotic  field  amplitude 
and  the  corresponding  normalized  asymp¬ 
totic  plasma  electron  density  A^^(r)/iVo  depicted  in 
Fig.  6(c).  It  should  be  noted  that  the  profiles  of  the  inten¬ 
sity  /(r,z)  and  the  electron  density  iV(r,z),  obtained  as 
the  results  of  the  dynamical  calculations  of  the  propaga¬ 
tion  for  z^900  /im,  differ  from  /^(r)  and  N^(r)  for 
s  =0.554  by  much  less  than  1%.  In  addition,  we  observe 
that  the  energy  of  charge  displacement  for  =2.5  /xm  in 
the  asymptotic  state  accounts  for  a  fraction  of 
9.45  XlO'*^  of  the  total  energy  of  the  beam  per  unit 
length.  We  observe  that  this  tendency  for  a  solution  of  a 
nonlinear  Schrodinger  equation  involving  a  saturable 
nonlinearity  to  converge  to  the  lowest  stationary  solution 
was  originally  discovered  by  Zakharov,  Sobolev,  and 
Synakh  [37]. 

For  the  range  of  parameters  studied,  the  calculations 
clearly  show  that  the  charge  displacement  has  a  very 
strong  influence  on  the  character  of  the  propagation  after 
the  first  focus  is  formed.  The  pulsing  intensity  structure, 
consisting  of  alternating  foci  on  the  axis  of  propagation 
and  peripheral  focal  rings,  which  is  the  usual  behavior  for 
the  purely  relativistic  self-focusing  [see  Figs.  5(a)  and 
5(b)],  is  converted  into  a  stabilized  and  uniform  channel. 
Collaterally,  a  stabilized  cavitated  channel  in  the  electron 
density  is  also  formed.  We  remark  that  the  periodic  in¬ 
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tensity  structure  characteristic  of  the  relativistic  self- 
focusing  occurs  for  great  values  of  z  and,  therefore,  may 
be  regarded  as  the  corresponding  asymptotic  solution  of 
the  purely  relativistic  case.  The  charge  displacement 
leads  to  the  formation  of  the  asymptotic  amplitude  distri¬ 
bution  represented  by  the  corresponding  lowest  eigen¬ 
mode  U^  ol^r)  instead  of  this  pulsing  structure. 

Computations  have  been  also  performed  for  the  propa¬ 
gation  of  incident  plateaulike  wave  forms  with  flat  in¬ 
cident  phase  fronts  in  both  homogeneous  plasmas  and 
plasma  columns.  Additional  calculations,  have  also  ex¬ 
amined  the  behavior  of  focused  Gaussian  and  plateaulike 
incident  wave  forms  as  well  as  defocused  Gaussian  in¬ 
cident  wave  forms  in  homogeneous  plasmas  (see  Figs. 
7-11).  It  is  found  that  the  main  features  described 
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FIG.  7,  The  self-channeling  of  a  pulse  with  a  hyper-Gaussian 
incident-transverse-intensity  distribution  [A^2“8  in  Eq.  (27)] 
and  a  flat  initial  wave  front  in  initially  homogeneous  plasma  for 
the  case  of  a  relativistic  and  charge-displacement  propagation 
with  /o  =  3X10‘^  W/cm^,  /’o  =  3  ^m,  A.  =  0.248  /xm,  and 
^o“^f.o^7.5X  10^°  cm~\  (a)  The  distribution  of  the  normal¬ 
ized  intensity,  (b)  The  distribution  of  the  normalized  electron 
density,  (c)  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  [/,(r)//o]‘^^  and  the  normalized  elec¬ 
tron  density  N,(r)/SQ  for  5  =0.515. 
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above,  namely  ( 1 )  self-channeling,  (2)  stabilization  of  the 
mode  of  propagation,  (3)  the  confinement  of  a  substantial 
fraction  of  the  incident  power,  and  (4)  the  formation  of 
paraxial  cavitated  channels  in  the  electron  distribution 
are  common  aspects  of  the  dynamics  over  a  wide  range  of 
conditions. 

The  character  of  these  findings  is  now  illustrated  with 
five  representative  examples.  For  the  propagation  of  a 
beam  with  an  incident  plateaulike  transverse-intensity 
profile,  given  by  Eq.  (27)  with  N2  =  8,  which  has  a  flat  ini¬ 
tial  wave  front  incident  on  an  initially  homogeneous  plas¬ 
ma  (the  values  of  the  parameters  have  been  adopted  from 
above  and  Pq/P^^  =  1Q.  168),  the  asymptotic  radial  profile 
of  the  intensity  distribution  contains  approximately  77% 
of  the  incident  power  [Fig.  7(a)].  The  corresponding 
asymptotic  amplitude  is  the  lowest  eigenmode  of  the  non¬ 
linear  Schrodinger  equation  with  s  ==0.515,  The  propa¬ 
gation  of  the  same  wave  form  in  a  plasma  column,  with 
the  initial  electron  distribution  defined  by  Eq.  (29)  with 
7^3  =  8  and  results  in  a  quasistabilized  intensity 

distribution  containing  34%  of  the  incident  power  [Fig. 
8(a)].  The  purely  relativistic  propagation  of  incident  pla¬ 
teaulike  pulses  [A^2“^  Eq.  (27)]  in  plasma  columns 
[V3  =8  and  =rQ  in  Eq.  (29)]  also  results  in  the  forma¬ 
tion  of  quasistabilized  regimes  which,  in  this  case,  arise 
dynamically  from  the  defocusing  action  of  the  refraction 


FIG.  8.  The  self-channeling  of  a  pulse  with  a  hyper-Gaussian 
incident-transverse-intensiiy  distribution  [^2=8  in  Eq,  (27)] 
and  a  fiat  initial  wave  front  in  a  preformed  column-shaped  plas- 
ma  [^3=8,  r*  in  Eq.  (29)]  for  the  case  of  relativistic  and 
charge-displacement  propagation  with  /o  =  3X10*'’  W/cm,’ 
To  =  3  /im,  >.=0.248  /4m,  and  Vq  =iV,,o  =  7.5 X  10^®  cm‘^  (a) 
The  distribution  of  the  normalized  intensity,  (b)  The  distribu¬ 
tion  of  the  normalized  electron  density. 


generated  by  the  transverse  profiles  of  the  electron  densi¬ 
ty  [24].  For  this  situation,  approximately  25%  of  the  to¬ 
tal  incident  power  is  confined  [Fig.  5(d)].  The  compar¬ 
ison  of  these  cases  with  the  two  examples  discussed  above 
involving  the  charge-displacement  mechanism  indicates 
that  the  increase  in  the  value  of  the  confined  power  stems 
principally  from  the  substantial  additional  focusing  ac¬ 
tion  arising  from  the  inhomogeneous  electron  distribu¬ 
tion  produced  by  the  ponderomotive  force. 

Overall,  the  computations  reveal  that  the  charge  dis¬ 
placement,  which  generally  results  in  electronic  cavita¬ 
tion,  plays  an  important  role  in  stabilizing  the  mode  of 
confined  propagation  that  develops  dynamically.  This 
stabilization  naturally  occurs  by  refraction  of  the  radia- 
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FIG.  9.  Self-channeling  of  an  initially  focused 
in  Eq.  (35)]  pulse  with  a  Gaussian  incident-transverse-intensiiy 
distribution  [V2  =2  in  Eq.  (35)]  in  initially  homogeneous  plasma 
for  the  case  of  relativistic  and  charge-displacement  propagation 
with  /o  =  3X10‘’  W/cm%  ro  =  3  /im,  X=0.248  /xm,  and 
^o“^ro~7.5X  10^°  cm“\  (a)  The  distribution  of  the  normal¬ 
ized  intensity,  (b)  The  distribution  of  the  normalized  electron 
density,  (c)  Radial  dependence  of  the  asymptotic  solutions  for 
the  normalized  amplitude  [/,(r)//o]‘^  and  the  normalized  elec¬ 
tron  density  V,(r)/Vo  for  s  =0.566. 


5842 


A.  B,  BORISOV  et  al 


45 


tion  into  the  central  paraxial  region.  Indeed,  the 
influence  of  the  charge  displacement  on  the  propagation 
is  so  strong  that  self-channeling  occurs  even  in  the  cases 
of  extremely  focused  and  defocused  incident  wave  forms. 
The  results  of  the  corresponding  calculations  for  ex¬ 
tremely  focused  incident  wave  forms,  for  Ry=Ryo/2 
[Eq.  (35)],  with  both  Gaussian  in  Eq.  (35)]  and 

plateaulike  [A'2  =  8  in  Eq.  (35)]  incident-transverse- 
intensity  distributions,  are  depicted  in  Figs.  9  and  10,  re¬ 
spectively  (the  parameters  of  beam  and  plasma  being  the 
same  as  the  other  examples  above).  In  these  two  cases, 
the  asymptotic  transverse  profiles  of  the  amplitude  and 
the  electron  density  are  found  to  correspond  to  the  lowest 
eigenmodes  of  the  nonlinear  Schrodinger  equation  with 
values  of  the  parameter  5^0.566  (Fig.  9)  and  5=^0.505 


(Fig.  10).  Note,  in  strong  contrast  to  the  situation  involv¬ 
ing  charge  displacement,  the  analogously  strongly  fo¬ 
cused  incident  wave  forms  for  the  purely  relativistic  case 
propagate  in  the  single-focus  regime  [Fig.  5(c)],  and  stable 
confinement  does  not  develop. 

Figure  1 1  illustrates  the  relativistic  and  charge- 
displacement  propagation  of  an  initially  strongly  de¬ 
focused  wave  form  =  — Ry  o/2  in  Eq.  (35)]  having  a 
Gaussian  initial  transverse-intensity  distribution  [N2=2 
in  Eq.  (35)].  After  the  initial  stage  of  the  defocusing,  this 
pulse  evolves  into  a  paraxial  structure  that  is  analogous 
to  those  described  above.  Furthermore,  in  this  example 
the  asymptotic  transverse  profiles  of  the  amplitude  and 
the  electron  density  are  found  to  correspond  to  the  lowest 
eigenmode  of  the  nonlinear  Schrodinger  equation  with 


FIG.  10.  The  self-channeling  of  an  initially  focused 
[Rf  =  Rf  o/2  in  Eq.  (35)]  pulse  with  a  hyper-Gaussian  incident- 
transverse-intensity  distribution  [A,  =8  in  Eq.  (35)]  in  initially 
homogeneous  plasma  for  the  case  of  relativistic  and  charge- 
displacement  propagation  with  /o  =  3X10‘‘^  W/cm%  ro  =  3  ^m, 
A.=0.248  fim,  and  Nq  =A<..o  =  ‘7.5X  10^°  cm~^  (a)  The  distribu¬ 
tion  of  the  normalized  intensity,  (b)  The  distribution  of  the  nor¬ 
malized  electron  density,  (c)  Radial  dependence  of  the  asymp¬ 
totic  solutions  for  the  normalized  amplitude  [/,(r)//o]‘^'  and 
the  normalized  electron  density  :V,(r)AVo  for  5  =0.505. 


FIG.  11.  The  self-channeling  of  an  initially  defocused 
[Rf—  —  Rf^o/2  in  Eq.  (35)]  pulse  with  a  Gaussian  incident- 
transverse-intensity  distribution  [^2=2  in  Eq.  (35)]  in  initially 
homogeneous  plasma  for  the  case  of  relativistic  and  charge- 
displacement  propagation  with  /q  =  3  X  10*’  W/cm^  Tq  =  3  /xm, 
A. =0.248  /xm,  and  Ao  =  A,o=7.5X  10^°  cm’’^  (a)  The  distribu¬ 
tion  of  the  normalized  intensity,  (b)  The  distribution  of  the  nor¬ 
malized  electron  density,  (c)  Radial  dependence  of  the  asymp¬ 
totic  solutions  for  the  normalized  amplitude  [/,(r)//o]*^^  and 
the  normalized  electron  density  A,(r)/Ao  for  s  =0. 800. 


45 


relativistic  and  charge-displacement  self-  . . . 


5843 


s^O.  800. 

For  these  values  of  the  dimensionless  parameters  as¬ 
sumed  in  these  examples  (Figs.  9-11)  (aj  =248.6192, 
that  is  Po/^cr“22.252  for  Gaussian  incident 
wave  form),  the  magnitude  of  the  quantity  Rf=Rf^o/2j 
the  factor  determining  the  degree  of  the  initial  focusing 
(defocusing),  is  0.124.  For  ^=0.248  /xm  and  ro  =  3  /xm, 
the  focusing  distance  R  ^kr^Rf  ^/l,  used  in  the  exam¬ 
ples  above,  has  a  magnitude  of  28.35  /xm.  The  fact  that 
the  self-channeling  of  such  extremely  focused  or  de- 
focused  distributions  occurs  (Figs.  9-11)  means  the  self¬ 
channeling  should  be  expected  to  be  nearly  independent 
of  the  initial  conditions  of  focusing  or  defocusing. 

A  critical  power  for  self-channeling  arising  from 
the  relativistic  charge-displacement  mechanism  can  be 
defined.  We  now  discuss  this  issue  with  a  model 
developed  in  an  earlier  study.  The  threshold  power  is 
defined  as  the  power  that  separates  the  asymptotic  behav¬ 
ior,  with  respect  to  large  distance  of  propagation  (z),  into 
two  distinct  classes.  Since  its  definition  rests  on  an 
asymptotic  property,  it  can  only  exist  in  the  limiting  case 
of  vanishing  losses  (/x=0  in  Ref.  [24]).  For  a  power 
P  <Pcr>  asymptotic  transverse-intensity  profile  tends 
to  zero  at  large  z.  In  contrast,  for  P  >  the  asymptotic 
profile  of  the  intensity  tends  to  the  lowest  eigenmode  of 
the  governing  nonlinear  Schrodinger  equation.  Addition- 

_ _ _ _ i 


ally,  the  normalized  value  of  this  critical  power  is  exactly 
twice  the  corresponding  normalized  value  of  the  critical 
power  for  the  Kerr  effect  self-focusing  found  in  cubic 
media  for  initiaUy  flat  wave  forms  [34,37].  The  expres¬ 
sion  corresponding  to  this  critical  power  [36]  for  the  rela¬ 
tivistic  and  charge-displacement  mechanism  is 

Per  =  ( Ve  ^ )  (P  )p  ^  /^p.o 

^1.6198X10^°((y/6>^,o)^  W  ,  (36) 

where  is  the  Townes  mode  [34]. 

It  has  also  been  shown  that  the  value  of  the  relativistic 
and  charge-displacement  self-channeling  threshold 
power,  in  cases  involving  initially  focused  or  defocused 
beams,  exceeds  P^^  and  depends  on  the  degree  of  the  ini¬ 
tial  focusing  or  defocusing.  For  a  given  magnitude  of  the 
curvature  of  the  wave  front,  the  value  of  the  threshold 
power  for  initially  defocused  beams  is  greater  than  that 
for  initially  focused  beams.  Finally,  it  has  been  shown 
both  analytically  and  numerically  that  for  self-channeling 
of  an  arbitrary  wave  form  to  occur,  it  is  sufficient  that  the 
Hamiltonian  of  the  purely  relativistic  case,  considered  as 
a  functional  of  the  initial  transverse  amplitude  distribu¬ 
tion,  should  be  negative  [36].  The  precise  statement  for 
this  condition  is 


r  dr  <0  . 


(37) 


We  now  comment  on  the  behavior  of  pulses  having  ini¬ 
tial  amplitude  distributions  U^ir)  close  to  higher  eigen- 
modes  Fj,„(r),  n>\.  These  higher  modes  generally  are 
associated  with  electronically  cavitated  channels.  The 
self-channeling  in  these  cases  could  result  in  asymptotic 
distributions  of  intensity  and  electron  density  corre¬ 
sponding  to  certain  higher  eigenfunctions 

particular,  this  conjec¬ 
ture  has  established  an  initial  amplitude  distribution 
f/o('')»  which  is  close  to  the  first  eigenmode  V^  ^ir)  with 
s  =0.544.  Direct  numerical  calculations  showed  that  the 
asymptotic  distributions  of  both  intensity  and  electron 
density  that  evolved  in  this  case  correspond  to  the  first 
eigenmode  N.jir)  with  s  essentially  equal 

to  0.544.  It  should  be  additionally  noted,  however,  that 
these  higher  eigenmodes  ( n  >  1 )  are  quite  possibly  unsta¬ 
ble  against  small  azimuthal  perturbations  that  destroy 
their  axial  symmetry. 

The  principal  result  of  this  section  is  the  finding  that 
the  combined  action  of  the  relativistic  and  charge- 
displacement  mechanisms  can  result  in  self-channeling 
with  the  formation  of  stabilized  paraxial  modes  over  a 
rather  wide  range  of  physical  conditions.  Moreover, 
these  spatially  confined  modes  are  generally  associated 
with  corresponding  cavitated  channels  in  the  electron 
density.  Finally,  the  characteristics  of  these  channeled 


I  ^ 

modes  have  asymptotic  behavior  that  is  described  by  the 
appropriate  lowest  eigenmodes  of  the  governing  non¬ 
linear  Schrodinger  equation. 

V.  CONCXUSIONS 

A  theoretical  approach  suitable  for  the  numerical  in¬ 
vestigation  of  the  two-dimensional  (r,z)  dynamics  of 
propagation  of  coherent  ultrashort  (t,  »t»t^)  relativ¬ 
istic  laser  pulses  in  cold  underdense  plasmas  has  been 
developed.  Four  basic  physical  phenomena  are  included 
within  the  scope  of  this  method.  They  are  (1)  the  non¬ 
linear  dependence  of  the  index  of  refraction  due  to  the 
relativistic  increase  in  the  mass  of  the  electrons,  (2)  the 
variation  of  the  index  of  refraction  resulting  from  the 
perturbation  of  the  electron  density  by  the  ponderomo- 
tive  force,  (3)  the  diffraction  of  the  radiation,  and  (4)  the 
refraction  caused  by  nascent  transversely  inhomogeneous 
electron-density  distributions.  The  equations  studied  in 
this  work  may  be  regarded  as  the  generalization  of  those 
treated  previously  in  other  studies  [18,19,33]  involving  in¬ 
itially  inhomogeneous  plasmas.  Further  studies  are  con¬ 
tinuing  with  an  extension  of  this  analysis,  which  includes 
the  azimuthal  coordinate  in  the  description  of  the  propa- 

i-jgation. 
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The  main  conclusions  stemming  from  the  calculations 
are  the  following. 

(1)  The  cooperative  effect  of  the  relativistic  and 
charge-displacement  mechanisms  leads  asymptotically  to 
stable  high-intensity  z-independent  modes  of  self¬ 
channeling,  and  a  major  fraction  of  the  incident  power 
can  be  confined  in  these  paraxial  modes.  Stable  cavita¬ 
tion  of  the  electron  density  is  a  general  feature  of  these 
spatially  confined  modes. 

(2)  The  z-independent  modes,  to  which  the  solutions  of 
the  equation  describing  the  relativistic  and  charge- 
displacement  propagation  tend  asymptotically,  are  recog¬ 
nized  as  the  lowest  eigenmodes  of  the  governing  nonlinear 
Schrodinger  equation. 

(3)  A  separate  study  of  purely  relativistic  propagation 
shows  that  beams  with  flat  incident  phase  fronts  exhibit  a 
pulsing  behavior  in  homogeneous  plasmas  but  can  under¬ 
go  quasistabilization  in  suitably  configured  plasma 
columns.  However,  sufficiently  sharply  initially  focused 


beams  generally  exhibit  the  development  of  only  a  single 
focus.  In  significant  contrast,  the  present  study  shows 
that  with  both  the  relativistic  and  charge-displacement 
mechanisms,  initially  focused  beams  also  generally  lead 
to  confined  modes  of  propagation.  Finally,  it  should  be 
noted  that  the  equation  for  relativistic  self-focusing  can 
be  considered  as  a  general  model  equation  of  propagation 
in  saturable  nonlinear  media. 

ACKNOWLEDGMENTS 

The  authors  acknowledge  fruitful  discussions  with  A. 
R.  Hinds,  R.  R.  Goldstein,  and  A.  McPherson.  Support 
for  this  research  was  provided  by  the  U.S.  Air  Office  of 
Scientific  Research,  the  U.S.  Office  of  Naval  Research, 
the  Strategic  Defense  Initiative  Organization,  the  Army 
Research  Office,  the  Department  of  Energy,  and  the  Na¬ 
tional  Science  Foundation  under  Grant  No.  PHY- 
9021265. 


‘Present  address:  Department  of  Physics,  University  of  Il¬ 
linois  at  Chicago,  Chicago,  Illinois  60680. 

[1] U.  Johann,  T.  Luk,  H.  Egger,  and  C.  K.  Rhodes,  Phys. 
Rev.  A  34,  1084(1986). 

[2]  T.  S.  Luk,  U.  Johann,  H.  Egger,  H.  Pummer,  and  C.  K. 
Rhodes,  Phys.  Rev.  A  32,  214  (1985);  K.  Boyer,  T.  S.  Luk, 
J.  C.  Solem,  and  C.  K.  Rhodes,  ibid.  39,  1 186  (1989). 

[3]  A.  1.  Akhiezer  and  R.  V,  Polovin,  Zh.  Eksp.  Teor.  Fiz.  30, 
915  (1956)  [Sov.  Phys.  JETP  3,  696  (1956)]. 

[4]  R.  Noble,  Phys.  Rev.  A  32,  460  (1985). 

[5]  A.  Magneville,  J.  Plasma  Phys.  44,  231  (1990). 

[6]  S.  V.  Bulanov,  V.  1.  Kirsanov,  and  A.  S.  Sakharov,  Fiz. 
Plasmy  16,  935  (1990)  [Sov.  J.  Plasma  Phys.  16,  543 
(1990)]. 

[7]  P.  Sprangle,  E.  Esarey,  and  A.  Ting,  Phys.  Rev.  Lett.  64, 
201 1  (1990);  P.  Sprangle,  C.  M.  Tang,  and  E.  Esarey,  Phys. 
Rev.  A  41,  4663  (1990);  A.  Ting,  E.  Esarey,  and  P.  Spran¬ 
gle,  Phys.  Fluids  B2,  1390  (1990). 

[8]  L.  D.  Landau  and  E.  M.  Lifshits,  The  Classical  Theory  of 
Fields  (Pergamon,  New  York,  1971). 

[9]  E.  S.  Sarachik  and  G.  T.  Schappcrt,  Phys.  Rev.  D  1,  2738 
(1970). 

[10]  J.  Kruger  and  M.  Bovyn,  J.  Phys.  A  9,  1841  (1976). 

[11]  D.  M.  Volkov.  Z.  Phys.  34,  250(1935). 

[12]  J.  N.  Bardsley,  B.  M.  Penetrantc,  and  M.  H.  Miitleman, 
Phys.  Rev.  A  40,  3823(1989). 

[13]  C.  Max,  J.  Arons,  and  A.  B.  Langon,  Phys.  Rev.  Lett.  33, 
209(1974). 

[14]  G.  Schmidt  and  W.  Horton,  Comments  Plasma  Phys. 
Controlled  Fusion  9,  85  (1985). 

[15]  H.  Hora,  Physics  of  Laser-Driven  Plasmas  (Wiley,  New 
York,  1981). 

[16]  P.  Sprangle  and  C.  M.  Tang,  in  Laser  Acceleration  of  Parti¬ 
cles  (the  Norton  Simon  Malibu  Beach  Conference  Center  of 
the  University  of  California,  Los  Angeles),  Proceedings  of 
the  Second  Workshop  on  Laser  Acceleration  of  Particles, 
edited  by  C.  Joshi  and  T.  Katsouleas,  AIP  Conf.  Proc.  No. 


130  (AIP,  New  York,  1985),  p.  156. 

[17]  P.  Sprangle,  C.  M.  Tang,  and  E.  Esarey,  IEEE  Trans.  Plas¬ 
ma  PS-15,  145(1987). 

[18]  Gou-Zheng  Sun,  E.  Ott,  Y.  C.  Lee,  and  P.  Guzdar,  Phys. 
Fluids  30,  526(1987). 

[19]  T.  Kurki-Suonio,  P.  J.  Morrison,  and  T.  Tajima,  Phys. 
Rev.  A  40,  3230(1989). 

[20]  E.  Esarey,  A.  Ting,  and  P.  Sprangle,  Appl.  Phys.  Lett.  53, 
1266(1988). 

[21]  W.  B.  Mori,  C.  Joshi,  J.  M.  Dawson,  D.  W.  Forslund,  and 
J.  M.  Kindel,  Phys.  Rev.  Lett.  60,  1298  (1988). 

[22]  P.  Gibbon  and  A.  R.  Bell,  Phys.  Rev.  Lett.  61,  1599  (1988). 

[23]  C.  J.  McKinstrie  and  D.  A.  Russel,  Phys.  Rev.  Lett.  61, 
2929(1988). 

[24]  A.  B.  Borisov,  A.  V.  Borovskiy,  V.V.  Korobkin,  A.  M. 
Prokhorov,  C.  K.  Rhodes,  and  O.  B.  Shiryaev,  Phys.  Rev. 
Lett,  65,  1753  (1990). 

[25]  J.  C.  Solem,  T.  S.  Luk,  K.  Boyer,  and  C.  K.  Rhodes,  IEEE 
J.  Quantum  Electron.  25,  2423  (1989). 

[26]  A.  B.  Borisov,  A.  V.  Borovskiy,  V.  V.  Korobkin,  A.  M. 
Prokhorov,  O.  B.  Shiryaev,  T.  S.  Luk,  J.  C.  Solem,  K. 
Boyer,  and  C.  K.  Rhodes  (unpublished);  A.  B.  Borisov,  A. 
V,  Borovskiy,  V.  V.  Korobkin,  A.  M.  Prokhorov,  O.  B. 
Shiryaev,  J.  C.  Solem,  K.  Boyer,  and  C.  K.  Rhodes  (un¬ 
published);  A.  B.  Borisov,  A.  V.  Borovskiy,  V.  V.  Korob¬ 
kin,  A.  M.  Prokhorov,  O.  B.  Shiryaev,  and  C.  K.  Rhodes, 
Kratkie  Soobsheniya  po  Fizike,  No.  9,  3  (1991)  (in  Rus¬ 
sian). 

[27]  A.  B.  Borisov,  A.  V.  Borovskiy,  V.  V.  Korobkin,  A.  M. 
Prokhorov,  O.  B,  Shiryaev,  X.  M.  Shi,  T.  S.  Luk,  A. 
McPherson,  J.  C.  Solem,  K.  Boyer,  and  C.  K.  Rhodes  (un¬ 
published). 

[28]  T.  S.  Luk  and  C.  K.  Rhodes,  Phys.  Rev,  A  38,  6180  (1988). 

[29]  G.  A.  Askaryan,  Zh.  Eksp.  Teor.  Fiz.  42,  1567  (1962)  [Sov. 
Phys.  JETP  15,  1088  (1962)]. 

[30]  T.  S.  Luk,  A.  McPherson,  K.  Boyer,  and  C.  K.  Rhodes, 
Opt.  Lett.  14,  1113  (1989). 


45 


relativistic  and  charge-displacement  self-  .  .  . 


5845 


[31]  A.  J.  Taylor,  C.  R.  Tallman,  J,  P.  Roberts,  C.  S.  Lester,  T. 
R.  Gosnell,  P.  H.  Y.  Lee,  and  G.  A.  Kyrala,  Opt.  Lett.  15, 
39(1990). 

[32]  P,  Maine,  D.  Strickland,  P.  Bado,  M.  Pessot,  and  G. 
Mourou,  IEEE  J.  Quantum  Electron.  24,  398  (1988). 

[33]  P.  Sprangle,  A.  Zigler,  and  E.  Esarey,  Appl.  Phys.  Lett. 
58,  345(1991). 

[34]  R.  Y.  Chiao,  E.  Garmire,  and  C.  H.  Townes,  Phys.  Rev. 
Lett.  13,  479  (1964). 

[35]  H,  A.  Haus,  Appl.  Phys.  Lett.  8,  128  (1966). 

[36]  A.  B.  Borisov,  A.  V.  Borovskiy,  V.  V.  Korobkin,  A.  M. 


Prokhorov,  O.  B.  Shiryaev,  and  C.  K.  Rhodes,  J.  Laser 
Phys.  1,  103(1991). 

[37]  V.  E.  Zakharov,  V.  V.  Sobolev,  and  V.  S.  Synakh,  Zh. 
Eksp.  Teor.  Fiz.  60,  136  (1971)  [Sov.  Phys.  JETP  33,  77 
(1971)]. 

[38]  A.  B.  Borisov,  A.  V.  Borovskiy,  V.  V.  Korobkin,  C.  K. 
Rhodes,  and  O.  B.  Shiryaev  (unpublished);  A.  B.  Borisov, 
A.  V.  Borovskiy,  V.  V.  Korobkin,  C.  K.  Rhodes,  and  O.  B. 
Shiryaev,  in  Short  Wavelength  Lasers  and  Their  Applica- 
tionSy  edited  by  V.  V.  Korobkin  and  M.  Yu.  Romanovskiy 
(Nova  Science,  Commack,  NY,  1992),  p.  261, 


Appendix  F; 


"Observation  of  Relativistic  and  Charge-Displacement  Self-Channeling  of  Intense 
Subpicosecond  Ultraviolet  (248  nm)  Radiation  in  Plasmas" 


61 


Volume  68,  Number  15 


PHYSICAL  REVIEW  LETTERS 


13  APRIL  1992 
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Experimental  studies  examining  a  new  relativistic  regime  of  high-intensity  short-pulse  propagation  in 
plasmas  have  been  performed  which  present  evidence  for  the  formation  of  a  stable  mode  of  spatially 
confined  (channeled)  propagation.  For  an  electron  density  of  —  L35x|0^'  cm“^  and  a  power  of 
—  3x  10"  W,  the  results  indicate  a  channel  radius  <  I  pm  and  a  peak  intensity  —10*’  W/cml  Com¬ 
parison  of  these  findings  with  a  dynamical  theory  yields  agreement  for  both  the  longitudinal  structure 
and  the  radial  extent  of  the  propagation  observed. 

PACS  numbers:  52.40.Db.  42.25.-p,  42.65.Jx 


A  fundamentally  new  regime  of  electromagnetic  propa¬ 
gation  is  expected  to  arise  in  plasmas  for  short-pulse  radi¬ 
ation  at  sufficiently  high  intensity.  Calculations  of  the 
propagation  in  plasmas,  incorporating  both  relativistic 
Il,2l  and  charge-displacement  mechanisms  [3-7],  indi¬ 
cate  that  the  combined  action  of  these  processes  can  lead 
to  a  new  stable  form  of  spatially  confined  (channeled) 
propagation.  This  Letter  (1)  reports  the  results  of  the 
first  experimental  study  probing  the  physical  regime 
relevant  to  the  observation  of  relativistic  and  charge- 
displacement  self-channeling  and  (2)  presents  the  initial 
comparison  of  these  experimental  findings  with  matching 
theoretical  calculations  performed  with  the  computation¬ 
al  procedures  described  in  Ref.  I7l. 

The  experimental  arrangement  used  in  these  studies  is 
illustrated  in  Fig.  1(a).  The  source  of  radiation  was  a 
subpicosecond  KrF*  (X— 248  nm)  laser  that  has  been  de¬ 
scribed  elsewhere  [8l.  It  delivered  a  linearly  polarized 
power  of  — 3xl0“  W  (—150  mJ,  pulse  duration  —500 
fs)  in  a  beam  with  a  diameter  of  —42  mm.  When  this 
radiation  was  focused  into  the  chamber  with  lens  L 1 
(//7),  a  focal  radius  ro*^3.5  pm  was  measured,  giving  a 
maximum  intensity  /o-^8.6x|0*^  W/cm^  The  medium 
was  provided  by  filling  the  chamber  statically  with  gas 
[He,  Ne,  Ar,  Kr,  Xe,  N2,  CO2,  or  a  mixture  of  Xe  (4%) 
and  N2  (96%)]  up  to  a  maximum  density  of  —  1.89 x  I0^° 
cm 

The  diffracted  248-nm  radiation  was  measured  as  a 
function  of  the  angle  (0)  with  respect  to  the  direction  of 
the  incident  radiation.  The  incident  laser  beam  was 
blocked  by  a  metal  disk  on  the  output  window  of  le 
chamber  and  lens  L2  imaged  the  region  near  the  focal 
zone  on  a  fluorescent  screen  5.  The  diaphragm  D  in  front 
of  lens  L2  restricted  the  collection  of  the  diffracted  light 
to  a  solid  angle  of  ^5®  while  simultaneously  increasing 
the  depth  of  field.  The  angle  between  the  axis  of  the  lens 
L2  and  the  axis  of  the  incoming  laser  radiation  could  be 
readily  varied  up  to  a  maximum  angle  of  6—15®.  Two 
flat  mirrors  coated  for  high  reflection  (—99%)  at  248 
nm,  both  having  a  spectral  bandwidth  of  — 10  nm,  served 


in  reflection  as  spectral  filters  (F)  for  the  diffracted  laser 
radiation  so  that  only  the  scattered  248-nm  radiation 
could  illuminate  the  screen.  An  attenuator  A  was  em¬ 
ployed  to  adjust  the  intensity  on  the  screen  and  the  im¬ 
ages  formed  were  recorded  through  the  visible  fluores- 


Axial  Distance  z  (nm) 

FIG.  I.  (a)  Experimental  apparatus  used  in  studies  of  propa¬ 
gation.  See  text  for  description,  (b)  Data  concerning  the  pat¬ 
tern  of  propagation  observed  with  a  single  pulse  in  N2  at  a  den¬ 
sity  of  —1.35x10^  cm“^  The  maximum  intensity  is  half  the 
detector  (CCD)  saturation.  The  radiation  is  incident  from  the 
left.  Inset:  Photographic  data  with  a  vertical  spatial  resolution 
of  —10  /im.  The  graph  illustrates  the  one-dimensional  axial 
profile  taken  along  the  direction  of  propagation  (z)  of  the  pho¬ 
tographic  data  (inset).  The  spacing  of  the  maxima,  5  —  200 
±  20  /im,  is  indicated. 


Volume  68,  Number  15 


PHYSICAL  REVIEW  LETTERS 


13  APRIL  1992 


cence  produced  with  a  microscope  and  a  charge-coupled- 
device  (CCD)  camera. 

The  characteristic  behavior  observed  is  well  illustrated 
by  the  data  recorded  with  N2.  The  measured  result, 
shown  in  Fig.  1(b),  corresponds  to  a  density  pn,“^L35 
X  10“®  cm  To  the  left,  in  the  photographic  inset,  a  rel¬ 
atively  large  cone  of  light  Rayleigh  scattered  from  the 
plasma  is  visible  at  all  angles  as  the  energy  propagates  to¬ 
ward  the  focal  point  of  the  lens,  while,  in  the  region  to  the 
right  of  the  conical  apex,  a  narrow  filament  developed. 
The  diameter  of  this  filament  is  not  greater  than  10  /im, 
the  measured  spatial  resolution  of  the  imaging  system. 
The  distribution  of  intensity  observed  along  the  filament 
exhibited  several  bright  features  attributed  to  diffraction 
because  they  could  not  be  seen  for  0>2O°.  Since  the 
axis  of  the  imaging  lens  corresponded  to  an  angle 
©■“7.5°,  the  scale  along  the  abscissa  of  the  photographic 
data  is  reduced  by  almost  eightfold,  giving  the  maximum 
length  of  the  filament  as  —1  mm.  The  graph  in  Fig.  1(b) 
represents  a  one-dimensional  axial  profile,  taken  along 
the  direction  of  propagation  (z),  of  the  observed  intensity 
pattern  (inset).  Three  peaks  (a,^,  y)  are  visible  with  a  a 
spatial  separation  of  5  “200  ±20  pm.  The  normal  Ray¬ 
leigh  range  for  the  focal  geometry  used  was  --200  pm. 

The  diameter  of  the  filamentary  channel  is  an  impor¬ 
tant  dynamical  variable,  which  we  estimated  by  measur¬ 
ing  the  maximum  angular  deviation  of  the  diffracted 
light.  The  experimental  value  (p  of  this  diffracted  cone 
was  0  —  20°,  a  magnitude  indicating  a  radius  0.9  pm 
though  the  relation  0“I.22X./r^  Filaments  of  this  gen¬ 
eral  nature  were  observed  at  densities  above  —  1.35  x  10^® 
cm  in  N2,  Ne,  Ar,  Kr,  CO2,  and  a  mixture  of  Xe  (4%) 
and  N2  (96%),  but  not  in  He  and  Xe,  two  materials  dis¬ 
cussed  further  below. 

Two  mechanisms  exist  that  could  modify  the  refractive 
index  of  the  medium  and  lead  to  the  observed  behavior. 
They  are  (I)  the  Kerr  effect  stemming  from  the  ions  and 
(2)  the  relativistic  and  charge-displacement  process  [7]. 
Since  the  pulse  duration  is  very  short  (  —  500  fs),  the 
motion  of  the  ions  is  negligible  [3l,  and  no  contribution 
can  arise  from  expulsion  of  the  plasma  from  the  high- 
intensity  zone.  An  implication  of  the  estimate  of  the 
channel  radius  (r^— 0.9  pm)  is  that  the  observed  propa¬ 
gation  is  associated  with  intensities  in  the  10 **-10” 
W/cm“  range.  Under  these  conditions,  available  experi¬ 
mental  evidence  [9,101  on  multiphoton  ionization  indi¬ 
cates  that  He  should  be  fully  ionized  and  the  C,  N,  and  O 
atoms  constituting  the  molecular  materials  would  retain, 
at  most,  only  \s  electrons. 

Consider  explicitly  the  case  of  N2,  which  has  estimated 
[9,10]  threshold  intensities  for  the  production  of 

and  of  l.6x|0'®,  6.4x|0**,  and  1.3x|0” 
W/cm“,  respectively.  Hence,  the  volume  of  the  channel 
would  be  largely  ionized  to  with  certain  localized 

high-intensity  regions  contributing  some  Two 

consequences  of  this  pattern  of  ionization  follow,  namely, 
(i)  the  Kerr  effect  arising  from  the  ions  is  small,  since  the 


polarizabilities  of  the  remaining  I5  electrons  are  low,  and 
(ii)  the  electron  density  (n^ )  initially  produced  in  the  fo¬ 
cal  region  is  nearly  uniform.  Therefore,  =  I.35x  10“‘ 
cm  for  the  data  on  N2  shown  in  Fig.  1  (b). 

A  critical  power  Per  for  self-channeling,  arising  from 
the  relativistic  and  charge-displacement  mechanism,  can 
be  defined  [11]  as 

Per  *  Ve  “)  go  (r)r  dr((o/(Op^o)  ‘ 

=  1.62xl0'®(ncr/ne.)  W,  (1) 

where  m*.,  e,  and  c  have  their  customary  identifications,  (o 
is  the  laser  angular  frequency,  ojp^o  is  the  plasma  frequen¬ 
cy  for  the  uniform  unperturbed  plasma  with  electron  den¬ 
sity  Her  is  the  critical  electron  density  (for  X  “248  nm, 
ngr“  1.82X  10^^  cm“^),  and  goir)  is  the  Townes  mode 
[121. 

The  critical  powers  associated  with  the  experimental 
conditions,  for  He  and  N2  at  a  medium  density 
p“L35xiO^®  cm“\  are  1.08x10*2  and  2,19x10**  W, 
respectively.  Therefore,  since  the  incident  power  was 
p=3x  10**  W,  no  filament  was  expected  in  He,  a  predic¬ 
tion  conforming  with  the  observation  of  none.  Moreover, 
the  diffracted  cone  of  radiation  was  also  absent  with  He. 
in  contrast,  P/Per— 1-37  for  N2,  a  condition  that  held 
generally  (P/Pcr>l)  for  all  materials  which  exhibited 
evidence  for  channel  formation.  We  note,  however,  that 
some  contribution  from  the  Kerr  effect  may  be  present, 
even  for  the  light  materials  (Ne,  N2,  and  CO2),  in  the 
early  stage  of  channel  formation  prior  to  the  development 
of  a  substantial  level  of  ionization,  and  that  the  heavier 
gases  (Ar,  Kr,  and  Xe)  may  involve  a  more  significant 
influence  from  the  Kerr  process.  A  specific  estimate  of 
the  nonlinear  index  change  arising  from  both  N^"**  and 
N®^  at  an  intensity  of  —10*’  W/cm^  indicates  that  their 
contribution  is  less  than  10”^  that  of  the  free  electrons; 
hence  the  ionic  contribution  can  be  neglected  in  N2  for 
the  conditions  studied. 

A  direct  comparison  will  now  be  made  between  the 
theoretical  analysis,  fully  described  in  Ref.  [7],  and  the 
experimental  findings  for  N2.  This  comparison  can  be  ac¬ 
complished  for  both  the  longitudinal  intensity  profile  and 
the  radial  extent  of  the  channel.  Figure  2(a)  illustrates 
the  intensity  profile  /(r,z)//o  calculated  with  physical  pa¬ 
rameters  corresponding  to  those  of  the  experiment  (i.e., 
P  — 3x|0"  W,  ro“3.5  pm,  n*.  “  1.35x  10^*  cm“\  and 
P/Pcr*1.37).  Importantly,  all  of  these  parameters  are 
based  on  independent  measurements  of  (1)  the  laser  pulse 
involving  determinations  of  the  energy  and  power  (P), 
(2)  the  focal  radius  (ro)  of  the  incident  radiation,  and  (3) 
the  characteristics  of  the  multiphoton  ionization  [lOl  gen¬ 
erating  the  electron  density  (n^).  Therefore,  this  compar¬ 
ison  does  not  involve  a  fit  with  a  free  parameter.  The 
normalized  electron  density  calculated  is  presented  in 
Fig.  2(b),  from  which  it  is  seen  that  electronic  cavitation 
occurs  only  near  the  positions  of  the  maxima  in  the  inten¬ 
sity  profile  [Fig.  2(a)l.  Curve  A  in  Fig.  2(c)  represents 
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the  one-dimensional  axial  intensity  profile  /(0.r)//o  cor¬ 
responding  to  the  calculated  distribution  shown  in  Fig. 
2(a).  The  spacing  (S)  of  the  maxima  is  seen  to  be 
5~I85  fim,  a  value  in  close  agreement  with  the  experi¬ 
mental  figure  (5-200±20 //m)  illustrated  in  Fig.  Kb). 
Furthermore,  analysis  has  shown  that  the  spacing  S  is 
quite  sensitive  to  the  power  Po  and  electron  density  n,., 
particularly  in  the  region  close  to  the  threshold  (see  Eq. 
(I)].  With  respect  to  the  results  illustrated  in  Fig.  2(c), 
an  increase  in  n,  by  less  than  5%  causes  a  reduction  in 
the  spacing  S  by  approximately  25%.  Therefore,  substan¬ 
tial  changes  in  the  physical  parameters  would  grossly 
alter  the  comparison  of  the  experimental  and  theoretical 
results. 

Theoretical  studies  (3-7,11]  indicate  that  the  charge 
displacement  plays  a  very  important  dynamical  role.  In 
order  to  test  this  hypothesis,  identical  calculations  were 
made  for  N?  for  the  purely  relativistic  case  (2l  which  ex¬ 
plicitly  neglects  the  charge-displacement  term,  namely, 
elimination  of  the  term  (c^/otloroU±(l+///r)''^  in  Eq. 
(24)  of  Ref.  [7],  Significantly,  the  resulting  axial  profile 
(curve  B  in  Fig.  2(c)]  exhibits  only  a  stng/e  relatively 
weak  maximum,  for  0:<  z  <  600  /rm,  an  outcome  sharply 
at  variance  with  both  the  full  theoretical  analysis  and  the 
experimental  observation.  Although  the  expected  charge 
displacement  is  highly  localized  (Fig.  2(b)],  this  compar¬ 
ison  reveals  the  strong  influence  it  has  on  the  propagation 
(3,7,11].  At  a  greater  incident  power  10).  a 

continuous  channel  in  the  electron  distribution  is  expect¬ 
ed  to  develop  17]. 

The  measurements  indicated  an  approximate  value  of 
r,=0.9  fjm  for  the  radial  extent  of  the  channel,  a  result 
that  can  be  compared  with  the  corresponding  theoretical 
figure.  Figure  2(d)  illustrates  five  radial  intensity  profiles 
/(r,r,)//o  of  the  distribution  pictured  in  Fig.  2(a).  Since 
the  measurement  of  this  angularly  scattered  radiation  did 
not  correspond  to  a  known  longitudinal  position,  this 
comparison  can  only  be  qualitative,  but  the  radial  distri¬ 
butions  shown  indicate  that  the  expected  value  lies  in  the 
interval  0.5<r:<  1.0  /im,  a  range  that  comfortably  in¬ 
cludes  the  experimental  value  r,. 

The  results  observed  with  Xe  deserve  additional  discus¬ 
sion,  since  those  experiments  did  not  give  evidence  for  the 
formation  of  a  channel.  In  significant  contrast  to  the  case 
involving  Ni,  the  electron  density  n,  produced  by  the 
multiphoton  ionization  (10]  in  Xe  is  expected  to  be  very 

FIG.  2.  Calculations  for  Ni  with  /»-3x  10"  W,  ro-3.5  pm, 
/I, -I.35X  10-'  cm”’,  and  /o-8.6x|0"  W/cm\  (a)  Normal¬ 
ized  intensity  l(r,:)/lo.  (b)  Normalized  electron  density 

N(r,z)/Nn  for  N;  with  No^n,.  (c)  Normalized  one-dimen¬ 
sional  axial  intensity  profiles  /(0,c)//o.  Curve  A.  full  theory  for 
data  in  panel  (a).  tf—ISS  pm.  Curve  B.  calculation  with 
charge-displacement  term  neglected,  (d)  Normalized  radial  in¬ 
tensity  profiles  Hr.Zi)/lo  corresponding  to  panel  (a).  Longitu¬ 
dinal  positions  ;i-l72  pm.  r:-245  pm,  r.i-358  pm.  r4-44l 
m.  and  r5“559  pm  and  r,“0.9  pm. 
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nonuniform  spatially.  For  intensities  spanning  10'‘-10'* 
W/cm^  the  corresponding  density  n,  would  vary  by  over 
a  factor  of  2.  Since  this  nonuniformity  would  tend  to 
reduce  Ihe  refractive  index  locally  in  the  central  high- 
intensity  region,  a  significant  defocusing  action  is  expect¬ 
ed  which  could  suppress  the  channel  formation. 

Finally,  we  note  (1)  that  the  intensity  distribution  is 
not  expected  to  depend  strongly  on  the  state  of  polariza¬ 
tion  113,14]  and  (2)  that  losses  to  the  plasma  may  be 
significant,  particularly  at  electron  densities  close  to 

The  first  experiments  examining  a  new  relativistic  re¬ 
gime  of  high-intensity  pulse  propagation  in  plasmas  have 
been  performed  and  the  findings  indicate  the  formation  of 
a  channeled  mode  of  propagation  over  a  length  consider¬ 
ably  greater  than  the  Rayleigh  range.  Specific  compar¬ 
isons  of  the  experimental  observations  with  a  dynamical 
theory,  which  explicitly  includes  both  the  influence  of  the 
relativistic  mass  shift  and  the  displacement  of  the  elec¬ 
tronic  component  of  the  plasma,  produce  excellent  agree¬ 
ment  for  both  the  longitudinal  structure  of  the  intensity 
profile  and  the  radial  extent  of  the  channel.  While  the 
present  channel  contains  several  foci,  a  continuous  chan¬ 
nel  is  predicted  to  develop  at  higher  power.  Finally,  the 
intrinsically  very  high  concentration  of  power  associated 
with  this  mechanism  of  channeled  propagation  provides 
an  efficient  and  general  method  for  the  production  of  con¬ 
ditions  necessary  for  x-ray  amplification  [15]. 
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Appendix  G:  "Ultrahigh  Power  Compression  for  X-Ray  Amplification:  Multiphoton  Cluster 
Excitation  Combined  with  Nonlinear  Channeled  Propagation" 
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ABSTRACT 


The  ability  to  apply  power  densities  controllably  at  or  above  vigorous  thermonuclear 
levels  (>  10^^  W/cm^)  in  materials  is  the  basic  issue  for  achieving  efficient  amplification  of 
X-rays.  Recent  experimental  and  theoretical  findings  concerning  (1)  the  multiphoton  production 
of  X-rays  from  clusters  and  (2)  high-intensity  modes  of  channeled  propagation  in  plasmas 
indicate  an  entirely  new  method  for  producing  the  conditions  necessary  for  strong  amplification 
in  the  multi-kilovolt  range.  These  two  new  nonlinear  phenomena  are  being  united  to  produce 
and  control  the  imperative  power  compression;  the  multiphoton  mechanism  serves  to  establish 

t 

the  condition  locally  while  the  confined  propagation  provides  the  required  spatial  organization. 
The  present  work,  which  experimentally  demonstrates  the  first  combined  expression  of  these  two 
complex  nonlinear  processes  through  direct  X-ray  imaging  of  Xe(M)  emission  (-1  keV)  in 
stable  self-trapped  channels,  (a)  reveals  the  exceptional  compatibility  of  their  mutual  scaling  for 
realizing  the  necessary  power  density,  03)  provides  confirming  evidence  for  the  action  of  a 
superstrong  coherent  multi-electron  intense-field  interaction  in  the  X-ray  generation  from  the  Xe 
clusters,  and  (y)  furnishes  new  detailed  information  on  the  dynamics  of  the  radial  intensity 
distributions  associated  with  the  channeled  propagation.  The  resulting  knowledge  of  the  scaling 
relations  underlying  these  phenomena  enables  the  optimum  conditions  for  amplification  to  be 
specified  up  to  a  quantum  energy  of  —  5  keV.  The  harmonious  use  of  these  new  nonlinear 
processes  is  expected  to  lead  to  an  advanced  generation  of  extraordinarily  bright  X-ray  sources 
in  the  multi-kilovolt  region  having  a  peak  brightness  of  —  10^^-1(P^  7S“^  (mrad)"^  (mm)"^ 
(0. 1  %  a  level  sufficient  for  biological  holographic  imaging  capable  of  providing  a  high 

resolution  visualization  of  the  molecular  anatomy  of  cells,  tissues,  and  organisms  in  the  natural 
state. 
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Stable  channeled  propagation  of  subpicosecond  (248  nm)  radiation  arising  from 
relatiyistic-charge-displacement  nonlinearities  [1]  has  been  observed  [2,3]  in  underdense  plasmas. 
Experimental  tests  [4]  of  a  concept  [5]  for  a  molecular  design  associated  with  multiphoton  X-ray 
production  induced  by  a  strong  nonlinear  interaction  have  shown  that  copious  Xe(M),  Kr(L),  and 
Xe(L)  emission  can  be  generated  with  corresponding  Kr  and  Xe  clusters  in  the  1-5  keV  region 
[4,6-8].  The  experiments  demonstrated  that  high  X-ray  yields  could  be  readily  produced  from 
gaseous  targets  and  estimates  stemming  from  the  observed  X-ray  spectra  indicated  the 
achievement  of  maximal  atom-specific  energy  transfer  rates  of  ~  0. 1-1.0  W/atom. 

The  present  work  (1)  demonstrates  a  compatible  unification  of  the  physical  conditions 
governing  the  combined  expression  of  these  two  highly  nonlinear  phenomena,  (2)  provides  data 
corroborating  the  presence  of  a  superstrong  multiphoton  interaction  in  the  X-ray  production  from 
the  clusters  [6-9],  and  (3)  gives  important  new  information  on  the  properties  of  the  radial 
intensity  profiles  produced  by  the  nonlinear  propagation.  The  outcome  is  a  congruous  picture 
of  the  mutual  scaling  of  these  processes  that  enables  a  general  specification  of  the  optimum 
conditions  necessary  for  amplification  up  to  a  quantum  energy  of  ~  5  keV. 

Calculations  have  shown  that  the  observed  characteristics  of  the  channeled  propagation 
[1]  can  be  put  in  close  correspondence  [2,3]  with  a  theoretical  picture  incorporating  relativistic 
[10,11]  and  charge-displacement  [12-18]  mechanisms.  In  particular,  those  studies  indicated  that 
intensities  as  high  as  —  10^®  W/cm^  were  dynamically  produced  in  the  narrow  channels  formed 
[2,3],  a  value  sufficient  for  the  copious  multiphoton  generation  of  multi-kilovolt  X-ray  emission 
from  Kr  and  Xe  clusters  [4-8].  Therefore,  the  use  of  a  medium  composed  of  Kr  or  Xe  clusters 
presents  the  opportunity  for  a  study  combining  the  dynamics  of  the  confined  propagation  and  the 
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properties  of  the  nonlinear  cluster  excitation  through  direct  X-ray  imaging  of  the  high-intensity 
zones  produced  by  the  channeling.  We  describe  below  the  results  obtained  by  using  Xe(M) 
emission  (-I  keV)  generated  from  Xe  clusters  in  channels  produced  with  subpicosecond  248 
nm  radiation. 

The  experimental  arrangement  used  to  conduct  spectroscopic  studies  of  rare  gas  clusters 
has  been  described  elsewhere  [6].  The  laser  [19]  used  had  a  wavelength  of  248  nm,  a  pulse 
width  of  —  300  fs,  a  power  of  —  0.7  TW,  and  was  focused  to  a  maximum  intensity  of  ~  10^^ 
W/cm"^  with  an  f/3  off-axis  parabolic  mirror  in  a  gas  target.  This  target  was  hydrodynaAically 
produced  by  a  cooled  high-pressure  pulsed-valve  fitted  with  a  circular  sonic  nozzle  having  a 
diameter  of  ~  0.5  mm.  In  order  to  eliminate  distortion  and  absorption  of  the  incident  beam  as 
it  approached  the  focal  zone  through  the  outer  regions  of  the  material  produced  by  the  nozzle 
[20],  a  barrier  with  a  ~  100  nm  hole  was  positioned  near  the  opening  of  the  nozzle  as  shown 
in  the  inset  in  Fig.  2(b),  The  focus  of  the  f/3  mirror  was  then  placed  at  the  position  of  the 
opening  in  the  barrier  so  that  the  initial  conditions  of  medium  density  and  laser  intensity 
corresponding  to  the  launching  of  the  wave  were  well  defined. 

Previous  spectroscopic  studies  [6]  have  demonstrated  that  bright  Xe(M)  emission 
spanning  the  9.5-16  A  region  is  readily  produced  from  Xe  clusters.  In  order  to  record  the 
spatial  pattern  of  the  propagation,  this  emission  was  observed  with  an  X-ray  "pinhole"  (CCD) 
camera.  The  camera  was  an  Hectrum  EDCIOOO  using  a  video  chip  having  a  196  x  165  pixel 
array  with  an  individual  pixel  size  of  13.5  /tm  x  16  /xm.  The  X-ray  images  were  obtained  by 
viewing  the  region  excited  by  the  focused  248  nm  radiation  directly  below  the  pulsed-jet  in  a 
direction  transverse  to  the  propagation  of  the  248  nm  beam.  For  images  of  the  Xe(M)  emission, 
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the  camera  was  equipped  with  a  filter  with  the  composition  (2000  A  A1  +  1  /xm  polycarbonate 
+  10  /xm  Be)  and,  consequently,  detected  all  radiation  with  a  quantum  energy  ^  700  eV. 
When  the  camera  was  equipped  with  a  ~  25  /xm  pinhole  and  positioned  to  give  a  magnification 
factor  of  1.18,  the  brightness  of  the  Xe(M)  emission  was  sufficiently  high  that  recordings  of 
single-shot  exposures  could  be  obtained.  Given  the  magnification  used  and  the  geometry  of  the 
pixel,  the  pinhole  camera  has  a  predicted  limiting  spatial  resolution  A  ~  30  /xm,  a  property 
verified  by  the  observed  characteristics  of  the  images  described  below. 

Examples  of  the  X-ray  images  recorded  for  two  different  densities  of  the  Xe  medihm  are 
illustrated  in  Figs.  1(a)  and  2(a),  respectively.  The  former  shows  a  strongly  longitudinally 
modulated  distribution  while  the  latter  demonstrates  that  a  highly  uniform  distribution  can  be 
produced  under  appropriate  conditions.  The  corresponding  axial  X-ray  intensity  distributions 
(radius  r  =  0  line-outs  of  these  data),  which  are  represented  in  Figs.  1(b)  and  2(b),  are  found 
to  be  in  good  correspondence  with  the  matching  theoretical  distributions  shown  in  Figs.  1(c)  and 
2(c).  Significantly,  as  in  the  earlier  study  [2],  all  the  parameters  involved  in  the  computed 
distributions  are  based  on  quantities  derived  from  independent  measurements,  so  that  the 
comparisons  with  the  experimental  results  are  sharply  defined  and  do  not  involve  a  fit  with  any 
free  parameters.  Importantly,  the  conditions  of  these  experiments  also  corresponded  to  the  zone 
of  stable  confined  propagation  previously  found  in  an  analysis  of  the  response  of  the  channeled 
power  distribution  to  small  azimuthal  perturbations  [16,18],  a  point  shown  explicitly  below. 
Furthermore,  all  the  computed  results  illustrated  in  Figs.  1  and  2  allow  for  the  variation  of  the 
density  of  the  medium  caused  by  the  hydrodynamic  expansion  from  the  nozzle  [21]. 

Of  chief  interest  in  the  present  context  is  the  radial  dependence  of  the  X-ray  images, 
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since  the  behavior  of  the  transverse  dimension  of  the  emitting  region  represents  a  sensitive  probe 
of  both  the  dynamics  of  the  propagation  and  the  mechanism  of  excited  state  production  in  the 
clusters.  Consider  the  image  shown  in  Fig.  2(a),  which  has  a  transverse  width  of  —  120  /xm 
at  its  midpoint  (z  =  1050  /xm).  This  width  is  considerably  greater  than  the  resolution  (A  ~  30 
fim)  of  the  camera  and  both  the  observed  and  calculated  values  [  ~  2  ^m]  associated  with  the 
confined  propagation  of  the  248  nm  radiation  [2,3].  Consequently,  the  measurement  shows 
unambiguously  that  X-rays  are  being  produced  well  outside  the  small  central  core  of  the  channel 
that  confines  the  propagation  of  the  248  nm  energy.  In  addition,  the  image  shown  in  Fi|.  1(a) 
clearly  demonstrates  that  the  radial  extent  of  the  X-ray  emission  can  have  a  significantly 
structured  longitudinal  (z)  dependence. 

The  existence  of  a  bright  emitting  zone  having  a  diameter  of  ~  100  /xm  and  the  observed 
longitudinal  signature  in  the  radial  structure  is  inconsistent  with  either  (a)  collisional  excitation 
arising  from  energetic  electrons  produced  in  the  core  of  the  channel  or  (b)  recombination  into 
excited  states  of  highly  ionized  ions  which  migrate  radially  from  the  central  region.  At  the 
outset  we  see  that  the  largest  observed  radial  dimension  in  Fig.  2(a)  is  greater  than  the 
corresponding  maximum  value  in  Fig.  1(a),  although  the  density  of  the  medium  in  the  former 
is  significantly  higher  than  in  the  latter.  This  trend  is  exactly  the  opposite  of  that  expected  for 
the  range  of  energetic  charged  particles,  either  electrons  or  ions,  in  matter.  Furthermore,  if  the 
origin  of  the  X-ray  emission  were  to  arise  from  collisional  excitation  by  electrons  expelled  from 
the  channel  through  the  effect  of  the  ponderomotive  potential,  then  it  should  be  possible  to  relate 
the  observed  radial  size  of  the  emitting  region  to  an  estimate  of  the  maximum  range  of  electrons 
[22]  in  materials  with  the  known  range-energy  [23]  relationships.  Specifically,  since  the  data  in 

6  A.  B.  Borisov  et  al. 


72 


Figs.  1(c)  and  2(c)  show  that  the  peak  intensity  in  the  channel  is  expected  to  fall  in  the  10^^- 
lO^O  -W/cm^  region,  for  electrons  having  kinetic  energies  on  the  scale  of  the  ponderomotive 
potential  (-  100-5(X)  keV),  the  expected  range  would  be  -  100  mg/cm^.  The  corresponding 
distance  would  exceed  several  centimeters  for  the  prevailing  experimental  conditions,  a  value 
in  clear  conflict  with  the  data.  Of  course,  only  a  rather  small  fraction  of  electrons  would  be 
expected  to  behave  in  this  fashion,  since  the  energy  associated  with  the  charge  separation  would 
otherwise  be  too  great. 

Energetic  ions  can  also  be  ruled  out  as  the  source  of  the  extended  region  of  emission. 
It  can  be  readily  seen  from  data  on  ion  ranges  in  materials  [24]  that  the  Xe  ion  kinetic  energy 
necessary  to  produce  the  observed  depth  of  penetration  (~  100  /ig/cm^)  is  —  0.5-1. 0  MeV. 
This  required  energy  is  far  above  the  value  that  could  be  produced  by  the  coulomb  explosion  of 
a  cluster  under  our  experimental  conditions.  For  example,  the  maximum  kinetic  energy  of  Xe 
ions  generated  in  the  disintegration  of  a  Xej25  cluster  irradiated  at  —  10^^  W/cm^  would  be 
^  60  keV. 

The  remaining  possibility  is  direct  excitation  of  the  clusters  throughout  the  observed 
radiating  zone,  a  mechanism  requiring  fulfillment  throughout  that  region  of  the  conditions 
specified  by  the  allowed  zones  governing  prompt  X-ray  emission  [4-6].  The  threshold  intensity 
for  the  generation  of  Xe(M)  emission  from  Xe  clusters  has  been  independently  determined 
experimentally  [6]  to  be  ~  3  x  10^^  W/cm^.  Significantly,  this  experimental  value  was  found 
to  be  substantially  lower  than  the  minimum  value  predicted  (~  4  x  10  ”  W/cm**^)  with  the 
original  formulation  of  the  multiphoton  cluster  interaction  [5],  a  conclusion  that  followed 
independent  of  the  distribution  of  cluster  sizes  involved  in  the  study  [6,9].  Since  the  full 
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ponderomotive  potential  of  a  free  electron  is  considerably  less  than  100  eV  at  an  intensity  of  ~ 
3  X  10^-^  W/cm*^  for  irradiation  with  a  wavelength  of  248  nm,  the  available  energy  of  a  single 
electron  is  well  below  the  threshold  energy  [25,26]  (M5  -  675  eV)  needed  for  M-shell 
ionization  in  Xe.  However,  a  generalization  of  the  cluster  interaction,  incorporating  an 
enhancement  in  the  coupling  strength  for  inner-shell  excitation  arising  from  the  coherent  motion 
[9,27,28]  of  the  (Z)  field-ionized  [29]  electrons  induced  by  the  external  driving  field,  reduced 
the  required  threshold  intensity  to  -  3.8  x  10^^  W/cm^,  a  value  agreeing  with  the  experimental 
determination  within  the  uncertainty  of  measurement  [6].  Physically,  the  extension  ’of  the 
allowed  region  to  lower  intensities  arises,  because  the  Z  (field-ionized)  coherently  driven 
electrons,  individually  having  a  mass  m,  act  like  a  auasi-particle  [6,8,9,27]  with  mass  Zm  and 
charge  Ze.  In  this  situation,  the  single  particle  energies  add  [27]  so  that  the  required  excitation 
energy  [5]  of  the  inner-shell  state  Cg  can  be  shared  in  the  interaction.  Consequently,  this  feature 
of  the  interaction  leads  to  a  Z-fold  reduction  in  the  necessary  single-particle  energy,  thereby, 
substantially  extending  the  allowed  zone  by  lowering  the  threshold  condition  [6]  to  the  boundary 

f 

Ig  shown  in  Fig.  3.  Since  the  radial  extent  of  the  X-ray  emission  from  the  channels  can  be 
directly  related  to  this  threshold  intensity  (Iq),  the  minimum  value  of  which  is  designated  as  I^/ 
in  Fig.  3,  images  of  the  emitting  region  can  give  critical  information  on  the  existence  and 
properties  of  the  coherent  interaction. 

An  evaluation  of  the  radial  distance  over  which  conditions  exist  corresponding  to  the 
allowed  zone  for  Xe(M)  emission  requires  knowledge  of  the  radial  profile  of  the  intensity 
distribution  of  the  248  nm  radiation  supported  by  the  channeled  propagation.  The  characteristic 
form  of  this  radial  intensity  distribution,  for  conditions  corresponding  to  the  Xe(M)  emission 

8  A.  B.  Borisov  et  al. 


74 


shown  in  Fig.  2(a)  at  z  =  1050  /xm,  is  illustrated  in  Fig.  4.  This  distribution,  evaluated  as  a 
function  of  longitudinal  position  (z),  can  be  used  to  determine  the  radial  locus  of  the  threshold 
intensity  for  Xe(M)  emission.  It  is  important  to  note  that,  since  the  confined  248  nm  radiation 
propagates  in  a  channel  that  acts  as  a  dielectric  waveguide,  a  substantial  field  can  exist  at  radial 
positions  far  exterior  to  the  narrow  high-intensity  core.  This  important  characteristic  is  plainly 
exhibited  by  the  radial  intensity  profile  illustrated  in  Fig.  4.  Indeed,  for  a  threshold  value  of  I^ 
—  3  X  10^^  W/cm^,  Fig.  4  indicates  a  maximum  radius  of  ~  27  /xm,  a  value  more  than  ten¬ 
fold  the  radius  of  the  central  high  intensity  zone  [2,3].  This  radial  characteristic  ’of  the 
propagating  field  greatly  enlarges  the  effective  volume  of  the  channel  for  X-ray  generation  and, 
consequently,  is  a  crucial  feature  of  the  joint  action  of  these  two  nonlinear  processes. 

Three  significant  comparisons  of  the  experimental  data  with  corresponding  theoretical 
findings  can  be  made  through  referral  to  the  radial  dependence  shown  in  Fig.  4.  The  three 
indicated  situations  concerning  the  threshold  intensity  for  Xe(M)  radiation  are  the  following: 
(a)  The  experimentally  determined  threshold  PthCa)  s  3  x  10^^  W/cm^]  is  shown,  a  value 
which  agrees  closely  with  the  theoretical  magnitude  I^»  shown  in  Fig.  3  that  arises  with  the 
inclusion  of  the  coherent  electronic  motions  [6,9].  This  choice,  as  noted  above,  corresponds  to 
a  maximum  radius  of  r(a)  -  27  ^m.  03)  The  threshold  intensity  that  would  apply  [5],  ij  the 
mherent  interaction  were  absent,  is  l^W)  =  4  X  10^^  W/cm^.  This  case  corresponds  to  a 
relatively  small  radial  extent  of  T(fi)  -  2  ^m.  (7)  The  value  PjJjCt)  =  9.5  x  10^^  W/cm^]  is 
also  shown,  the  threshold  determined  by  matching  the  observed  radial  size  [r(7)  -  60  /xm]  of 
the  X-ray  emitting  region  to  the  computed  distribution. 

The  theoretical  radial  loci  for  the  former  two  threshold  intensities  (a  and  /3)  are  compared 
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to  the  corresponding  experimental  data  in  Figs.  1(d)  and  2(d).  Particularly  significant  is  the 
strong  dependence  experimentally  observed  at  z  -  1250  /xm  (locus  A)  in  Fig.  1(d)  and  the 
matching  feature  that  is  clearly  present  in  the  computed  locus  (B)  for  the  threshold 
incorporating  the  coherent  electronic  process.  In  contrast,  the  locus  (C)  in  Fig.  1(d),  which 
corresponds  to  the  threshold  Pj|j(j8)]  that  would  apply  without  the  coherent  process,  does  not 
exhibit  this  structure  and  differs  grossly  in  magnitude  and  in  form  from  both  the  experimental 
measurements  (A)  and  that  calculated  with  the  threshold  involving  the  coherent  coupling  (B). 
Finally,  as  seen  from  Fig.  4,  a  channel  radius  of  r(7)  —  60  /xm,  the  magnitude  matching  the 
measured  value,  could  only  result  if  the  threshold  intensity  Pth(7)]  were  reduced  to  ~  10^'* 
W/cm^.  However,  since  the  ponderomotive  potential  is  —  1  eV  at  this  intensity  and 
measurements  have  shown  that  electron  heating  plays  a  negligible  role  in  producing  the  observed 
X-ray  emission  [7],  there  exists  no  physical  basis  for  a  threshold  this  low  and  this  possibility  is 
ruled  out.  Therefore,  it  is  not  possible  to  reconcile  either  the  shapes  or  the  magnitudes  of  the 
calculated  radial  loci  with  the  experimental  measurements  without  the  enhanced  coupling 
characteristic  of  the  coherent  mechanism  [6,9]. 

The  results  illustrated  in  Fig.  1(d)  and  2(d)  show  that  the  calculated  radial  loci  for  the 
case  involving  the  coherent  motion  Ptii(o()  in  Fig.  4]  fall  consistently  below  the  measured  values 
by  a  factor  of  ~  2.  Since  the  field  occurring  in  the  region  exterior  to  a  dielectric  waveguide 
is  sensitive  to  the  detailed  structure  of  the  radial  index  profile  [30],  small  inaccuracies  in  the 
computation  of  the  electron  hydrodynamics  and  the  resulting  radial  electron  density  distributions 
could  cause  the  observed  differences.  In  the  current  formulation,  the  theory  [1]  calculates  the 
electron  density  by  assuming  the  instantaneous  establishment  of  a  force  balance  between  the 
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radially  outward  ponderomotive  force  and  the  oppositely  directed  electrostatic  restoring  force. 
Since  the  true  electronic  motions  are  not  instantaneous  [31],  this  computational  procedure  has 
the  proclivity  to  maximize  the  displaced  charge  and,  consequently,  to  overestimate  the  refractive 
confinement  of  the  channel.  The  outcome  is  the  collateral  tendency  to  underestimate  the  field 
intensity  outside  the  radial  zone  carrying  the  main  channeled  power.  The  effects  of  electron 
collisions  and  an  electron  temperature  Tg  >  0,  which  are  also  neglected,  would  produce  a 
similar  influence. 

The  present  experiment  is  a  particularly  discriminating  probe  of  the  dynamics  'of  the 
electron  density,  since  the  radial  extent  of  the  Xe(M)  emission  is  highly  sensitive  to  variations 
in  the  radial  intensity  profile  of  the  channeled  248  nm  radiation.  As  shown  in  Fig.  4,  there 
generally  exists  a  broad  second  maximum  (Ijjj)  in  the  intensity  profile.  When  Ijjj  —  Ijjj,  the 
situation  corresponding  to  the  case  Ith(a)  involving  the  coherent  interaction,  a  modest  change 
in  the  value  of  Ijjj  can  lead  to  relatively  large  changes  in  the  radius  of  the  emitting  zone. 
Specifically,  a  decrease  in  the  value  of  by  a  factor  of  only  —  3  leads  to  a  fourteen-fold 
reduction  in  the  radius  of  the  emitting  region;  the  transverse  dimension  of  the  emitting  zone 
collapses. 

It  has  been  possible  to  establish  experimental  conditions  under  which  this  radial  collapse 
can  be  controlled  and  Fig.  5(a)  illustrates  a  demonstration  of  this  phenomenon.  Since  the  plasma 
density  ng  will  slowly  pass  through  a  maximum  as  the  propagating  charmel  transverses  the 
spatially  inhomogeneous  hydrodynamically  produced  target,  and  since  the  critical  power  P^j. 
depends  inversely  on  ng  [1],  it  is  possible  to  produce  the  conditions  <  1,  ijq  ^  1  with 
<  Itjj(a),  and  %<  I  sequentially  along  the  path  of  propagation.  Under  these  circumstances, 
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the  observed  X-ray  emission  will  suddenly  undergo  radial  collapse  as  the  channel  forms,  exhibit 
a  narrow  extended  zone  of  propagation  because  Ijj^  <  Itji(a),  and  abruptly  resume  the  expanded 
profile  when  the  critical  condition  (tiq  ^  1)  for  confined  propagation  again  fails  to  be  ^tisfied. 
This  pattern  of  emission  conforms  to  the  observed  distribution  shown  in  Fig.  5(a).  Figure  5(b) 
illustrates  the  radial  profile  of  the  Xe(M)  emission  at  a  point  (z  =  956  fim)  in  the  narrow  region, 
the  characteristic  width  of  which  matches  the  estimated  spatial  resolution  limit  of  the 
measurement  Ajjj  —  A  —  30  /tm. 

For  the  sensitive  conditions  of  this  study,  the  difference  of  a  factor  of  —  2  betwden  the 
calculated  and  observed  radii  shown  in  Figs.  1(d)  and  2(d)  indicates  that  the  level  of 
approximation  used  to  calculate  the  electron  distribution  is  sufficiently  good  that  it  leads  to  a 
reasonable  assessment  of  the  intensity  profile  of  the  propagating  mode.  Although  further  study 
is  required  to  verify  these  hypotheses  concerning  the  electron  hydrodynamics,  we  tentatively 
attribute  the  difference  in  the  calculated  and  observed  loci  to  the  specific  approximations  used 
in  calculating  the  distribution  of  the  displaced  charge  that  is  dynamically  generated  by  the 
channeled  propagation. 

The  information  on  the  scaling  relationships  governing  the  cluster  excitation  and  the 
confined  propagation  revealed  in  these  studies  has  exceptionally  important  implications  for  their 
effective  cooperative  action.  The  results  show  how  a  judicious  design  can  harness  these  two 
very  different  nonlinear  processes  to  produce  harmoniously  conditions  ideal  for  X-ray 
amplification.  Specifically,  in  order  to  convert  a  substantial  quantity  of  energy  efficiently,  we 
require  (1)  that  the  range  of  intensities  characteristic  of  the  propagating  channel  conform  to  the 
span  of  intensities  specified  by  the  corresponding  allowed  zone  for  core-state  excitation  of  the 
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clusters  [4-6]  and  (2)  that  this  overlap  in  conditions  involving  the  propagating  energy  encompass 
a  large  spatial  volume.  The  simultaneous  fulfillment  of  these  two  requirements  constitutes  a 
specification  of  the  optimum  conditions  for  X-ray  amplification. 

In  order  to  illustrate  the  procedure  involved  in  establishing  this  optimum,  we  will  refer 
to  the  information  concerning  the  Xe(M)  emission  presented  in  Figs.  3  and  4.  It  is  known  from 
the  analysis  [1]  of  relativistic  and  charge-displacement  self-channeling  that  a  critical  power  P^j. 
and  a  normalized  radius  of  the  lowest  eigenmode  Pg^g  exist  which  enable  the  peak  intensity  in 
the  channel  Igjj  to  be  expressed  as 


‘dk 


jV «  ^  /  0 


(1) 


where  s  Pgjj/Pgj.  and  Pgjj  represents  the  power  in  the  channel  [32].  Furthermore,  since 
the  normalized  radius  of  the  eigenmode  [3]  is  nearly  constant  (pgjg  =1.7)  over  the  range  of 
principal  interest  (1.4  <  <  10),  Eq.  (1)  reduces  to  the  simple  expression 

7 .  -  7  X  10^*^  Wlcm^  (2) 

^2 

with  X  given  in  units  of  micrometers.  The  computed  value  of  Igj^  for  the  conditions 
corresponding  to  the  data  shown  in  Fig.  2  (X  =  0.248  pm,  =  1.36)  is  indicated  on  Fig.  4. 
Furthermore,  it  was  generally  found  that  for  conditions  of  interest  the  ratio  [1]  was 

restricted  to  the  narrow  range  0.5  <  ^  0-7- 

Since  we  seek  ideally  to  bound  the  upper  intensity  produced  in  the  channel  with  the 
maximum  of  the  allowed  zone  shown  in  Fig.  3,  we  desire  the  condition 

U  •  4.xW-  ® 

Moreover,  since  the  relevant  Imax^J)  ^ax^^l  ^  ^  readily  determined  [5]  for  any 
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shell(j)  for  any  atom  from  a  tunneling  picture  involving  known  atomic  parameters  [29],  Eq,  (3) 
can  be  simply  and  generally  evaluated.  For  a  given  atomic  shell,  this  results  in  a  relationship 
between  the  wavelength  (X)  and  the  atomic  number  of  the  atom. 

The  analysis  of  the  radial  dependence  of  the  observed  X-ray  emission  discussed  above 
showed  that  a  large  increase  (~  10^)  in  the  volume  is  obtained  if  the  condition 

<  I.  w 

holds,  the  actual  situation  illustrated  in  Fig.  4  for  the  Xe(M)  emission.  The  threshold  intensity 
lp>  can  be  evaluated  from  the  known  constraints  on  the  allowed  zone  [5,6,33],  a  profcedure 
which  involves  atomic  spectroscopic  data,  the  interatomic  spacing  r^  in  the  cluster,  and  the 
wavelength  X.  The  corresponding  value  of  Ijjj  is  determined  by  the  dynamics  of  the  propagation, 
a  calculation  [1-3]  which  incorporates  the  electron  density  n^,  the  wavelength  X,  and  the  incident 
power  Pq.  Although  highly  nonlinear,  fast  numerical  methods  can  readily  determine  the  optimal 
zone  of  the  multidimensional  problem  stated  by  Eqs.  (3)  and  (4)  subject  to  the  sole  constraint 
that  the  plasma  is  underdense,  a  condition  necessary  for  propagation  to  occur.  The  outcome, 
for  a  specific  atomic  shell  and  a  given  incident  power  P^,  is  a  specification  of  the  wavelength 
X,  the  plasma  density  n^,  and  the  cluster  structural  parameter  r^. 

Explicit  calculations  have  revealed  important  characteristics  of  this  procedure  for 
optimization.  We  consider  the  special  situation  involving  a  fixed  wavelength  (X  =  248  nm)  and 
a  power  Pq  ~  1  TW,  since  these  experimental  parameters  can  be  practically  produced  [19]  and 
data  for  these  conditions  exist.  The  two  principal  features  of  the  scaling  that  emerged  are  (1) 
the  existence  of  a  short  wavelength  limit  X^  that  can  be  amplified  under  the  optimized  conditions 
and  (2)  the  strong  tendency  for  the  optimum  conditions  to  occupy  a  small  zone  in  the  available 
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space  of  parameters.  Concerning  the  former,  a  limiting  value  of  X/X^  is  expected  to  exist 
arising  from  the  natural  physical  bounds  on  plasma  density  and  intensity  associated  with  the 
production  of  stable  channels  [16,18]  with  radiation  at  wavelength  X.  Exploration  of  this 
question  has  led  to  the  limit  X/X^^  —  1(P.  A  good  example  of  a  system  rq)resenting  this  limiting 
behavior  concerns  the  generation  of  U(M)  emission  [34-36]  in  the  4-6  teV  range  from  U 
clusters. 

Analysis  of  the  conditions  representing  optimization  following  the  procedure  outlined 
above,  which  concentrated  on  transitions  in  the  ~  0.7  keV  to  ~  5  keV  range,  led  to  thd  result 
shown  in  Fig.  6.  The  optimal  zone  is  circumscribed  by  the  small  rectangle  near  the  center  of 
the  stable  channeling  region  [16,18].  Points  A  and  B  correspond  to  optimization  of  Xe(M) 
emission  from  Xe  clusters  (r^j  =  4.4  A)  with  radiation  at  248  nm  having  Pq  =  1  TW  and 
=  1.5  TW,  respectively.  Consideration  of  the  variation  of  the  atomic  and  molecular  parameters 
for  other  systems  showed  that  the  resulting  optimum  condition  was  weakly  dependent  on  these 
variations  over  the  relatively  narrow  ranges  physically  possible.  Hence,  the  optimal  zone  is 
highly  localized  and  the  rectangular  region  in  Fig.  6  is  drawn  to  indicate  its  estimated  extent. 
The  chief  result  is  that  the  optimal  conditions  for  amplification  up  to  ~  5  keV  occupy  a 
relatively  small  zone  that  is  centrally  located  in  the  stable  region. 

The  conditions  corresponding  to  the  experimental  data  shown  in  Fig.  1  and  Fig,  2  are 
designated  in  Fig.  6  by  points  El  and  E2,  respectively,  both  of  which  are  at  locations 
significantly  removed  from  the  optimal  region.  However,  comparison  of  the  images  shown  in 
Figs.  1(a)  and  2(a)  clearly  shows  that  significantly  greater  emission  occurs  for  the  latter  (E2). 
This  result  is  fully  consistent  with  the  analysis  given  above,  since  as  one  travels  from  El  to  E2 
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in  Fig.  6,  the  movement  is  directly  toward  the  specified  optimal  zone. 

In  conclusion,  the  ability  to  controllably  apply  power  densities  comparable  to 
thermonuclear  levels  is  a  fundamental  requirement  for  achieving  the  amplification  of  X-rays. 
Two  new  nonlinear  phenomena  involving  (1)  the  multiphoton  production  of  X-rays  from  clusters 
and  (2)  channeled  propagation  are  being  united  to  produce  and  control  the  required  power 
density.  Experiments  showing  the  first  combined  expression  of  these  two  complex  nonlinear 
processes  have  provided  confirming  evidence  for  the  action  of  a  superstrong  coherent  multi¬ 
electron  intense-field  interaction  in  the  X-ray  generation  from  the  clusters  and  fumishSd  new 
detailed  information  on  the  dynamics  of  the  propagation.  A  crucial  result  is  the  demonstration 
of  the  exceptional  compatibility  of  the  mutual  scaling  of  these  phenomena  for  realizing  the 
conditions  necessary  for  efficient  amplification  in  the  multi-ldlovolt  spectral  region.  Most 
importantly,  the  knowledge  obtained  in  these  studies  shows  explicitly  how  the  optimum 
conditions  for  amplification  can  be  easily,  generally,  and  narrowly  specified  up  to  a  quantum 
energy  of  —  5  keV.  The  combined  use  of  these  nonlinear  processes  is  expected  to  lead  to  an 
advanced  generation  of  extraordinarily  bright  X-ray  sources  in  the  multi-kilovolt  region  having 
a  peak  brightness  of  -  lo31-lo33  (mrad)'^  (mm)'^  (0.1  %  BW)"^  a  level  sufficient  for 
the  performance  of  microholographic  assays  of  living  cells,  tissues,  and  organisms. 
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FIGURE  CAPTIONS 


Fig.  1:  Data  and  calculations  associated  with  images  of  Xe(M)  emission.  Stagnation 

pressure  65  psia  and  plenum  temperature  233  K.  (a)  Pinhole  camera  single¬ 
exposure  recording.  The  color  scale  is  defined  by  black  (zero),  red  through  violet 
ascending  intensity,  and  white  (maximum).  The  Rayleigh  range  (28.5  ftm)  is 
indicated,  (b)  Axial  (radius  r  =  0)  line-out  of  X-ray  image  shown  in  panel  (a). 
The  distance  5  =  200  /xm  separating  two  adjacent  peaks  in  the  longitudinal  (z) 
direction  is  indicated,  (c)  Computed  normalized  intensity  profile  of  channeled 
248  nm  radiation  corresponding  to  experimental  conditions  for  data  illustrated  in 
panel  (a).  =  9.4  x  10^^  W/cm^,  =  0.6  TW,  and  n^  =  6.5  x  10^° 

cm*^.  The  initial  condition  for  the  calculation  used  a  focal  diameter  of  3  fim. 
The  distance  6  s  214  fim  separating  two  adjacent  peaks  in  the  longitudinal  (z) 
direction  is  indicated,  a  value  that  is  close  to  the  experimental  magnitude  shown 
in  panel  (b).  (d)  Comparison  of  the  transverse  structures  of  the  observed  Xe(M) 
emission:  A  Observed  radial  locus  corresponding  to  data  shown  in  panel  (a).  B 
Computed  radial  locus  corresponding  to  Ith(«)  =  3  x  10^^  W/cm^.  C 
Computed  radial  locus  corresponding  to  I^hOS)  s  4  x  10^^  W/cm^.  See  text  for 
details. 

Fig.  2:  Data  and  calculations  associated  with  images  of  Xe(M)  emission.  Stagnation 

pressure  95  psia  and  plenum  temperature  233  K.  (a)  Pinhole  camera  single- 
exposure  recording.  The  color  scale  is  defined  by  black  (zero),  red  through  violet 
ascending  intensity,  and  white  (maximum).  The  Rayleigh  range  (28.5  /tm)  is 
indicated,  (b)  Axial  (radius  r  =  0)  line-out  of  X-ray  image  shown  in  panel  (a). 
The  inset  shows  the  relative  positions  of  the  0.5  mm  nozzle  and  the  barrier  with 
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the  ~  100  opening,  (c)  Computed  normalized  intensity  profile  of  channeled 
248  nm  radiation  corresponding  to  experimental  conditions  for  data  illustrated  in 
panel  (a).  =  9.4  x  10^^  W/cm^,  =  0.6  TW,  and  n^  =  9.6  x  10^0 

cm**’.  The  initial  condition  for  the  calculation  used  a  focal  diameter  of  3  /im. 
(d)  Comparison  of  the  transverse  structures  of  the  observed  Xe(M)  emission:  A 
Observed  radial  locus  corresponding  to  data  shown  in  panel  (a).  B  Computed 
radial  locus  corresponding  to  IthCa)  s  3  x  10^^  W/cm^.  C  Computed  radial 
locus  corresponding  to  IjjjCS)  =  4  x  10^^  W/cm^.  See  text  for  details. 

f 

Fig.  3:  Representation  of  the  allowed  zone  (shaded)  for  production  of  prompt  Xe  M-shell 

emission  from  Xe  clusters  as  a  function  of  intensity  of  irradiation  and  cluster  size 
for  excitation  with  a  wavelength  of  248  nm.  The  upper  boundary  Iniax^^  “ 
established  by  the  intensity  required  to  fully  ionize  the  Xe  N-shell  as  found  with 
a  tunneling  picture  [Ref.  (27)].  The  lower  boundaries  [Iq  and  Ijj(248  nm)]  are 
determined  by  a  classical  formulation  of  the  work  done  by  an  intense  external 
wave  on  an  initially  ionized  free  electron  in  the  cluster  Refs.  [4-6].  The  effect 
of  coherent  electron  motion  is  included  [Ref.  (6)].  The  point  jS'  is  specified  by 
n^/  s  3  and  I^/  s 3.8 X 10^^  W/cm^. 

Fig.  4:  Transverse  distribution  of  channeled  248  nm  intensity  distribution  shown  in  Fig. 

2(c)  for  longitudinal  position  z  =  1050  ftm.  The  three  levels  (i  =  a,  j8,  7)  for 
the  threshold  intensities  [I^(i)]  are  indicated.  Corresponding  radial  parameters 
r(a)  s  27  /xm,  if/S)  =  2  /xm,  and  r(7)  =  60  are  also  shown.  See  text  for 
discussion.  The  peak  intensity  in  the  channel  Ich.  defined  in  Eq.  (2),  and  the 
intensity  representing  the  broad  second  maximum  of  the  distribution  are 
designated.  The  characteristic  intensities  Imax^  ^/3'»  transferred  from  Fig. 
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3,  are  also  indicated.  The  magnitudes  are  such  that  lyiCa)  =  The  intensity 
rq)resenting  the  initial  condition  for  the  calculated  distribution  is  Iq  =  9.4  x 
lO^^W/cm^. 

Fig.  5:  (a)  Image  of  Xe(M)  emission  with  pinhole  camera  having  a  25  urn  pinhole  for  a 

1:2  (Xe:He)  mixture  with  a  stagnation  pressure  of  90  psia  and  a  plenum 
temp^mure  of  293  K.  Unlike  the  images  shown  in  Figs.  1(a)  and  2(a),  the 
nozzle  had  a  diameter  of  0.8  mm  and  the  barrier  illustrated  in  Fig.  2(b)  was 
removed.  The  arrow  specifies  the^  longitudinal  position  (z  =  956  /tm)  of  the 
transverse  distr^tion  shown  in  panel  (b).  The  color  scale  is  defined  by  black 
(zero),  red  through  violet  ascending  intensity,  and  white  (maximum),  (b) 
Transverse  Xe(M)  intensity  profile  at  z  =  956  /xm  with  a  measured  width  s 
30  /xm,  a  value  in  agreement  with  the  estimated  resolution  limit  A  ~  30  ;xm. 

Fig.  6:  Stability  map  for  the  propagation  in  initially  homogeneous  plasmas  of  pulses 

having  perturbed  Gaussian  incident  transverse  intensity  distributions  as  computed 
in  Ref.  (16).  =  Ro^p/c  with  Rq  =  radius  of  the  incident  Gaussian  profile, 

Qp  the  plasma  frequency,  and  c  the  speed  of  light.  77^  =  Po^Pcr’  incident 
peak  power  normalized  by  the  critical  power  defined  in  Ref.  1.  The  optimal 
zone  is  designated  by  the  small  rectangle.  Points  A  and  B  correspond  to 
optimization  of  Xe(M)  emission  with  248  nm  radiation  with  Pq  =  1  TW  and  Pq 
=  1.5  TW,  respectively.  Point  El  (tjq  =  1.32,  =  6.84)  corresponds  to  the 

experimental  data  shown  in  Fig.  1.  Point  E2  (tj^  =  1.94,  p^  =  8.31) 
corresponds  to  the  experimental  data  shown  in  Fig.  2.  The  eigenmode  curve  is 
shown  in  addition  to  the  areas  representing  stable  and  unstable  (shaded) 
propagation. 
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Abstract 


The  spatial  control  of  the  energy  deposited  for  excitation  of  an  x-ray  amplifier  plays  an  important 
role  in  the  fundamental  scaling  relationship  between  the  required  energy,  the  gain  and  the  wavelength. 
New  results  concerning  the  ability  to  establish  confined  modes  of  propagation  of  short  pulse 
radiation  of  sufficiently  high  intensity  in  plasmas  lead  to  a  sharply  reduced  need  for  the  total  energy 
deposited,  since  the  concentration  of  deposited  power  can  be  very  efficiently  organized. 


I.  Discussion  of  Research 

Recent  theoretical* and  experimental^  studies  have  led  to  a  fundamental  development  concerning 
the  generation  of  high-brightness  x-ray  eources.  These  results  affect  our  ability  to  controllablv 
apply  very  high  power  densities  in  materials,  the  basic  issue  for  the  creation  of  bright  and  efficient 
sources  of  radiation  in  the  x-ray  range.  The  main  significance  of  this  work  is  the  establishment 
of  the  scaling  law  concerning  the  energy  requirements  for  x-ray  amplification  in  the  kilovolt  range 
shown  in  Fig.  1.  Importantly,  the  parameters  represented  in  Fig.  1.  and  which  define  the  relationship 
presented,  are  based  on  both  theoretical  and  experimental  information. 

The  critical  governing  issue,  which  determines  the  scaling  relationship  between  the  required 
excitation  energy  (E)  and  the  amplifier  gain  (G)  of  x-ray  lasers,  is  the  spatial  control  of  the 
deposited  energy.  The  information  presented  in  Fig.  2  shows  that  optimizing  the  gain  (G)  per 
unit  energy  (E)  calls  for  the  guided  mode  of  propagation  in  order  to  optimally  control  the  deposition 
of  the  energy.'*  Overall,  in  comparison  to  traditional  forms  of  excitation,  for  a  fixed  x-ray  energy 
output  {Ey)  and  wavelength  (Xj^),  a  reduction  of  several  orders  of  magnitude  in  the  necessary 
energy  (E)  results,  as  shown  in  Fig.  1,  if  this  form  of  confined  (channeled)  propagation  can  be 
achieved.  Therefore,  if  this  scaling  holds,  a  relatively  small  and  useful  laboratorv-scale  technology 
becomes  feasible. 

Recent  experiments,^  which  are  supported  by  carefully  developed  theoretical  analysis,*"^  have 
demonstrated  the  basic  physics  of  a  new  form  propagation  exactly  of  the  type  necessary  for  the 
implementation  of  x-ray  lasers  of  a  fundamentally  new  regime  of  electromagnetic  propagation  is 
expected  to  arise  in  plasmas  for  short-pulse  radiation  at  sufficiently  high  intensity.  Dynamical 
calculations  of  the  propagation  in  plasmas,  incorporating  both  relativistic*-®  and  charge-displacement 
mechanisms, indicate  that  the  combined  action  of  these  processes  can  lead  to  a  new  stable 
form  of  spatially  channeled  propagation.  Specifically,  these  experimental  studies  which  have 
examined  a  new  relativistic  regime  of  high-intensity  short-pulse  propagation  in  plasmas,  present 
evidence  for  the  information  of  such  a  stable  mode  of  spatially  confined  (channeled)  propagation. 
For  an  electron  density  of  -  1.35  x  10**  cm'*  and  a  power  of  -  3  x  10**  W,  the  results  indicate 
a  channel  radius  <  1  pm  and  a  peak  intensity  -  10**  W/cm*.  Comparison  of  these  findings  with 
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Fig.  1:  Scaling  relationship  between  required  excitation  energy  (E)  and  quantum  energy 

characteristic  of  the  amplifier.  Parameters:  total  gain  exponent  G  =  100,  energy  efficiency 
T\^  =  10'*,  channel  diameter  6  =  3  pm,  x-ray  (hw^)  cross  section  for  stimulated  emission 
Oy,  x-ray  cross  section  for  stimulated  emission  for  radiatively-broadened  transition  OxR- 
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Fig.  2;  Spatial  distribution  of  energy  of  excitation  (E)  for  an  x-ray  (ft^x)  afriplifier.  Parameters 
are  the  same  as  in  Fig.  1  with  X  the  wavelength  of  the  excitation  energy,  assumed 
longitudinally  delivered,  and  with  p.  p*.  and  Pg  representing  the  particle,  inversion  and 
electron  densities,  respectively.  The  analysis  shows  that  optimization  of  G/E  requires  a 
guided  mode  of  propagation  so  that  high  concentrations  of  power  can  be  organized  into 
high-aspect-ratio  spatial  volumes. 
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a  dynamical  theory^  yield  close  agreement  for  both  the  longitudinal  structure  and  the  radial 
extent  of  the  propagation  observed.  These  results  represent  a  profound  change  in  the  field  of-  x- 
ray  laser  research  because  they  alter  drastically  the  fundamental  scaling  relationships  among  the 
relevant  physical  variables. 

The  implications  of  this  development  for  general  applications  to  x-ray  imaging  and  the  micro¬ 
characterization  of  condensed  matter  are  extremely  important  and  propitious.  In  terms  of  the  x- 
ray  source,  they  are  (1)  that  a  properly  controlled  energy  deposition  rate,  sufficient  for  the 
production  of  stimulated  x— ray  sources  up  to  a  few  kilovolts  in  quantum  energy,  can  now  be 
achieved  with  an  excitation  energy  of  ~  1  J,  (2)  that  an  x-ray  output  energy  of  ~  1  mJ  per  pulse 
is  achievable  with  laboratory-scale  technology,  and  (3)  that  an  x-ray  beam  diameter  (-2-3  pm) 
arises  as  a  natural  consequence  of  the  physics.  These  parameters  represent  an  exceptionally  high 
peak  brightness  figure  that  permits  a  new  and  completely  unexplored  range  of  physical  measurements 
to  be  made.  Indeed,  a  high— brightness  source  of  this  nature  is  ideal  for  the  microimaging  of 
condensed  matter.  In  particular,  an  x-ray  source  with  these  parameters  is  perfectly  matched  to 
the  requirements  for  holographic  imaging  of  biological  material in  terms  of  an  its  relevant 
properties,  specifically,  wavelength  (10  -  40  A),  pulse  energy  (~  i  mJ),  pulse  length  (~  10'“^  s), 
beam  diameter  (-2-3  pm),  and  divergence  (-  1  mrad). 


II.  Conclusion 


Recent  experimental  and  theoretical  results  on  electromagnetic  propagation  at  high  intensities 
in  plasmas  lead  to  very  favorable  scaling  relationships  for  high-brightness  x-ray  amplifiers.  Such 
x-ray  sources  will  have  important  applications  in  the  holographic  imaging  of  biological  materials. 
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